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RESULTS OF THE FIRST CONGRESS 


OF THE INTERNATIONAL FEDERATION 


OF AUTOMATIC CONTROL (IFAC) 


Translated from Avtomatika i Telemekhanika, Vol. 21, No. 9, pp. 1229-1239, 


September, 1960 


The sessions of the IFAC First Congress, which was 
prepared and organized by the National Committee of 
the Soviet Union for Automatic Control, were held 
from June 27 to July 2 on the premises of the Lomono- 
sov Moscow State University. The Congress was a very 
important event in the history of automatic-control] de- 
velopment. Its import can be judged by the following 
figures. The number of participants at the Congress in- 
cluded 1190 delegates from 29 countries, and almost as 
many guests, During the sessions of the 21 Congress 
sections, 270 reports were heard and discussed; over 400 
statements were pronounced in the discussions of these 
reports. 

The Congress was opened by the address of the 
IFAC president, Prof. A. M. Letov, on June 27 at the 
first plenum, The IFAC past-president, G. Chestnut, 
welcomed the Congress delegates in his speech at this 
session. The first Vice-President of the Council of 
Ministers of the USSR,Comrade A, N.Kosygin addressed 
the Congress delegates _in the name of the Soviet 
Government, The Congress delegates were greeted by 
the President of the Academy of Sciences USSR, Acade- 
mician A, N, Nesmeyanov, Prof. E, Gerecke (Switzer- 
land), Prof. R. Oldenburger (USA), and Comrade Wu Ju- 
yen (Chinese People's Republic). 

In his closing address at the first plenum, the pres- 
ident of the National Committee of the Soviet Union 
for Automatic Control, Academician V. A, Trapezni- 
kov, read the report “Automation and Mankind." 

The subsequent work of the Congress took place 
during the sessions of different sections. The problems 
in automatic-control theory were considered in nine sec- 
tions; five sections were devoted to problems in provid- 
ing the technological basis of automation; the work in 
six sections concerned the problems in applying auto- 
mation theory and technology in automating industrial 
processes; one section was devoted to general problems 
pertaining to automatic control as a scientificotechni- 
cal branch of science. 

The closing plenum was held on June 2. During 
this session, the following persons surveyed the scientif- 
icotechnical problems discussed at the Congress; the 
president of the IFAC Scientificotechnical Committee 
on Theory, Academician B, N. Petrov (USSR), the pre- 
sident of the IFAC Scientificotechnical Committee on 


Automation Technology, Prof. G. Boromiss (Hungary), 
and the president of the IFAC Scientificotechnical Com- 
mittee on Practical Application, Prof. D. Mosley (USA). 
The president of the IF AC Consulting Committee, Prof. 
D, Eckman (USA), made an announcement concerning 
the activities of this committee. The president of the 
IF AC Scientificotechnical Committee on Terminology, 
Prof. E. Gerecke (Switzerland), dwelled on the problems 
in the work of this committee, The Congress was closed 
by the final address of the IFAC president, Prof. A. M. 
Letov. 


SURVEY OF THE SCIENTIFICOTECHNICAL 
PROBLEMS CONSIDERED AT THE 
CONGRESS* 


Automatic-Control Theory 





The theory of automatic control was treated in 
139 reports, The reports themselves as well as the en- 
suing discussions were a vivid demonstration of the high 
level attained in the development of the theory of auto- 
matic control, which is based on a widespread use of rig- 
orous mathematical methods and is concerned with the 
solution of basic problems in the analysis and synthesis 
of modern automation systems, 

Theoretical problems were discussed in the work of 
nine sections, 

Twenty-three reports were discussed in the Section 
dealing with the Theory of Continuous Linear Systems. 
Reports devoted to problems in the invariance theory 
and the theory of multiple feedback and autonomous 
control systems created the greatest interest and the 
liveliest discussions, The discussion of these problems 
showed that, at the present time, the theory provides 
the basic principles for the design of automation sys- 


* This survey was compiled from the above-mentioned 
reviewing reports by the three presidents of the IFAC 
Scientificotechnical Committees: Academician B. N. 
Petrov (USSR, Committee on Theory), Prof. G. Boromiss 
(Hungary, Committee on Technology), and Prof. D. 
Mosley (USA, Committee on Practical Application), 
This survey does not review the work of the Section for 
General Problems in Automatic Control (terminology, 
training of specialists, etc.). 
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tems which respond little to undesirable disturbances and 
accurately follow the control inputs. 

In the work of this Section, much attention was 
paid to the analysis of multiple feedback control sys- 
tems and general methods for their synthesis, methods 
for improving the dynamic properties of control systems 
with one and several controlled parameters, and also 
problems in the theory of systems with lag elements. 

A number of reports were devoted to improved and 
generalized methods for determining the stability and 
the performance of automatic-control systems. 

Fourteen reports were discussed in the work of the 
Section dealing with the Theory of Continuous Nonline- 
ar Systems, Exact as well as approximate methods for 
the investigation of such systems were surveyed. 

A number of reports were concerned with the devel- 
opment of the describing-function method. It was shown 
that this method, which gained widespread use in engi- 
neering practice, can be extended, in matrix representa- 
tion, to the case of multiperiod oscillations, subharmon- 
ic and superharmonic resonance phenomena, oscillation 
synchronization, and parametric amplification. An ap- 
proach to the solution of some problems in the synthesis 
of nonlinear systems subject to random disturbances, bas- 
ed on a combination of the describing-function and the 
statistical-linearization methods, was indicated. 

It was also shown that, by using different variants of 
the describing-function method, it is possible to study 
the effect of the limited power of feec sources in control 
systems and to find conditions for suppressing self-oscil- 
lations by external periodic inputs, and to determine the 
conditions for suppressing self- oscillations in systems with 
backlash, One of the reports presented a variant of the 
describing-function method where Nichols loci are used. 

Two points of view were expressed in the discussion 
on the application of the describing-function method. 
One of them states that, regardless of the extreme effec - 
tiveness of the describing-function method, it should be 
used with a certain caution, since the fact that no self- 
oscillations are detected by means of this method does 
not confirm the absence of any complex periodic mo- 
tions in a system. According to the other point of view, 
the simple and efficient describing-function method 
yields accurate results in practically all problems en- 
countered in practice, and can be developed so that it 
yields subharmonic solutions. 

Besides the above problems, the following problems 
were considered in this Section: effect of programmed 
variation of the regulator gain on the control perfor- 
mance, problems in calculating the static characteristics 
of complex nonlinear systems, and a new method for cal- 
culating the integral quadratic criterion of nonlinear sys- 
tems. 

In the Section dealing with the Theory of Discrete 
Systems, twenty-five reports were discussed; problems in 
the theory of relay (off-on), sampled-data, and digital 


* 


automation systems were considered, Along with prob- 
lems in the analysis and synthesis of discrete systems in- 
tended for the solution of ordinary control problems, 
much attention was paid to problems in designing opti- 
mum and self-adjusting discrete systems. 

The reports devoted to the theory of relay systems 
presented exact and approximate methods for the analy- 
sis of subharmonic oscillations,as well as methods for the 
design of relay systems with optimum quickness of re- 
sponse which are subject to random disturbances, 

Reports relating to the theory of sampled-data sys- 
tems were the most numerous in this Section, A survey 
of mathematical methods for the investigation of such 
systems was presented; among the problems discussed 
were those encountered in the design of multimode and 
multiple feedback sampled-data systems, problems in 
the design of sampled-data systems with variable para- 
meters, methods for statistical calculation of sampled- 
data systems, and methods for iavestigating nonlinear 
sampled-data systems, The problems in the design of 
optimum and self-adjusting sampled-data systems whose 
structure comprises feedback with respect to the perfor- 
mance criterion were also considered, A number of re- 
ports were devoted to the theory of relay sampled-data 
(digital) systems. The simplest systems of this type were 
analyzed by using the phase-plane method and by intro- 
ducing lags with random variation. The recently devel- 
oped branches of the theory of digital automatic systems, 
including the determination of the effect of quantization 
with respect to a certain level, methods for determining 
periodic modes, the synthesis of systems with lag com- 
pensation, and the synthesis of optimum systems, were 
presented, The liveliest discussion developed with re- 
gard to problems in the application of digital computers 
in sampled-data control systems and the application of 
the phase-space and operator methods in investigating 
sampled-data systems. The necessity for developing 
methods for the analysis and synthesis of nonlinear sam- 
pled-data systems and sampled-data systems subject to 
random disturbances was emphasized in the discussion, 

Eighteen reports were discussed in the work of the 
Section dealing with Stochastic Problems. One of the 
main problems in this field is that of the statistical anal- 
ysis and synthesis of nonlinear systems. Three different 
trends in work on this group of statistical automation prob- 
lems became apparent through the reports and the ensuing 
discussions, The first trend pertains to the category of 
problems which include the development and generaliza- 
tion of the statistical-linearization method and its appli- 
cation,as well as the formulation of an exact theory for 4 
given class of nonlinear systems, which would be based 
on the theory of Markov's random processes, and the de- 
velopment of methods for the synthesis of systems which 
can be reduced to linear systems. During the discussions 
great interest on the part of scientists from different coun- 
tries was shown in the statistical-linearization method, 
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which is conceptually close to the harmonic linearization 
method and the deterministic theory of nonlinear sys- 
tems. In addition to what was said in the course of the 
discussions regarding the application of the Markov the- 
ory of random processes in the analysis of automatic 
systems, Soviet scientists reported a number of interest- 
ing results which were not known outside of our country 
before. In the discussion concerning problems in the 
statistical theory of nonlinear systems capable of reduc- 
tion to linear systems, a number of interesting results, 
promising successful application in the design of self- 
adjusting and self-optimizing systems, were announced, 

The second trend pertains to problems in the theory 
of statistical synthesis of linear systems, The discussion 
revealed the appearance in this field of two new com- 
petitive approaches in the USA which have been well 
developed in the papers of American as well as Soviet 
scientists, These two approaches are the application of 
the dynamic programming method for the determination 
of the optimum program for the controlling section of a 
system and the development of methods for the direct 
determination of differential or difference equations of 
optimum systems. 

The third trend pertains to problems in finding the 
algorithms of optimum systems, It was demonstrated in 
the discussion that the majority of problems in detecting 
and processing information in the presence of noise and 
random disturbances can be considered from a unique 
point of view and solved by a unique theoretical meth- 
od, 

In addition to the reports and the discussion con- 
cerning the above three basic trends, the work in this 
Section also included the discussion of problems in ex- 
perimental determination of characteristics of control 
systems and systems to be controlled,as well as of statis- 
tical characteristics of random processes, problems in 
statistical calculation of linear systems, problems in de- 
termining errors in the conversion of continuous infor- 
mation into discrete information, as well as problems in 
applying the information-theory methods to the deter- 
mination of the transmission capacity of automatic sys- 
tems and their components, 

Twelve reports were discussed in the Section deal- 
ing with the Theory of Optimum Systems; general meth- 
ods for investigating optimization problems,as well as 
the resvlts obtained in solving individual problems, were 
considered, The general methods for the ifivestigation 
of optimum systems are mainly connected with the de- 
velopment of methods for the solution of nonclassical 
variation problems, The report dealing with these prob- 
lems was devoted to the maximum principle and the 


theory of optimum automatic-control systems and aroused 


great interest and a lively discussion. The participants 
in the discussion remarked that this principle provides a 
powerful tool for the investigation of nonlinear optimiz- 
ation systems, while the restrictions imposed on the con- 


trol functions are of a general nature; at present, it is 
one of the most general principles, thus making it possi- 
ble to solve a large number of problems by means of it. 

Much attention was paid to problems met in the syn- 
thesis of optimum systems. The inverse time methods 
and methods for construction of isosurfaces in the system 
phase space, both of which are based on the application 
of the maximum principle, were discussed, Practical 
examples of system synthesis performed by the above- 
mentioned methods were given, 

The connection between the maximum principle 
and the dynamic programming method was analyzed 
during the discussion, It was noted that the dynamic 
programming method, when applied to optimum con- 
trol problems, yields functional equations which are 
closely related to the maximum principle equations. 
However, it was indicated that the application of the 
dynamic programming method involves considerable 
calculation difficulties. 

Great interest was exhibited in the formulation and 
solution of some problems in optimum control of sys- 
tems with distributed parameters (which are described 
by partial differential equations),as well as of general 
and particular problems in the theory of optimum linear- 
control systems. 

Twenty reports were presented to the Section deal- 
ing with the Theory of Self- Adjusting Systems. The 
items discussed in this Section were biological and me- 
chanical systems capable of adaptation to changing con- 
ditions,with the aim of securing such an interaction be- 
tween the system to be controlled and the surrounding 
medium, for which, in a given sense, the optimum state 
of the system to be controlled is secured, Systems where 
adaptation is secured by automatic scanning as well as 
by changing the system state with respect to an open- 
loop or closed-loop cycle (without automatic scanning) 
were considered, 

Problems in the following branches of the theory of 
automatic control were considered; the statistical the- 
ory of automatic optimization systems, dynamics of 
self-adjusting systems, bioelectric control, the study of 
characteristics of physiological control systems, and the 
theory of self-adjusting autopilots. 

Interesting reports were presented on the theory and 
application of computers for the approximate determin- 
ation of characteristics of systems subject to automatic 
and self-adjusting control, 

The reports concerning problems in biological cy- 
bernetics brought about an especially lively discussion, 
which showed that Soviet, Yugoslav, and American sci- 
entists contributed much to the development of princi- 
ples for the design of systems in many ways similar to 
living organisms, Six reports were devoted to the results 
obtained in studying biological self-adapting systems. 
The results obtained by a group of scientists in studying 
problems of bioelectric control attracted deserved atten- 
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tion; these results make it possible to use biocurrents in 
a living organism for the control of artificial organs sup- 
planting the lost ones. 

At present, the methods and procedures used in au- 
tomatic-control theory are applied in the investigation 
of biological systems to an ever-increasing extent. This 
trend in cybernetics, which is beginning to yield inter- 
esting results, was reflected in the reports and in the en- 
suing discussions, It should be noted that the discussion 
in the work of this Section was perhaps the liveliest, in- 
dicating the especially great interest of the scientific 
world in this important and promising direction of the 
automatic-control theory. 

It was indicated in a number of reports and in the 
following discussions that it was necessary to investigate 
and develop self-adjusting systems capable of operating 
under conditions of rapid drift of characteristics and pa- 
rameters of the systems to be controlled. These prob- 
lems are closely connected with problems in the stabili- 
ty of the self- adjustment process, and, therefore, the im- 
portance of exact mathematical formulation of stability 
conditions in self-adjusting systems was repeatedly em- 
phasized in the discussions, 

Seven reports were discussed in the Section dealing 
with the Theory of Structures and Signal Flow Diagrams. 
A number of important problems in the theory of relay 
devices and the theory of signal flow diagrams were dis- 
cussed; these problems included: development of meth- 
ods for the synthesis of relay devices and problems in the 
design of relay structures with assigned reliability; prob- 
lems in the application of a new logic language for the 
description of discrete time sequences; problems in the 
mechanization of the analysis and synthesis of relay de- 
vices; problems in noise stability and the efficiency of 
information transfer in remote measurement systems. 
However, due to the fact that no reports from a number 
of eminent scientists from abroad working in this field, 
scientists from the USA in particular, were submitted to 
the Section, many other important problems in the mod- 
ern theory of structures and the theory of signal flow dia- 
grams were not discussed. This gap was partially filled 
by scientists from the Soviet Union and abroad, who took 
part in the discussion. 

Considering that the structural theory of relay de- 
vices is one of the important branches of the modern au- 
tomatic-control theory, it was planned to hold in the 
nearest future a special symposium on the theory of re- 
lay devices and related scientific fields. 

The discussion in the Section dealing with Special 
Mathematical Problems comprised thirteen reports per- 
taining to three special mathematical problems in au- 
tomatic-control theory. 

The first problem pertains to the determination 
and investigation of steady-state conditions in automa- 
tic control systems which are described by differential 
equations with continuous right-hand sides. 
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A number of important aspects of this problem 
were treated in the reports; a more rigorous and exact 
mathematical definition of the steady-state conditions 
concept was given; the results obtained in investigating 
the dependence of the periodic solutions of a system of. 
differential equations on its parameters were described. 

The second problem pertains to the investigation 
of nonlinear systems with discontinuous characteristics. 
With regard to these systems, problems in determining 
steady-state conditions were discussed. 

The third problem pertains to the stability of au- 
tomatic-control systems. The main problems treated 
in the reports,and in the lively discussion they gave 
rise to,pertain to the further development of methods 
of determining stability by using the second Lyapunov 
method. 

Eleven reports were discussed in the Section deal- 
ing with Simulation and Experimental Investigation 
Methods. The application of simulation methods, where 
not only analog, but also digital and analog-digital 
(combined) devices are used for the investigation and 
design of complex control systems, is a separate and 
important problem. 

The solution of this group of problems will help in 
mechanizing certain human mental processes and con- 
siderably increase their efficiency. By using simulation 
methods, the experimental techniques are rendered 
more accurate and prerequisites for a more efficient 
use of equipment for the testing of actual systems are 
secured. The development of equipment and methods 
for digital simulation was widely reflected in the re- 
ports and the ensuing discussion; with regard to pro- 
grammed-control systems, new methods for the repro- 
duction of second-order curves by means of digital 
computers were proposed. The application of combin- 
ed digital-analog computers (in particular, for the re- 
production of nuclear reactor kinetics) was also consid- 
ered. 

It was noted that, besides. improving the accuracy, 
the application of digital-simulation methods consider- 
ably increases the maximum allowable range in which 
the variables change. 

A method which was proposed for the solution of 
partial differential equations by means of electronic 
models aroused great interest; this method is based on 
the connection between the equations to be solved and 
the processes which appear in practical circuits subject 
to random disturbances. 

One of the reports mentioned the development of 
an original simulating device which operates with high- 
frequency alternating current and which consists of pas- 
sive circuits combined with follow-up multiplier sys- 
tems; this device has a relatively high accuracy and is 
reversible, which considerably enlarges its scope. 


It was noted in the discussion that the new ideas 
in the application of electronic models for the optimiz- 
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ation of the performance of systems used for the auto- 
matic control of industrial processes offer great oppor- 
tunities for the design of optimum systems, as well as 
for obtaining more accurate information on the sys- 
tems to be controlled. 


Automation and Remote- 





Control Technology 








Fifty-five of the reports submitted to the Congress 
were devoted to problems in construction, investiga- 
tion, and design of automation and remote-control 
equipment. These reports were discussed at the ses- 
sions of five Sections. 

Nine reports were discussed in the Section dealing 
with Electric and Magnetic Components of Control Sys- 
tems. Problems in the theory of magnetic amplifiers 
(three reports) were most thoroughly treated. A de- 
tailed theoretical analysis of the effect of dynamic 
properties of magnetic cores in the presence of eddy 
currents and magnetic viscosity on the characteristics 
of the basic magnetic amplifier types was presented 
for the first time in one of these reports. These results 
make it possible to determine the theoretical charac- 
teristics of different amplifier types and to find ways 
for their improvement. 

Moreover, among the problems considered in the 
work of this Section were problems in the design of 
logic elements and their use in control circuits; prob- 
lems pertaining to measurement methods; and prob- 
lems in the design and investigation of some special 
electric automation devices (photoresistors and quick- 
action step motors). 

In the report devoted to logic components, atten- 
tion was drawn to the fact that a system of logic ele- 
ments can be without a special memory component, 
since, in the design of circuits requiring the application 
of such a component, one can use standard logic com- 
ponents and circuit feedbacks. 

Ten reports were discussed in the Section dealing 
with Electrical Computers and Simulating Devices and 
Programming Machines and Components. According to 
their character, the problems considered in these reports 
can be divided into three groups. 

The first group was dealt with in reports devoted to 
problems in the design of new computers. The follow- 
ing items were considered in these reports; the design of 
a new computer for algebraic functions of a complex 
variable; the design of a synthesizer in the shape of a 
nonuniform electronic lag line; the design of a discrete- 
ly acting computer of the electric-grid type intended 
for the solution of unsteady-state field equations. Also 
among the items discussed was a system for the control 
of industrial processes containing a correcting computer 
having a remote measurement selector for centralized 
control. (A ternary code is used in the selector of this 
computer.) Finally, this group included a report on the 


application of automatic optimizers for automatic syn- 
thesis and the solution of variation problems. 

The reports of the second group related to comput- 
ers for the automation of industrial processes. The 
SKZhAT electronic digital computer, designed for the 
control of such processes as chemical synthesis, distilla- 
tion, and blast-furnace processes, belongs to this class of 
computers. 

The reports of the third group treated problems in 
the improvement of computer components. There was a 
communication on the automatic zero stabilization in 
operating amplifiers, for which a number of circuits for 
reducing the zero drift were proposed and a comparative 
evaluation of different circuits was given. The problem 
of the design of nonlinear computer units, based on semi- 
conductor resistors made of silicon carbide, was discuss- 
ed. A report devoted to the logical synthesis of func- 
tional generators was read. 

The reports and the ensuing discussion showed the 
great importance of computer techniques in the solution 
of modern problems in automatic control . 

Eight reports were discussed in the Section dealing 
with Data Transmitters, Components, and Automation 
and Remote-Control Systems for Industrial Processes. 
The main trends in the development of remote control 
and measurement techniques were discussed in the re- 
ports presented to this Section. In the field of remote 
measurement, the main development trend consists in 
the application of contactless digital and code-pulse 
systems; in the field of remote control of a large number 
of distributed systems, this trend can also be reduced to 
the application of contactless equipment. In the work 
of this Section, much attention was paid to problems in 
the operational accuracy and reliability of remote mea- 
surement and control systems, Among the items discussed 
were performance criteria for remote control systems 
based on different characteristics, such as the specific 
signal content index, the information transfer quality in- 
dex, the information channel utilization factor, and the 
technical efficiency. 

Ten reports were discussed in the Section dealing 
with Pneumatic Devices in Automation and Computer 
Techniques. Problems in the design of the following 
pneumatic and hydraulic automation components were 
considered; pneumatic throttles of the nozzle-flapper 
type which control devices having jet pipes, servovalves 
with discharge feedback, and long lines which connect 
the control valve to the drive. Moreover, there was a 
report on the general principles of the theory and design 
of pneumatic automation elements. This generalized 
report aroused great interest and a lively discussion. 


Among the reports presented to this Section, the 
most interesting were those devoted to the development 
of pneumatic equipment designed for mathematical and 
logic operations, It was shown that pneumatic devices 
can be used successfully for the realization of analogous 
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simulating devices and that the most diverse relay cir- 
cuits for algebraic and time logic functions can be con- 
structed, 

A report containing an analysis of the dynamics of 
the copying machine hydraulic drive and a report devot- 
ed to the investigation of automatic-control systems for 
complex hydrodynamic machines aroused great interest. 
A report considering problems of the dynamic suspension 
of machine parts in a viscous liquid brought about a live- 
ly discussion. 

Eighteen reports, divisible into three groups, were 
discussed in the Section dealing with Automatic-Control 
Equipment and Devices. 


Sensing elements were considered in the reports of 
the first group. The results obtained in the development 
of a thermomagnetic compensation gas analyzer to be 
used for the accurate measurement of high oxygen con- 
centrations were discussed. The application of radioac- 
tive isotopes in automatic-control systems for mining 
extraction and sinking machines was demonstrated. 

The reports of the second group were devoted to 
improved measurement methods, 

Methods for the automation of sampling and chem- 
ical analysis and for the automation of measurements 
in plotting the dynamic characteristics of nuclear reac- 
tors were considered, A report presenting a method for 
determination of the hunting angle of the earth's satel- 
lites by means of a gyroscope and by means of devices 
using nuclear radiation was discussed. 

The third group of reports was devoted to the results 
achieved in the design of individual intrument parts. 
These achievements include, for instance, the produc- 
tion of components with a higher temperature stability 
and a method for determining optimum vibrations,which 
exclude the effect of Coulomb friction. 


Application of Theory and 
Technology for the Automation 
of Industrial Processes 











Problems in practical application of the automation 
theory and technology were treated in sixty-seven re- 
ports. These reports were discussed at the sessions of 
six Sections, 


Ten reports were discussed in the Section concerned 
with Automation of Machine Construction. In the work 
of this Section, most attention was paid to problems in 
the ¢ssign and practical realization of automatic-con- 
trol systems for different machines. These problems 
comprised the selection of the simplest and most effi- 
cient devices for recording different types of informa- 
tion; the selection of computer and programming de- 
vices; the determination of the optimum operating con- 
ditions for a given machine, The groups of problems 
considered in the work of this Section encompassed the 
main trends in the development of programmed control 
techniques for machine tools, as well as a number of spe- 
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cial problems connected with the operational reliability 
of machine tools, 

In the discussion following the reading of reports, it 
was noted that digital-control methods play an especial- 
ly important role in machine construction. 

Twenty-two reports were discussed in the Section 
occupied with the Automation of Power Systems. 

An analysis of the reports submitted to this Section 
indicates thata number of new and important scientif- 
ic and practical results have been achieved. The most 
important are the following: Equations which determine 
the main characteristics of the most important power- 
generating systems to be controlled (boiler units, nuclear 
reactors, and heat exchangers) have been derived, Ex- 
perimental data characterizing the properties of typical 
systems to be controlled (industrial furnaces, power-gen- 
erating devices, etc.) have been obtained, New criteria 
for the optimum operating conditions of systems to be 
controlled, which take into account the effect of the 
controlled variable deviations ard of changes in the en- 
ergy carrier flow on the efficiency of devices, have been 
proposed. A number of original schemes for the automa- 
tion of the most important systerns to be controlled (boil- 
er units, nuclear reactors, metallurgical furnaces, and 
fractionating columns) have been described. New views 
have been expressed on the principles of operating mode 
optimization to be applied to controlled systems, and 
first attempts have been made to determine the scope of 
different types of self-adjusting systems. Specifications 
for the characteristics and parameters of certain systems 
to be controlled, based on control conditions, have been 
formulated and substantiated. 

Nevertheless, the discussion showed that, at the pre- 
sent time, no substantial progress toward solving the fol- 
lowing important problems in industrial automation has 
been achieved: general principles for the design of auto- 
matic-control systems for complex industrial combines 
have not yet been formulated, and methods for improv- 
ing the operational reliability of automatic-controi sys- 
tems have yet to be devised. 

Eight reports were presented to the Section concern- 
ed with Automation in the Chemical and Oil-Refining 
Industries, The main attention in the work of this Sec- 
tion was centered on two problems. The first pertained 
to the determination of dynarnic properties of systems to 
be controlled in the chemical industry, and the second 
problem involved the application of computers for the 
analysis, calculation, and control of technological pro~- 
cesses, The considerable number of reports devoted to 
the study of dynamics of systems to be controlled and 
the animated discussion called forth by them convincing- 
ly prove that the development of computer techniques 
does not lessen the importance of determining the dy- 
namic properties of systems to be controlledin analytical 
form. The study of characteristics of systems to be con- 
trolled still remains one of the most important stages in 
designing complex automation systems, 
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However, the problems encountered in the applica- 
tion of computers to automatic-control systems for tech- 
nological processes were considered only in their most 
general form, indicating that computer techniques have 
not yet been introduced in industrial automation on a 
large scale, 

In the reports and the discussion, attention was 
drawn to the necessity of directing efforts toward the 
solution of the following important problems; determin- 
ation of dynamic properties of systems to be controlled 
and of entire technological chains; development of op- 
timization algorithms for technological processes; estab- 
lishment of efficiency criteria for a justified application 
of computers and the development of general principles 
of complex automation. 

Ten reports were discussed in the Section dealing 
with Automated Electric Drives and Electrical Machines. 
The greatest interest was aroused by problems met in 
seeking new methods for accelerating transient processes 
in automated electric drive systems, problems in the de- 
sign of new components and devices, and methods for 
the synthesis of nonlinear automated drive systems. 

The reports and the discussion emphasized the ne- 
cessity for a further development and a more widespread 
application of methods for the analysis and synthesis of 
automated drive systems with optimum quickness of re- 
sponse, the necessity for a more widespread application 
of computers in systems for the control of electrical 
drives, the need for devising new compact and reliable 
automation components and devices, and the necessity 
of developing improved measuring instruments of high 
accuracy. 

Thirteen reports were discussed in the Section occu- 
pied with Automation of Metallurgical Processes. Among 
the problems discussed were the determination of static 
and dynamic characteristics of systems to be controlled, 
problems met in the design of new automatic control 
systems for complex metallurgical units, and problems 
in the design of new, improved components for automa- 
tion systems. Equations were adduced which fairly ac- 
curately described the static and dynamic properties of 
a number of important metallurgical devices (rolling 
mills, boiling-layer furnaces, heating wells, etc.); a 
large number of pertinent experimental data also were 
given, 

It was announced that new control systems, based on 
self- adjustment principles and the extensive use of com- 
puter techniques, were developed and tested in industry. 
Methods for the complex automation of individual units 
and of entire shops and plants, where process optimiza- 
tion requirements, determined by a number of techno- 


logical and efficiency specifications are taken into ac- 
count, were outlined in these papers. Information was 
provided on the development of new original and effi- 
cient components in automatic-control systems for in- 
dustrial processes (electron-optical equipment for the 
measurement of dimensions and rate of rollings, pneu- 
matic meters with vane-shaped profiles, thermoelements, 
etc.). 

The discussion that developed in the work of this 
Section showed that, for the further development of au- 
tomation in metallurgy, it will be necessary to devise 
methods for the mathematical description of character- 
istics of systems to be controlled, to develop unique 
principles for the design of automatic-control systems, 
and to establish performance criteria for automatic-con- 
trol systems to be used in complex metallurgical de- 
vices, It was also indicated that methods of automatic 
control with respect to disturbances, statistical methods 
for determining disturbances acting on a system, and 
self-adjusting systems for devices whose operating mode 
is optimized by means of models should be developed 
and utilized on a larger scale; that it is necessary to fur- 
ther improve and increase the reliability of components 
in automatic-control systems for metallurgical devices, 
which is especially pertinent to data transmitters, final 
mechanisms and regulating components, since they op- 
erate under heavy-duty conditions in metallurgical 
plants, 


Four reports were discussed in the Section dealing 
with Unclassified Problems, The greatest interest was 
shown in reports devoted to problems in control and reg- 
ulation in the cement and sugar industries. It was re- 
marked in the discussion that these reports contain many 
valuable results obtained in the field of automation in 
sugar refining and cement production; however, it was 
also noted that the complex processes occurring in blast 
furnaces have not been thoroughly analyzed in develop- 
ing automatic-control systems to be used in the cement 
industry. 

In conclusion, it should be noted that the numerous 
discussions which developed during the work of the Con- 
gress were entirely businesslike, that they contained a 
great deal of well-meant criticism, and that, conse- 
quently, they were fruitful. 


The work of the Congress, the slogan of which was 
"The theory should be widely applicable, the technology 
should have the maximum reliability, and the automa- 
tion systems should be highly efficient,* showed that 
this slogan very accurately characterizes the modern 
trends in automatic-control development, 
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Some fundamental problems in communication theory and control theory are compared, 
The problem of designing an optimum (in the statistical sense) closed-loop dual-control system, 
is formulated. Its solution, as well as examples and some generalizations, will be given in parts 


I, Ml, and IV. 


Introduction 

A general block diagram of signal transmission, as 
investigated incommunication theory is shown in Fig. 
1, The transmitted signal x* proceeds from the trans- 
mitting device A to the communication channel H*. 
The mixing of signal and interference (noise) h* now 
takes place. The resultant signal y* represents a mix- 
ture of the transmitted signal and the interference. The 
resultant signal proceeds to the input of receiver B. The 
optimum receiver problem consists in obtaining the sig- 
nal x at its output,such that it is, in a specified sense, 
closest to the transmitted signal x* or to some transfor- 
mation of the signal x*. The mathematical side of the 
problems related to such systems has been the subject of 
important investigations by A. N. Kolmogorov [1], N. 
Wiener [2], C. E. Shannon [3], and A. Wald [4]. This 
type of system was investigated in the works on com- 
munication theory by V. A. Kotel"nikov [5], D. Middle- 
ton, D, Van Meter [6], and others. The cited works dif- 
fer in their various approaches to the problem, but are 
all basically concerned with the investigation of the 
scheme represented by the block diagram in Fig. 1. The 
results obtained in the above-cited works, and in parti- 
cular the Kolmogorov-Wiener theory, have proved use - 
ful in formulating the statistical theory of automatic- 
control systems. This theory has been expounded in the 
books of V, S. Pugachev [7], J. H. Laning, Jr.,and R. H. 
Battin [8], and others, The fullest consideration has been 
given to the theory of linear systems, If a system is lin- 
ear, then whatever the <losed-loop system, it is easy to 
obtain an open-loop system equivalent to it. That is 
why the automatic-control systems,which, as a rule, are 
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closed-loop systems, enable one to use the scheme as in 
Fig. 1, provided that the system is linear. Some com- 
plications and difficulties arise when the interference 
does not appear at the input of the system, but at the in- 
put of the controlled object —the latter being inside the 
closed-loop network, This also creates difficulties which 
are not, however, of a fundamental nature. More seri- 
ous difficulties arise due to the bounds to which the pow- 
er of the system's signals are subjected. The problem 
becomes more involved when the controlled object is 
nonlinear or it 1s required that an optimum control sys- 
tem, which often proves to be nonlinear, is designed. It 
is not always possible to proceed from an open nonlinear 
system to an equivalent closed-loop one; furthermore, 
this is an extremely involved process. In such a case, 
the open-loop scheme depicted in Fig, 1 cannot be used 
in practice. A number of attempts have been made to 
reduce approximately nonlinear systems to equivalent 
linear ones (see, for example, the paper of I. E. Kazakov 
{9}. Such studies are of considerable practical value 
but do not, in their present state, provide a means of es- 
timating how close the obtained approximation is to the 
true solution; neither do they enable one to synthesize 
the optimum system. 

In order to be able to solve optimum problems of 
the control theory, a fundamentally different approach 
is required, Firstly, a different block diagram replacing 
that depicted in Fig. 1 is needed, Before selecting a 
common scheme to be used in automatic-control 
theory, it seems advisable to have a preliminary survey 
of certain basic concepts of the theory. 

Figure 2, a shows the controlled object B with x as its 
output, u as the controller, and z as the disturbance (in- 
terference). When the system has several inputs and out- 
puts, one can regard x, u, and z as vector quantities. 

The output x depends on u and z. This dependence 
can be described either by a linear or a nonlinear opera- 
tor, and in a particular case of memoryless systems, only 
by a function. The interference z is generally a function 
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of time. Thus, since a change in the system's character- 
istics can be considered a particular result of interfer- 
ence (e.g., the parametric effect), then hereafter every - 
thing in the system's characteristics that changes with 
time will be attributed to interference. If, for example, 
x depends on u as in 


v == (ag + fo(t)]u2 + fa, +f, ()) a + ag + fe (t), (D 
then the vector 


z= {Zo, 245 29), 


Zo = fo(t), = fit), 22 = /2 (8) (2) 


gives the interference or disturbance, and the formula 


# = (ay + 2] u? + [a, + 2,]u + a, + 2, (3) 


represents a particular operator. 

If the object has memory and if, for example, its 
motion can be described by a differential equation of 
the oth order, then its state is considered as one of the 
characteristics of the object as described by the value of 
the vector x in the n-dimensional phase space. 
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Complete information about a system thus consists 
of information about its operator, interference (noise), 
and the state of the system. The controlling device may 
be given and is considered as known, The open-loop 
systems are automatic-control systems of a simple type. 
A block diagram of an open-loop system is shown in 
Fig. 2,b and is of the same character as the one given in 
Fig. 1. The exciting quantity x* enters the input of the 
regulating member A, determining how the output quan- 
tity x should vary. The output u of the regulating de- 
vice enters the input of the object B under control. The 
output x of the controlled object B does not proceed in 
this case to the regulating member A. 

The required rule of change in x can only be im- 
plemented when full information about the controlled 
system is available, i.e., when its operator and state x 
are known, at least at the initial moment of time, as well 
as the interference z. The latter should be known a pri- 
ori at all moments of time, including future moments. 
The required rule of change in x must be one of the ad- 
missible ones,such that it can be implemented for the 
given class of initial states of the system and for a class 
of controlling motions u staying within acceptable 
bounds. 

The above conditions, and, in particular, the a pri- 
ori full knowledge of the controlled system, cannot be 
satisfied in practice. This is why an accurate implemen- 
tation of the required control rule cannot be obtained. 
Sometimes the interference z is not known a priori, but 
it is possible to measure it with a device C (see Fig.2,c) 
and to introduce the outcome of the measurement into 
the controlling member A. One can then find in the lat- 
ter the required controlling rule u. Such a scheme is al- 
so an open-loop one. But the scheme depicted in Fig. 
2,c differs in some ways from the scheme with complete 
a priori information about the system, as now future 
magnitudes of the interference z remain unknown. Be- 
cause of this, the exact implementation of the required 
rule of variation of the controlling quantity x is not al- 
ways feasible. 

When the state x of the system is not known, then, 
generally speaking, it is not possible to implement the 
required rule of the change in x. To be able to attain 
the required variation in x, or one near to it, a feedback 
network is needed to feed the output quantity x to the 
input of the controlling member A (see Fig.2,d). Having 
compared x and x*, the controlling member generates 
the regulating action u,bringing x to its required value. 
The block diagram ofFig.2,d isa closed-loop scheme, 
and is of the utmost importance in the automatic-control 
theory. 
A closed-loop network offers far-reaching possibili- 
ties not available in an open-loop system. For example, 
it may be possible for a class of objects B of control to 
obtain a process x close to the one required even 
when the interference z remains unknown and incapable 
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of measurement. Let, for example, the interference z, 
together with the controlling action u, be applied at the 
input of the controlled object B, the latter representing 
an inertia member. If the quantity x* =x is to be attain- 
ed, the regulating member can be implemented in the 
form of an amplifier of high gain k >> 1, the difference 
x* ~x being sent to its input. 

It is not difficult to see that the requirement x® = x 
will be satisfied approximately, whatever the continuous- 
ly varying interference z, provided k is sufficiently high 
and the bounds of variability of u are such that the in- 
terference z can be compensated by u. 

This principle of neutralizing the interference can 
be generalized and applied to cases considerably more 
involved, combining the system's accuracy with its sta- 
bility. A detailed analysis of the applicability of this 
principle was carried out by M. V. Meerov in his mono- 
graph [10). 

When the interference z can be measured, it is pos- 
sible to implement a combined system (see Fig. 2, e) of 
measurement of the state of the controlled system x as 
well as of its interference. Such systems are of consider- 
able practical value. We shall not, however, concern 
ourselves with them but shall limit the study to the 
"pure" type of closed-loop systems, considering ther as 
being of primary importance. 

The input quantity x* may be previously unknown; 
usually, neither do we have any prior knowledge of the 
interference z. Consequently, these processes become 
random, and, in a favorable case, the a priori informa- 
tion is limited to our knowledge of their statistical char- 
acteristics. Such processes may be regarded as belonging 
to a class of curves x* (AX) and z(p) where A and p are 
parameter vectors(A4,. . .» A,)and (444, ..., Hy)» respec- 
tively, with their probability distributions either known 
or unknown. 

In communication channels connecting the blocks of 
a system, the errors of measurement or noise can be re- 
garded as subsidiary random processes as well, with eith- 
er known or unknown characteristics. Thus, the analysis 
of a control system and the synthesis of the regulating 
member can be regarded as problems of a statistical na- 
ture. The problem should be solved for an over-all block 
diagram in which all the above features of an automatic- 
control system are reflected, Such a block diagram is 
depicted in Fig. 3; it is the subject of the present paper 
as well as that of further papers in this series. 

The input quantity x* proceeds to the input of the. 
controlling member A through channel H® where it be - 
comes mixed with noise h*. Thus, the quantity y* en- 
tering the input of A is generally not equal to the actu- 
al value of the input quantity x*. There also exists a 
class of systems with the external input x* altogether 
absent. Generally speaking, however, it cannot be neg- 
lected. A similar mixing takes place of the state x of 
controlled object B and noise h in channel H; quantity 
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y entering A will not, as a rule, be equal to x. The reg- 
ulating action u proceeds next from A to controlled ob- 
ject B having previously passed through channel G where 
it was mixed with noise g. The quantity v proceeding 
to the controlled object is not, as a rule, equal to u. 


Dual Control 





One cannot neutralize, in a general case, the inter- 
ference z by a regulation u if the interference z is not 
known. Its direct measurement is not, as a matter of 
fact, often possible. In such a case, an open-loop sys- 
tem is useless. But the closed system in Fig. 3 shows 
how z can be indirectly determined by measuring the 
input and the output, in and out of object B, by studying 
its characteristics, The input of controlling member A 
enters both the input v and the output x of the object or, 
in any case, the quantities u and y related to v and x. 
The examination of the quantities u and y provides in- 
formation on the characteristics of object B. It should 
be understood that this information is never complete, as 
the noises g and h render an exact measurement of B's 
characteristics impossible; if the actual form of the ob- 
ject’s operator is not known either, a full determination 
of its characteristics would not be possible even in the 
absence of noise, unless the determination time is infin- 
itely great. The lack of complete information on the 
disturbance z can assume the form of ana posteriori prob- 
ability distribution of its parameters. Although the lat- 
ter does not provide precise values of the parameters, it 
is more accurate than an a priori distribution, as the 
former reflects the real character of the interference. 

If the random process can be measured directly, one 
is able eventually to specify its statistical characteristics 
more accurately. The method which provides such im- 
provement with the aid of dynamic programming was 
discussed in examples by R. Bellman and R, Kalaba [11, 
12) and also by M. Freimer (13). One is able to find 
the characteristics of the process x* more accurately in 
the open part of the block diagram in Fig. 3, or in a 
similar scheme in Fig. 1. 
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This formulation of the problem is characteristic for 
an open system. In a closed system its formulation be- 
comes totally different. It is shown that some processes 
in the system of Fig. 3 may occur which have no counter- 
part in open-loop systems. Whereas open systems can on- 
ly be studied by passive observations, the study may de- 
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velop into an active one in closed systems. In order to 
improve the investigation one may vary the signals u 
(or v) which act on the controlled object B, The object 
is, as it were, “reconnoitered” by signals of an enquiring 
character whose purpose it is to promote a more rapid 
and more accurate study of the object's characteristics 
and of the methods of controlling it. 

However, the controlling movements are necessary 
not only to study or to learn the characteristics of the 
object or the ways of controlling it, but alsoto implement 
the regulation, to direct the object to the required state. 
Thus, the controlling effects in the general block dia- 
gram in Fig. 3 must be twofold: they must, to a certain 
extent, be investigating as well as directing. 

The control whose regulating effects are of this 
twofold character will in the sequel be called dual con- 
trol; the papers in the present series will be devoted to 
the theory of dual control. 

Dual control is particularly useful and even indis- 
pensible in cases where the operator and the interference 
z in the object B are complex, and the object is thus dis- 
tinguished either by its complexity or by the variability 
of its characteristics, Some typical examples of systems 
with dual control are to be found in automatic search 
systems, in particular, in automatic optimization sys- 
tems (see, for example, [14 and 15]}). In these systems, 
the investigating or “trial” part can usually be separated 
easily from thecontrolling or “operating” part of the sig- 
nal, either by the difference in their frequency ranges or 
because they interweave in time. Such a separation, 
however, need not always take place; an effect can be 
twofold in character by virtue of being partly diagnostic 
and partly regulating. 

Thus, in dual-control systems, there is a conflict 
between the two sides of the controlling process, the in- 
vestigational and the directional. An efficient control 
can only be effected by a well-timed action on the ob- 
ject. A delayed action weakens the control process, But 
the control can only be effective when the properties of 
the object are sufficiently well known; one needs, how- 
ever, more time to become familiar with them. A too 
“hasty” controlling member will carry out the operation- 
al movement without making proper use of the results of 
trial investigations performed on the object. A too “cau- 
tious” system will bide its time unnecessarily long and 
process the received information without directing the 
object to its required state at the right time. In each 
case, the control process may not prove the best one and 
may not even prove to be up to the mark. Our problem 
is to find out, one way or another, which combination 
of these two sides of the regulation would prove to be 
most suitable. The operations must be so selected as to 
maximize a criterion of the control's quality. 

As shown above, the incomplete information about 
the object will be expressed by the presence of the prob- 
ability distributions of potentially possible characteristics 


of the object. The regulating member compares, as it 
were, the various hypotheses on the object, with proba- 
bility of its occurrence being attached to each hypothe- 
sis, These probabilities vary with time. There may be 
a control method such that the most probable hypothesis 
will always be selected and, therefore, assuming that it 
is valid, the optimum control method will be attained. 
Such a control system is not generally optimum in the 
absolute meaning of the word as the complete informa- 
tion on the object has not been utilized. The probabili- 
ty distribution of the different hypotheses extracted from 
the experiments is distorted as the probability 1 was as- 
cribed to one of them and the probability 0 to others. 

A better control method will be one whereby the proba- 
bilities of all the hypotheses would be taken into ac- 
count, 

The probability distribution of hypotheses will vary 
with time, the higher probabilities concentrating more 
and more in the region of those hypotheses which ap- 
proach the true characteristics of the object. The pace 
of concentration and, therefore, the success of the sub- 
sequent regulating movements, depends on the charac- 
ter of the preceding regulating movements, on how well 
they have “sounded” the object. Thus, two factors 
should be taken into account by the controlling member 
which decides the specific amount of regulating move- 
ment at any given moment of time: 

(a) The loss occurring in the value of the quality 

criterion due to the fact that the outcome of the 

operation at a given moment, and at subsequent 
moments of time, will cause a deviation of the ob- 
ject either from the required state or from the best 
attainable one, The average value of this loss shall 
be called the action risk. 

(b) The loss occurring in the value of the quality 

criterion due to the fact that the magnitude of the 

controlling action has not proved the best to obtain 
information on the characteristics of the object; in 
view of this, the subsequent actions will not be the 
best posssible ones either. The average value of 
this loss shall be called the in ation risk, 

It will be shown that for a certain class of systems, 
the total risk will be equal to the sum of the action and 
investigation risks, 

All systems of automatic search (see [14])are char- 
acterized by trial actions, Dual control, therefore, is 
applicable to all systems of automatic search and, in 
particular, to automatic optimization systems. It can 
also be applied to other types of closed-loop systems 
which do not belong to the automatic search class at all, 
To illustrate the difference between the two types of 
dual-control systems, a few examples will be given, 

Figure 4,a showsa system which operates as follows: 
the main regulating member A implements the control 
of object B, either in an open- or in a closed-loop net- 
work (the closed one is indicated by a dashed line), The 
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regulating movements u are of the investigational and 
action type simultaneously, The quantities u and x, 
from the input and the output of object B, respectively, 
enter an additional controlling member A‘, The latter 
receives the characteristics of object B from the results 
of the investigations; subsequently, in accordance with 
an algorithm given in advance and fed into the device 
from outside, the parameters of controlling member A 
are so computed that its controlling action is optimal. 
Having the results of the computation, the additional 
regulating action w establishes the computed optimum 
parameters in the main controlling member A. This pro- 
cess may repeat itself periodically. 

Such systems contain investigational movements; 
however, the automatic search is absent. The parameters 
of member A are established, not via automatic search, 
but from an algorithm given in advance, from a function 
of the determined characteristics of object B, There is 
no investigational component in the operation w. The 
channel w is not usually found in a closed network, as the 
change in the coefficients of A has no effect on the co- 
efficients of B, 

A block diagram of an automatic optimization sys- 
tem is presented inFig.4,b, Here the action u is dual in 
character, investigating object B as well as directing it 
to its optimum mode of action. The latter corresponds 
to an extremum of quantity Q dependent on x. The op- 
timum mode is found by means of automatic search. 

The latter is conducted in such a way that the informa- 
tion received from the investigating action u and from 
the output x of the system is analyzed in controlling 
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member A. This permits determination of the regulating 
part of the same action u on the same object B whose in- 
put and output were investigated, so that the same quan- 
tity x which was being investigated can be changed in 
the right direction. 

The combination of the investigating and directing 
operations has, so far, not constituted the whole search, 
but only one of its distinct features. There is no search 
in the system depicted inFig.4,a; it takes place, how- 
ever, in the system of automatic optimization shown 
inFig.4,b. But both systems are of the dual-control 
type. 
InFig.4,b the controlled object is inside a dashed 
rectangle denoted by O. In this case, it does not differ 
from object B. The dual control, however, can also be 
applied to control the entire automatic system, consid- 
ered as a complex object. For example,inFig.4,c the 
complex object O inside the dashed rectangle comprises 
the controlling member A and the object B of control. 
The auxiliary controlling member A‘ investigates the 
process x and, with the aid of the controlling process w, 
can vary the algorithm of regulation implemented by 
member A, The w processes are twofold in character. 
The investigation of changes in the algorithm of mem- 
ber A and their effect on the process x results in regulat- 
ing processes w, bringing the algorithm of member A to 
such a form that the process x will either prove ad- 
missible , favorable, or optimum, depending for what 
purpose the system will be used. Here an automatic 
search takes place in the closed network of processes 
wox7u. 

Statement of the Synthesis Rroblem 
of an Optimum System of Dual Control 








The problem of designing an optimum, in the spec- 
ified meaning of the word, controlling member A, as 
shown in Fig. 3, is formulated below. It is advisable 
when formulating the problem to make use of certain 
concepts of the theory of games and of A, Wald's [4] 
theory of statistical decisions (see D. Blackwell and M. 
A, Girshick [16], and also Chow [17]). In solving the 
variational problem as stated later in the present series 
of papers, use is made of the concepts of R. Bellman's 
dynamic programming (see, for example, [18]). In the 
subsequent parts of this series, the mathematical expo- 
sition may appear somewhat cumbersome but is actual- 
ly quite simple. The main contents of the papers deals 
with further development of the concepts of automatic 
control briefly described above. 

Cousider the scheme presented in Fig. 3. The fol- 
lowing limitations of the statement of the problem are 
introduced, 

1) A discrete-continuous system is investigated in 
which the time but not the level is quantized, All mag- 
nitudes occurring in the system are considered at dis- 
crete moments of time t=0, 1, 2,..., monly. Any 





is 





magnitude at the sth moment of time will carry the in- 
dex s. Thu, the considered quantities are x*, » Xgo Yoo 
V9 Bes etc. 


Such limitation enables one to simplify the compu- 
tation. Moreover, in many cases this actually occurs. 
The transition to the continuous time can in some cases 
be accomplished in an intuitive manner by making the 
time interval between the discrete values approach zero 
(see Part III). One meets with considerable difficulties 
in more fully examining the passage to the limit. 

2) The time interval, or the number of cycles n 
within which the process is being investigated, is as- 
sumed to be a fixed constant. In certain cases no major 
difficulties arise when proceeding: to the limit with 
n+ ©, A wider generalization relating to a variable 
number n of cycles not known beforehand would be of 
interest, but will not be tackled in the present paper. 

3) A Bayesian problem, in which a priori densities 
of random variables are given, is considered, Other 
formulations, for example, minimax, are also of con- 
siderable interest, but far more difficult to solve. This 
problem could also be formulated in relation to the 
concept of the so-called “inductive probability® (see, 
for example, the paper of L. S. Schwartz, B. Harris, and 
A. Hauptschein [19]). 

We assume that h%, hg, gg, are sequences of inde- 
pendent random variables with identical distribution 
densities P(h% )» P(h,), P(g,). Further, let z,.=2(s, ft), 
and xf =x(s,X) where and A are random parameter 
vectors with coordinates Hy and A »Tespectively: 


= (pr,.--, pm)» A= (Ag,--0yAghe (4) 


The a priori probability densities P(g) and P(A") 
are given. 

4) The object B is assumed to be memoryless, in 
other words, the values x, of its output depend only on 
the values of the input quantities z, and y, at the same 
moment of time: 


~.= F, (Zs, Us). (5) 


The functions Fy and z, are assumed to be finite 
and single-valued, continuous and differentiable. 

A generalization relating to objects with memory 
and with x, depending on x,, Z,, ¥, (t <8) will be given 
in Part IV, It should be pointed out that memoryless ob- 
jects are of great practical value. Namely, if the input 
data (initial conditions or values of parameters) are giv- 
en for a certain model, and one is able to carry out ex- 
periments using this model and also to register the re- 
sults, then such an object becomes equivalent to a mem- 
oryless one. 


5) A simple criterion W of quality is introduced, 
Let the partial loss function corresponding to the 


gth time moment be of the form 


W, = W (s, Zs; z°). (6) 


Moreover, let the total loss function W for the total 
time be equal to the sum of partial loss functions (such 
a criterion shall be called a simple one): 


sn 


W = > W (s, 22, 2s). (7) 


8-=0 


The smaller the mathematical expectation of W, 
the better is the system, It shall be called optimum 
when its average risk R (i.e., the mathematical expec- 
tation M of the quantity W) is minimal, The amount of 
risk is given by the formula 


R= M(W) = 
(8) 


= u{> W (s,2%,2)} = >) M{W,) = 5 R,. 


s=0 


Each R,=M {W,} will be called a partial risk due 
to the sth cycle. 

There may be many types of simple criteria, for 
example, 


W, = a(s)(z,—2;]*. (9) 


Criteria of practical importance need not always be sim- 
ple,and generalizations relating to other criteria would 
therefore be of interest. 

The formulation of the optimum strategy problem 
in terms of risks is not the only one in existence. There 
exist a number of studies in which closed systems are in- 
vestigated from the point of view of the information 
theory (see, for example, R. L. Dobrushin's paper [20]). 
As the primary aim of a control system does not lie in 
transmitting information but in designing required pro- 
cesses, the formulation of the problem in the language 
of statistical decisions fits in better with the intrinsic 
nature of the problem. 

6) All the quantities occurring in the sth cycle will 
be regarded as scalar. Our object, therefore, has only a 
single input v and a single output x. The exposition be- 
comes more involved with generalizations relating to 
objects with several inputs and outputs (see Part IV). 

7) We assume that the manner by which the signal 
and the noise are combined in H*, H, or G blocks is 








known and invariable, and that the blocks are memory- 
less, Thus, 
Varo (u, 8), 
y= y'(h’, 2"), 
y=y (h, z). 


(10) 


Therefore, the conditional probabilities P(y*|x*) 
and P(y|x) and P(v| u) make sense. 

8) We assume that the controlling member A gen- 
erally possesses a memory and that, moreover, for the 
sake of generality, the algorithm of its action is a ran- 
dom one, i.e., the part A exhibits random strategy. 

We introduce the vectors (0 =s =n): 





U, = (Uo, Uy, ... Us), 
y; = (y;, y sees y;), (11) 
I,= (Yo Yist-+s y,)- 


The controlling member can now be characterized 
by the probability densities 


P, (u,) =, (us, usa, ¥;, y,,) (@S#s<”)- (12) 
The problem consists in finding a sequence of func- 
tions F such that the average risk R (see [8]) becomes 


minimal. 
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This article presents an analysis of a continuous control system where the extremum of the 
characteristic of an inertialess system to be controlled is searched for and maintained in the 
presence of intensive random noise in the frequency band which includes the frequency of the 
periodic trial input. The stability conditions for this system are determined. 


The spectrum of uncontrolled external signals acting An inertialess system of the simplest form with high- 
on an extremal controlsystem usually contains a low - frequency noise at the input of the system tobe controlled 
frequency band, which is determined by the drift of ele- is considered below. The block diagram of this sys- 
ments in the system, and a low-frequency band deter- tem is given in the figure. 
mined by noise. Between these bands lies a frequency If the characteristic of the system to be controlled 
band which is almost free from noise. However, this is given by the quadratic parabola y=—x’, the system 
third band may be insufficiently broad to serve for as- operation will be described by the following equations: 
signing the frequency of trial inputs[1, 2]. Therefore, 1) The equation of the system to be controlled, 
the case wherein the trial input frequency lies within y = —(z+asinet + 2)’, 
the noise band is of especial interest. However, if we 
assume that the high-frequency noise appears only at where a sin wt is the trial input, and z is the noise; 
the output of the system to be controlled, no basic dif- 2) the demodulator equation, 
ficulties are encountered in analyzing the equations of 
the closed-loop system [2]. For a controlled system rH 
whose characteristic is given by a quadratic parabola, _ bo : - 
the control system is described by linear equations with a. \ spans @) 


t 
constant coefficients, For an arbitrary characteristic, the 


problem of determining the main deviation of the system 
operating conditions from the extremum can also be 
solved by using the statistical linearization method. 
However, if high-frequency noise appears at the system ( 
0 


3) the final mechanism equation, 


input, the coefficients in the equation for the closed- Pee: Few dt. 


loop system become random functions of time. The in- 


vestigati f such ati is a licated task. 
a ee ee een We shall limit our consideration to the case wherein 


the spectrum z(t) of the normal random noise has a nar- 
row band and the noise can, consequently, be represented 
in the form 


z(t) = A(t) sin wt + C (¢) cos of. 
(2) 
Here, A(t) and C(t) are slowly changing (in com- 
parison with the trial input) random processes, By taking 
into account condition (2), system (1) is reduced [2] to 
the equation 





1) System to be controlled; 2) 
generator of trial oscillations; 3) ae {a+ A(t)|jz=0 
demodulator; 4) final mechanism. :: 


(3) 
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or 


“= +E(t)a=0*, (4) 
while 


me = bia -+ ma), B: (T) — Ba (t), (5) 


where Mas B,, m,» and B are the mathematical expec- 


tations and the autocorrelation functions of the A(t) and 
£ (t) processes, respectively. 

Thus, the presence of high-frequency noise at the 
input of the system to be controlled causes the appear- 
ance of a parameter which changes in random fashion in 
time in the equation which determines the variable x of 
the control mechanism, According to this method [3], 
after correcting the error contained therein, we shall ob- 
tain the following equations for the mathematical expec- 
tation and dispersion: 


t ty 


\ B: (t) dx| (t — ta) x 


my (t) == Ly exp \— | m: 


\e. (t) dr|, (6) 


0 


xX exp |- 


i—te 


3; (t) = ziexp{—-2| me! -— 
x exp | aes (t) dr| — 
uv 
t—t, 


— 7, exp {—2| mz — \ Bz (t) dx| \« — te} x 


a 
x exp|--2 \ «Be (x) dr. (7) 
J J 
Here, X = Xo for t = to. 
The values of m, and o, for t= obviously depend 
on the sign in front of the quantities: 


foe) 
h = m; —\ Be (t)dt, p= m: — 2\B: (t)dr.(8) 


ota 


eo 
Since T= JB, (tT) dt > 0 and p = A-T, the following 
0 


variants can be realized: 
I A>0, p> 0; 
Il A>0, p=0; 
Ill A>0, p<0; (9) 
IV. A=0, p<0; 
Vv. A<0, p<. 


In this connection, the following cases can occur 
under steady-state conditions; 


2 | Bz (x) de| (t — tx 


I. m,(%)=0, 0% (%)=0 for any initial conditions; 


II, m,(*)=0 for any initial condition, and 
co 
63 (co) = 25 exp |—4 \ tB: (t) dx |; 
0 


III, m,(*)=0 for any initial condition; for xy» 0, 


oF (t) increases boundlessly with time or, in other words, 
the system is unstable with respect to dispersion; 


IV. m, (%)=x9exp[— {rB, (T)dT], and the system 


is unstable with respect to dispersion; 
V. The system is unstable with respect to the aver- 
age value as well as with respect to dispersion. 

Hence, it is obvious that the dispersion is limited 
under steady-state conditions if y =0. For a system 
described by Eq. (3), this condition is of the following 
form: 


b(a + m,) — 2b*( Ba(t)dt>0 (10) 
0 


or, for b > 0, 


@ > 2bt 0% —m,. 11) 


Here, q = B (0) is the daperticn of the random 


process nels and " =(1/B (On, (T)dr is the time of 


the correlation gate A(t). 

These conclusions were verified by experiment,T 
For b=1, Eq. (3) was reproduced by means of an elec- 
tronic simulating device, The random noise was secured 
by a low-frequency noise generator. In this, the cor- 


relation time was T, = 0.25 sec, and on =8v, 


A decrease in the a parameter, which corresponded 
to the amplitude of trial oscillations, led to “stability” 
losses, In this, as was to be expected from relation (7), 
which was written for Eq. (3), the “stability” disturbance 
took place aperiodically, i.e., for a small deviation 
from the a=4 limit, given by condition (11), the disper- 
sion o changed very slowly. The limiting value itself, 
which was found by experiment, coincided almost ex- 
actly with the theoretical value. 


* Equations of the type (4) were used in paper [3]. 
However, due to the error in calculating the character- 
istic function, the results obtained in this paper could 
not be used directly in this case. The error in [3] was 
detected by V. I. Osorin in the reference journal Me- 
khanika 10, 11232 (1959), 


T E. P. Gil*bo performed the experiment. 














Thus, the stability loss in scanning and in main- 
taining the extremum, caused by the presence of high- 
frequency noise at the input of the system to be con- 
trolled, is entirely possible. This limits the possible 
amplitudes of trial oscillations and, consequently, the 
maximum allowable scanning loss. 

Regardless of the simplicity of the problem as for- 
mulated at the outset, its investigation made it possible 
to obtain a definite physical result which can be applied 
to all continuous extremal systems operating in accord- 
ance with the gradient method. It consists in the fact 
that the consideration of random noises leads to the 
necessity of a more complicated approach to the sta- 
bility problem in such systems than was used, for in- 
stance, in [4], where the entire problem was reduced to 
the investigation of linear systems with constant coef- 
ficients. 


If the noise level at the input of the system to be 


controlled is high, the following of the extremum may 
be impossible due to the noise parametric action. 


1, 
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This article presents an analysis of extremal relay control systems under self-oscillation 
conditions with fluctuations at the regulator input. A method of statistical linearization of such 
systems, based on the reduction of externa] random disturbances to parametric action on the 
regulator, is proposed. This method makes it possible to obtain an expression for the equivalent 
gain in the form of a random step function of time according to the describing-function method. 
The average value of this function is taken as the statistically equivalent transfer coefficient of 


the nonlinear component. 


This method can be applied to ordinary nonlinear self-oscillation systems. 


In extremal relay systems (ERCS), the self-oscilla- 
tion mode is the operating mode. Under operating con- 
ditions, an actual system is usually subject to the action 
of external random disturbances. The present work ana- 
lyzes the behavior of ERCS under self-oscillation condi- 
tions,when € (t) fluctuations are applied to the regulator, 
input. The investigation is based on the statistical-lin- 
earization method in conjunction with the describing— 
function method applied in the form developed for ex- 
tremal systems in [1). 

The use of the statistical-linearization method [2] 
in the investigation of systems allowing stable self-oscil- 
lation states was first proposed in [3]. In [2] and [3], by 
operating with the first moments of random functions, 
the linear (equivalent in the probability sense) relation 
between the input and output moments of the function 
at the input and output of the nonlinear component is 
determined on the basis of a certain criterion, while, in 
[3], this is achieved by using the describing-function me- 
thod. This method cannot be applied to extremal relay 
systems, This is due to the fact that, for random distur- 
bances at the extremal regulator input, the parameter of 
this regulator (the relay threshold) becomes a function 
of the instantaneous value ofthe external random dis- 
turbance, Therefore, for nonlinearities of this type, the 
consideration cannot be limited to the random function 
moments at the regulator input. 

Another statistical-linearization method is proposed 
in the present paper. It is based on a nonlinear reduc- 
tion of the random function to a form suitable for intro- 
ducing the regulator fluctuating parameter in the expres- 
sion for the equivalent gain of the nonlinear element, 
which is obtained according to the describing-function 






method, Thus, if, in [3], the equivalent coefficient is 
obtained as a function of the first moments of the exter- 
nal random disturbance, here it represents a random 
function of time. Its average value is found with respect 
to the distribution density of the determined fluctuating 
parameter function. 

The results obtained in investigating ERCS in ac- 
cordance with the above-mentioned idea for the system 
linearization are given below.* 


1. Formulation of the Problem 
and the Method of Solution 








Let &(t) represent a steady-state normal random 
process with zero average value and intensity 0. , which 
is smaller than the amplitude Y of self-oscillations in 
the system, and let the basic components of the fluctua- 
tion spectrum &(t) be comparable to the frequency w of 
self-oscillations. Such a statement of the problem is of 
practical value since the high-frequency fluctuation 
components can always be suppressed by means of a 
weak filter in the control circuit without affecting the 
dynamic properties of the system. 

As is known, ERCS have two basically nonlinear in- 
ertialess components. In approximating the static char- 
acteristic y= f,(x) of the system to be controlled by a 
quadratic parabola, due to the smallness of o, the sys- 
tem to be controlled can be represented in the form of 


a linear element with the equivalent transfer locus 


* In connection with another problem, a statistical-lin- 
earization method, conceptually close to the method 

mentioned above and based on different physical assump- 
tions, was proposed by A. A. Pervozvanskii in [8]. 











K(jw) in investigating self-oscillations [1]. The second 
nonlinear component —the controlling device with 
z=f,(y)— has the ratchet-relay characteristic [4] shown 
in Fig. 1. 

This characteristic can be described by the expres- 
sion 


+k for y $0, y > 0; 
z= + k_ for y +Ky> 0, . (1) 
< 
Sa kp for y +k,<0, vn ms 


Here, k_ is the relay gain (the choice of the sign 
sequence defends on the initial conditions) and K,y=« —y* 
is the regulator fluctuating parameter, which is expressed 
in terms of K (quantization step of the output quantity 
with respect to level) and y*, which depends on the sys- 
tem type. 

Let the input quantity of the control device be ex- 
pressed by the sum 


y = Y sin ot + € (2). (2) 
Then 


y" =Y+ & (¢),_ rT (3) 


for ERCS with a floating quantization grid, and 


y = [Y sin ot + § (2) max (3a) 


for ERCS with a fixed grid of the output quantity quan- 
tization with respect to level in the system to be con- 

trolled [5]. According to (1), we shall write the condi- 
tion for the pertinent switching moment in the system: 


y (tn) + %1 (tn) = 0. (4) 


By substituting (2) in (4), we can reduce the random 
component € (t) of the control device input to the fluc- 
tuating parameter kK, and thus represent the effect of an 
external random disturbance on the system only in the 
form of a parametric action on the regulator. 
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Since the actuation phase B of the relay enters the 
expression for the relay regulator equivalent gain, the 
fluctuating parameter of this regulator can be conven- 
iently transformed into the fluctuating actuation phase 
68 of the relay for a regulator sinusoidal input. As is 
shown below, the final results can thus be geometrically 
interpreted in a clear manner, Thus, the second stage in 
solving the stated problem consists in the nonlinear trans- 
formation of the relay parameter random function into 
the random function 58, with the ultimate aim of deter- 
mining the distribution density of the new random para- 
meter 66. 

Assuming that 58 = const #0 during the period T of 
the regular input first harmonic, the equivalent gain 
S(Y, 58) for the given period can be obtained. It should 
be noted that this operation is valid if the self-oscilla- 
tions in the system are settled within a time interval 
which is less than the period T, As is shown below, in 
ERCS this is practically always the case. 

In the final stage, the average value 5(Y, 68) is de- 
termined, after which, by plotting the graphic solution 
in the usual form given by the describing-function me- 
thod, we can find the average values of the self-oscilla- 
tion regime characteristics in ERCS when the above €(t) 
fluctuations act at the regulator input. 


2. Distribution of the ERCS Relay 
Actuation Phase under Self-Oscillation 
Conditions in the Presence of Fluctuations 


We shall first find the static characteristics of the 
phase distribution for an ERCS with a floating grid of 
output quantization with respect to level. 

The phase distribution will be expressed in terms of 
angles wt, where w is the periodic regime frequency at 
the output of the extremal system to be controlled. 

In the absence of fluctuations, the relay actuation 
phase 8, will be determined by the moment when the 
Yosin uyt signal intersects a certain level K, which is 
measured from the Y, value, i.e., 











Y,sin By = x — Yo, (5) 


where [5] Y<« =2Y, 

For a random external input, considering that the 
€(t) fluctuations enter the control signal linearly, we 
shall seek the phase B® of the sinusoidal component at 
the moment when the €(t,)+ Y%sin6* signal intersects 
the level measured from Y, + €(t,): 


& (ts), _e* +Y,sinB’ =x—[¥o+ E(t). © 


By eliminating (5) from (6), we obtain 








Y, (sin 8’ — sin Bo) + & (ti) +°E (42) = 9. (7) 











In order to determine the statistical characteristics 
of the random function 68 according to Eq. (7), we shall 
use the solution which has been proposed in [6]. 

Let us define more accurately the above restrictions 
which have been imposed on €(t). We shall consider 
that €(t) is a centered function with a Gaussian distribu- 
tion and the autocorrelation function oF R(T), where 


a < Y (8) 


and the correlation time 





B — Bo n 
<<a (9) 


throughout the entire working range of the frequencies w. 

Condition (9) permits us to consider that the first in- 
tersection, the one in which we are interested, is also the 
only one; this condition also determines the independence 
of the €(t,) and €(t,) values of the random function £€(t). 

Under the assumption that the fluctuations are small, 
we shall introduce a small parameter 1, by denoting in 
correspondence with (9): 


(4) +84) =o) =e0 (EL), ao 
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and we write the initial equation (7) in the form 
Y, (sin B’ — sin By) + wt" (=) =0. (11) 


According to [6], we shall seek the solution of this 
equation in the form of a power series of J 


B’ = Bo + vB + w*Be +. -- (12) 


By substituting 6 in (11), we obtain 
Yo {sin (Bp + #Bi + B*Be + .. .)— 
— Sim Bo) + uo" |S (Bo-+4Br+ 8°82 +-..)] = 
We expand the functions of the 8* variable in a 


Taylor series with respect to By: 


Yo {(wB: + H*Bs-|- - - .) 608 By — 
— +B, +H%Be -+.-.)*sin By +. . -_ i. 
0. 


+ her (B) [=e + phe +--)]+- 


By equating the coefficients in front of the equal 
powers of the : parameter, we arrive at the following 
system of equations: 


BY rie yew = 0, 
BeYo cos By — BF . 





Y,sin Bo + 1 


whence we find the successive approximations for 
0< By =(/2) (the region wherein self-oscillations exist 


[5]): 








» (Bo 
onlay 
hae Y,cosB, ’ 
“(Be (hy eran 
B= w2 Y? cos? By + FTW, cos Bo)® 


Consequently, the solution of Eq. (11) with an accuracy 
to #” can be represented by the expression 


By 
¥ 
Bh = B= — py (a) 


(Be po (B 
ALAC OL C a 
TH (@e¥ 9 COS Bo)? ' 2 (¥_ cos Bo)* J: 





According to (10), due to the independence of €(t,) 
and €(t,) and the steady-state property of ¢(t), we final- 
ly write (13) in the form 





=e 5 oO’ i CY, sin By 

> candies Yy Cos By + (@o¥ cos Bo)* + 2 “(Yo cos Bo * 
(14) 

The 68(t) function has a physical meaning only for 
instants of time t, which coincide with the times of ob- 
serving the 6 phase of the first harmonic of the regular 
input, and it is equal] to 68(t,) for any instant of time 
within the t,St<t, 4, interval. 

Let us find some values of the statistical character- 
istic 68. By averaging (14) with respect to t,,, we obtain 
the following average value: 


a eae 
B= 2 Tyycorsyp as 


where oF = 2a, since, according to the condition g7=0, 


EE'=0 
By squaring and averaging (14), we shall find [6] 
the dispersion of 58: 
D (88) = 


pe oF ( cy oe of (Y sin Bo)* 
(Yo cos Be)? (Yq cos By)* 2(¥9cosBy)® ' 














Or a ae 











where uf = —(d°R/dr*), | 9 denotes the root-mean-square 


value of the frequency of the energy spectrum of fluctu- 
ations ¢. If condition (8) is satisfied, we can use only 
the first term in the above equation: 


a 
ree 


1 
D (88) = yi cosh” (16) 


The characteristic function 58 serves for determin- 
ing T [6] the distribution density w(58): 





w (88) = =|1 — ops aaa + 
— Cer ¥? cos 
(6 279 sing cos Bo| Yo 005 Bo , * = 
+ B 232 . Se Vin’ 


or, by taking into account the smallness of Or, 
38 ' * 2 


Y, cos cos By ->(= ) YQcos* hy 
w (5B) = (1 —_68 tg Bo) —* ‘- ; wane 


(17) 


For ERCS having a fixed quantization grid, the cal- 
culation of the sarne statistical characteristics of the re- 
lay actuation phase distribution 48 is considerably 
more complicated, Therefore, we shall calculate only 
the average value AB for this case. 

By analogy with (7) and on the basis of (3a), we ob- 
tain the initial equation in the following form: 


Y, (sin 8’ —sinB,) + 
+ [¥, sin oof; + & (ti)Imax t+ (42) —0 (18) 


After denoting s= Y sinBy and 9 =[Ypsinwot+ 
. &(t)] ax’ we shall use solution (14) for ¢ = &+ 9: 


"(E+ 0)? 


25"* 


AB = is 4+ E+06 +7) 


ws 3” 





Hence, for — = 0 and €*=0, and mutually indepen- 
dent values of —, &', n, and n‘, we obtain: 





cF) 


2? Toe)! (29) 


i.e., the problem of determining AB is reduced to the 
determination of 7 and 7)”, 
We shall determine 9 as 





n=s(r+or—s(r) +E), (ag 


where y corresponds to the extremum of g, and dy is 
determined from the condition 


s' (y+ or) + 8 (TE) =0. (21) 


Then 


n=s+O—s()=sOy), — (22) 


ny? = s 6 S as ;2 . 32 
n® = [s (Sy) + E(y + by)}° = s8* (67) + 32, (23) 


since s(5y) is a function of €", i.e. 
E. 


, it is independent of 
Thus, having determined, as before, the distribution 


density 6, from solution (21), we can calculate ® and 
and, consequently, also AB. In the final form, 


Ag = — sinfe _ (Se) 3g | 


y2 cos* Bo Wo 2¥2 cos Bo 
or 
mew (2) 4 
ae e— (=) 2Y2 Cox By (24) 


The determination of higher statistical characteris- 
tics AB involves considerable difficulties in calculation 
and requires the knowledge of higher €(t) moments. 


3. Determination of the ERCS Regulator 
Equivalent Gain in the Case of 
Parametric Random Action 











We shall determine the time t,, necessary for the 
self-oscillations to reach the steady-state in an ERCS 
with respect to the self-oscillation period T in the sys- 
tem, From the condition for the existence of self-oscil- 
lations in ERCS [5], it follows that T> T., =(20/w,,), 
where w,, is the maximum value of the self-oscillation 
frequency, which is determined by the equation 


20, (“ee + 0 (Gee ) => : (25) 


Here, 9, and 6, are the phase responses of linear 
factors in the equivalent diagram of the extremal sys- 
tem to be controlled before and after the nonlinear stat- 
ic characteristic, respectively. 

In the majority of cases of practical interest, the 
inertial portion of the system to be controlled can be 
represented by a connection in series of first-order fac- 
tors. In this, the transient process time can be approxi- 


t By taking only the first term in the expansion of the 
characteristic function ejvoB with respect to 68. 
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mately determined with respect to the interval 
where the real part of the transfer locus of the system 
inertial portion is positive: 


8, (Wp) -+- 3 (@p) = > (26) 


For any number of first-order factors, 
2x 


From a comparison between (25) and (26), it fol- 
lows that w,, is always smaller than Ws i.€., 


2x 2n ’ (27) 
t — a T. 
w <@ <0, < 


Thus, it can be said that self-oscillations in ERCS 
attain the steady-state in a time which is less than one 
period in all cases of practical interest, This fact makes 
it possible to determine, by using the describing function 
method, the equivalent complex gain S(Y, 58) of the 
control device as a function of 58(t,), assuming that 
58 = const during the averaging period t, , ,—t,. In this, 
S(Y, 68) will always represent a random step function 
of time instants t,. If we know the 58 distribution, we 
can determine the average value of the S(Y, 58) func- 
tion. 


For ERCS of the type under discussion, S(Y, 48) in 
exponential form for 68 #0 is expressed by [1] 


S(Y, 88) = akp ay ele 


or, for By= 1 — [1] and 56 = da, by 


i 


RY se 
The average value is given by 
NW, )=\re, dc.) w (8a) d (8a). 
—to ) 


After substituting in (29) the expressions for the in- 
tegrands from (28) and (17) and after integration (see 
{7]}), we obtain 


7 
ny op sin dy Yj cos* «, 
NW, 6) = fe oli 2 tale ee 


We shall expand exp (— a Yo cos” %) in a power 
series of (a/%¢ cos* 9) and, by using only the linear ap- 
proximation, we shall reduce (30) to the form: 

N(, 6a) = N(Y, 0) + bN(Y, 8a), (31) 
where N(Y, 0)=(— 1 %/4k,,,) e)™0 is the inverse of the 
equivalent gain in the absence of a random parametric 
action, and 6N(Y, 6a) is the average value of the cor- 
rection brought forth by fluctuations. 

By developing the N(Y, 0) and 6N(Y, 5c) expres- 
sions for = arc sin{(x -Y9)/¥J[1], after simple transfor- 
mations we finally obtain the following expression suit- 
able for graphic ——- 


N(Y,0) = — iz j_ [V 2xY. — x? + 7 (x — Y,)|, 
(32) 
ON (Y, da) = 
no} (x —Yo)* + ¥? cad (x — ¥o) | (33) 
V (2x¥o— x*)® V (2xY¥o— x?) 
The conditions for the existence of periodic solu- 
tions in the system are determined by the equation 








K (jo) = N(Y¥,0) + 8N(Y, 8a), (34) 
where K(jw) is the equivalent transfer locus of the extre- 
mal system to be controlled [1]. The graphic construc- 
tion in Fig. 2 gives an example of the fulfillment of 
condition (34), as a result of which the average values 
of the self-oscillation frequency and amplitude in the 
system are determined. 

As was indicated above, for ERCS with a fixed 
quantization grid, the derivation of the expression for 
the distribution density 48 is difficult, and the final 
solution of the problem cannot be obtained. We shall 
limit our analysis here to the following qualitative con- 
siderations. 

Consider the expression for the average deviation of 
the relay actuation phase, 


7 Eg 4 oF = oe 
B = Aa = da — ar i aa (35) 
z 
Let us introduce the R(T) = e /ayt" cos WT func- 
tion. 
Then, ma S ge 
t=—Ge) “+r eo 


where @, as is known, denotes the average value of the 
fluctuation energy spectrum frequency, and B denotes 
the average spectrum width. 

Ex xpression (35), together with (36), indicates a" 
Aa <éa always holds and that, for some values of wi, 




















Aa can change its sign. This means that for random ex- 
ternal disturbances applied to the regulator input of an 
ERCS having a fixed quantization grid, the self-oscilla- 
tion frequency for a given 
ing system and that, if Aa <0, it can be even lower 
than the self-oscillation frequency up of an undisturbed 


system. 
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Fig. 2. 


4. Certain Generalizations and Compara- 
tive Characteristics of the Statistical- 
Linearization Method under Consideration 











The statistical-linearization method which has been 
considered in application to ERCS can be even more 
readily applied to ordinary self-oscillation systems hav- 
ing relay characteristics of any type if, naturally, the 
restrictions imposed on €(t) are maintained. Since a 
random external disturbance €(t) applied to the relay 
component input in a self-oscillation system is always 
reduced to the random distribution of the relay-actuation 
phase, the equivalent gain S(Y, 5c) of a nonlinear ele- 
ment can be found as a step function with a random am- 
plitude distribution according to the usual rules of the 
describing-function method, 

The average value of S(Y, 5c) is determined once 
and for all with respect to the found expression for dis- 
tribution (17); this determination is reduced to the cal- 
culation of integrals of the type 


= —a?/ah 
\ eh x"dx == " ce : ) ’ 
da” (37) 
—0o 
which are calculated according to the method indicated 
in (7). 
For instance, let us find the S(Y, 5a) expression for 


a relay with the insensitivity zone 2a, As is known, in 
this case, 


4 
S(Y, 0) = => cosa, 
; (38) 


is lower than in the preced- 


where a= arcsin (a/Y), and the relay gain k. is assumed 
to be equal to unity for the sake of simplicity. 
Then 


S(Y, 6x) = + cos (a -+ da) 


and 


co 


S(Y, da) = \ S(¥, da) w(8a)d (da), (39) 


—o0 


where 


1 (2222 cost « 
w (da) = (4 —datga) = ee xr (o)e 
(see (17)], 


Let us consider some particular cases, 
1, Let a=0 (the relay has no insensitivity zone). 
For this case, we obtain 


S(Y, ba) = ~*, cos ba, 


v(do)— AL rely, 


(40) 
(41) 


i.e., we obtain a Gaussian distribution. 
By calculating the integral (39) for a=0, we ob- 
tain 


o i 
5Y, =e Mote, (4n 


where 8 =(¥//3o), 

2. For a#0, the w(5c) distribution is met- 
rical, leading to a complex expression for S(Y, 6c), 
Omitting the intermediate derivation steps, we shall 
write the result obtained incalculating integral (39) in 
the following form: 


S(Y, $a) = t cosae Ponte — 


ey 1 
— (tga +f) 5 Be BF costa (43) 


The second term in the right-hand side of Eq. (43) 
is due to the w(5a) terms which distort the symmetric 
6a distribution. 

Let us compare the obtained expressions for 
S(Y, 6c) with the results given in [3] for the same non- 
linearity types. 
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Fig. 3. 


Figure 3 shows the YS(Y, 6a) dependence plotted 
as a function of the 6 = ¥//2o parameter according to 
Eq. (42) and the curve, shown by the dotted line, plotted 
for the same case according to the equation given in 
[3]. The same diagram shows the YS(Y, 5) curves for 
a/Y =0.5 and a/Y = 1/72, which were obtained according 
to Eq. (43) without taking into account the distortion of 
the symmetry of w(5a), Similar dependences based on 
equations used in [3] are not provided here, since, for 
o->0, one. can use only a finite number of terms of 
the double infinite series,which expresses the equivalent 
gain value in this case in expanding the indeterminacies 
for their calculation. 


A comparison of characteristics for a relay without 
an insensitivity zone (Fig. 3) shows that both static-line- 
arization methods yield approximately equal results at 
least in the region of small values (in correspondence 
with the assumptions of the method proposed). The ex- 
pression S(Y, 6a) which was found here was given in un- 
expanded form and not in the form of an infinite series 
as in [3], thereby considerably reducing the amount of 
work involved in calculation. In conclusion, it should be 
noted that the statistical linearization method used in 
[3] cannot reveal the presence of a complex component 
in the expression for the equivalent gain for a relay with 
an insensitivity zone, which is important, for instance, 
in analyzing the stability of periodic regimes. 
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The present article describes a method of random scanning in application to problems in extremal 
control of multiparameter systems, A device which realizes this method is described. 


1. Formulation of the Problem 





In the majority of cases, the problems encountered 
in the synthesis of self-adjusting systems, which are 
currently gaining widespread use, are reduced to the 
optimization of a certain performance function Q of 
these systems, i.e., to the search for conditions which 
secure the minimum value Q® of this function. 

Let Q be a function of the parameters X,, Xp, ... 
Xn 


Q = QO (%, Ze, ..- Fn); (1) 
while Q® is the smallest of all possible values of this 
function: 


» » Sn) 
(— co << 4, << oo, i==1,..., A). 


(2) 


It is necessary to find values of the x4, x}, ..., oe 
parameters for which the Q function is minimized, i.e., 
values which satisfy expression (2), 

It is obvious that in order to solve this problem it 
is necessary to ascertain the existence of at least one 
minimum of the performance function, The number of 
minima is, however, not limited, It will be shown be- 
low that the proposed scanning system discriminates 
between the values of the minima and statistically 
tends to secure the smallest of these minima. 

At present, there are several automatic optimiza- 
tion methods, i.e., methods of searching for the optimizing 
series of parameters xf, .. Soitake From among these meth- 
ods, the method of alternating parameter changes 
(the Gauss- Zeidel method), the gradient method, and 
the method of quickest triggering (the Kantorovich me- 
thod) should be mentioned, These methods have been 
applied and are presently used in a large number of au- 
tomatic optimization systems (see, for instance, [1]). 


The scanning organized in correspondence with the 
above methods has the character of predetermined 
scanning, i.e., of scanning where the changes in the 
., X, parameters are rigorously predetermined 
by the type of function Q and the organization of the 
selected scanning method, 

However, predetermined scanning is not the only 
possible method, W. R, Ashby [2] came forward with 
his idea of random scanning, which was used as the ba- 
sis of the present paper. 


Xq, Xa, - 


2. The Random Scanning Method 





In terms of the control theory, Ashby's idea can be 
represented in the following manner, 

Assume that the system to be controlled can exist 
in a certain multiplicity set P={p} of states, finite or 
infinite. In this case, the system is considered to be 
adjusted if it is in one of the states contained in the 
multiplicity subset PCP, The only requirement impos- 
ed on this multiplicity subset is that it be of equal pow- 
er to the multiplicity set P, 

Thus, the problem of the adjustment of the system 
to be controlled is reduced to the selection of a method 
and organization of searching for one of the p>’ elements 
of the multiplicity set P among all possible elements 
p in the multiplicity set P, Ashby suggested that p* be 
sought not with respect to a given rule, as in the case 
of predetermined scanning, but by means of random se- 
lections to be performed by a “selection amplifier" 
(3). 

Ashby's selection amplifier operates in the follow- 
ing manner (see Fig. 1), A certain state p, selected at 
random, is communicated tothe system to be controlled 
SC by means of a generator of random events, GRE, 
which secures a regime of random selections from the 
source of all possible states SS, i.e., from the multipli- 
city set P, In this, the selector unit SU, after receiving 
information on the reaction of the controlled system to 
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the variant supplied by GRE, either sends a command 
to GRE to produce the next random selection in the 
case of p#p° or, in the case of p=p , disconnects GRE 
and thereby locks the system in this state. 

It is obvious that, in operation with respect to such 
commands, the selection unit will screen out the p#p" 
states until it finds, by chance, the required state, 
which it will preserve continuously until any external 
or internal factors transfer the system from this state 
into another p#p* state. After this, the above-con- 
sidered mode of random searching for the state p* will 
again become operative. W.R. Ashby's well-known 
homeostat operates on the basis of this principle [2 and 
3). 

It should be noted that Ashby's system is character- 
ized by three specific features; it has a source of all 
possibilities (the multiplicity set P of elements), a de- 
vice securing a regime of random selections from a 
general totality P (generator of random events), and a 
device which supplies the criterion according to which 
the required element is selected in the system to be 
controlled (selection unit), Another important feature 
is the presence of two sharply different operating con- 
ditions; the temporary, transient random scanning mode 
(for PtP ) and the other, the stable adjusted mode (for 
p# p*),whichwas to be found. It can readily be seen 
the the described method is essentially the realization 
of the trial and error method, which is widely used in 
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Fig. 1. 


3. Application of Random Scanning 
to Optimization Problems 








In order to apply the described random scanning to 
the problem of minimizing the Q(x, Xz, ..., x) func- 
tion, the “state” p of the “system to be controlled" Q 
must denote the velocities of its motion along the co- 
ordinates, i.e., p=(x1,%2,. . .,X,), father than its co- 
ordinates in a parametric coordinate system. It is obvi- 
ous thax p is the velocity vector of a point representing 
the state of the system to be controlled in an n-dimen- 
sional parametric space {x,, Xz, ..., x,,}. Under these 
conditions, the control system depicted in Fig. 1 is ap- 
plicable. The above-mentioned operating modes of the 
system are also preserved with the exception that the 
system may not be immediately transferred to the Q* 
state under the adjusted conditions (p= }* ) for n> 1; 
however, in this case, it will “approach" the required 
state, i.e., the performance function Q will have a 


smaller value at the end of the adjustment conditions 
duration than at the beginning. 

As the most natural criterion for the selection of 
states supplied at random, we can use, for instance, the 
value of the derivative 


dQ _ aE aQ - (3) 


If all possible states p =(X;,X,. . ., Xn) which the state 
source SS can supply pertain to the multiplicity set P, 
the multiplicity set P*must contain those states of the 
system to be controlled which satisfy the inequality 


<0, (4) 


e., the condition for the diminution of the performance 
function, This inequality represents the selection criter- 
ion, It can readily be seen that the probability that the 
velocity p satisfying this criterion will be selected as 
close to }, since the “quantity” p* € P* is slightly less 
than one half of all p € P if the hypersurfaces Q(x;, x2, 

. + «» Xp) = const in an n-dimensional parametric space 
are sufficiently smooth, 

Thus, the algorithm of random scanning in the given 
case consists in the fact that the p states, which lead to 
(dQ/dt)=0, are transient and are immediately replaced 
by other states also chosen at random. The states lead- 
ing to the inequality (dQ/dt)< 0 are stabilized until this 
inequality changes its sign due to the change in the de- 
rivative. It is obvious from (3) that the change in the 
derivative in the fixed state, i.e., for P= const, occurs 
due to changes in partial derivatives (0Q/0x;)(i=1, 2, 

.+» 0) during the motion of the point representing the 
system to be controlled in t! parametric space. The 
transition from the adjustment value of f to another val- 
ue takes place at the moment at which the derivative 
changes its sign, i.e., when (dQ/dt)=0. It can readily be 
seen that, as a result of operation with respect to the in- 
dicated algorithm, the quantity Q can be reduced to its 
minimum value Q” after a certain finite number of cy- 
cles in selecting p. 

However, the realization of such an algorithm re- 
quires the calculation of the (dQ/dt) derivative in corre- 
spondence with the selection criterion (4), which, in 
the endreduces the noise stability of this control system, 
especially for small values of the derivative modulus, 
Therefore, we chose a different method that did not re- 
quire the calculation of the derivative and was simple 
and reliable, although it did reduce the response quick- 
ness in adjustment to a certain extent. 

We proposed an arrangement whereby the Q value of 
the performance function is continuously compared with 
a certain level n, which changes at a, constant rate q< 0. 
Let the values of Q and q coincide at the initial moment. 








Then, 9 (t)=Q; + qt, where Q is the value of the Q func- 
tion at the beginning of the ith stage, and ¢ is the time. 
Let us denote by Q(p, t) the performance function 
value at the instant of time t for the fixed rate p of 
changes in the system parameters. If, for the rate value 
at the beginning of the ith stage and for sufficiently 
small values of t> 0 the inequality 


Q (pm, t) —n(t) <0, (5) 


holds, the Q; , , value of the performance function at the 
instant of oun’ t= te when the 9(t) and Q(Py, t) values 


are compared again, is smaller than the Q value at the 
initial moment. At the instant of time t,, the rate val- 
ue is switched from p, to the B, ,, value, and a new stage 
in the process is initiated. 

Consequently, if inequality (5) is satisfied after the 
beginning of the ith stage, the performance function 
value is smaller at the beginning of the next stage. 

Thus, inequality (5) can serve as the selection cri- 
terion, Then, all values of 6 satisfying this inequality 
pertain to the multiplicity subset >. The p values that 
do not satisfy (5) are immediately rejected. 

Criterion (5) does not require the calculation of the 
derivative, thus considerably simplifying its realization; 
however, it is somewhat narrower than criterion (4). Ac- 
tually, inequality (5) for t= 0 can be written in the dif- 
ferential form: 


aQ 
<4. (6) 


The obtained inequality determines the selection of 
rate p at the initial moment, The smaller the q value, 
the smaller the domain of § values which satisfy this in- 
equality, Therefore, by comparing (4) and (6), we as- 
certain that since q< 0, the domain of rate values satis- 
fying (5) is always smaller than the domain of rate val- 
ues satisfying (4). 

Figure 2 shows the variation of Q for two values of 
B;. The thin lines in the same figure indicate 9 levels 
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for different values hi G2> qs. It can readily be seen 
that, for q=qg, the p ) value does not belong to P*, 
while Ve €P* in an a interval, With an increase in q, 


the multiplicity set eo becomes wider, Thus, for q=q;, 
both rate values pos? ep”, 
Therefore, we can consider that 
P* = P*(q) = 


has been established. 

For q=0, inequalities (4) and (6) coincide, and the 
selection of rates at the initial momentwill be equal for 
both criteria. However, at the switching moment, cri- 
terion (6) will yield Q, | ,=Q;, i.€., there will be no ad- 
justment. For large q< 0, the ee will be weak, 
since it follows from (5) that the Q;-Q; , = 74 differ- 
ence is small, regardless of the fact that te selected 
rate led to Q< Q Pe “se during the operation process 
(see the q= as case for 5° » in Fig. 2). 

Thus, de is aaie wis an increase in q leads to the 
enlargement of the multiplicity subset P*and, on the 
other hand, lowers the adjustment efficiency. Therefore, 
there exists such a q=q° value which optimizes the con- 
sidered control system with respect to adjustment speed, 

The 9 level is essentially one of the memory or- 
ganization forms where 9 > 0 is “memorized” instead of 
Q. In this case, the “memory” is shifted in time inde- 
pendently of the Q value until the moment at which 
Q=. The presence of the 9—Q difference, the sign of 
which determines the device operating conditions, con- 
siderably improves its noise stability. 

It should be noted that this memory form is appli- 
cable also in predetermined scanning, which is organized 
according to the method of alternating parameter 
changes or according to the quickest triggering method. 
In this case, there is no necessity of determining tue 
minimum of the performance function Q in the selected 
direction with respect to criterion (4), and it is suffi- 
cient to introduce the described “memory” by means of 
the comparison signal. Such an improvement, as was 
shown above, considerably improves the control system 
noise stability, although it somewhat lowers the scan- 
ning rate. 

Figure 3 shows the block diagram of a device cap- 
able of realizing extremal control by the described me- 
thod of random scanning. 

In this diagram, the GRE generators secure a re- 
gime of random selections from the SS source of all 
possible states and act on the rate of change in the pa- 
rameters of the system to be controlled, SC. In this, 
the changing quantity Q is continuously compared with 
a certain signal 9 generated by the generator of com- 
parison signals, GCS. The result of the Q2 y compari- 
son determines the system operating mode, which is 
assigned by the control unit, CU. The operating modes 
are given in the table. 


























TABLE 1. 
Operating modes | Cc Q n 
preparation | Q> const e>0 
scanning | Q<n var q<0 











During the preparation stage,GRE's are disconnect- 
ed, and the controlled system is not excited, i.e., the 
tracing point is immobile and Q=const. For Q=9, the 
preparation operating mode is replaced by the scanning 
mode and during this stage the system parameters change 
at the constant rate f, the direction and modulus of 
which are determined by GRE at random, At the mo- 
ment when Q= 9, the scanning mode is replaced by the 
preparation mode, etc, 

It can readily be shown that the indicated device 
corresponds to the above-considered algorithm of ran- 
dom scanning with the selection criterion (5). Here, 
the q-parameter has a rigorously defined physical mean- 
ing: it is the rate of change in the comparison signal 9 
under scanning conditions. It is obvious that, in the ab- 
sence of hysteresis in the considered control system, the 
operation is possible only for 


Onin < q. (8) 


In the opposite case, the Q<q inequality would 
never be satisfied. 


4. Effects of Hysteresis 


It has been shown that the described device operates 
under two regimes: “scanning” and “preparation,” and 
that the moment of transition from one regime to an- 
other is determined by the coincidence of Q and 9 sig- 
nals, However, in actual devices, due to a number of 
causes (inertia, hysteresis of the comparison system, 
etc.), the switching of regimes takes place for different 
relations between 9 andQ. Thus, the moment of transi - 
tion from the preparation regime to the scanning regime 
or the reverse is determined by the equations, 
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respectively, where 6, and 4, are usually sufficiently 
small positive quantities; 5; = 5(c) and 5,=5,(q, p). 

It can readily be seen that, for disconnected GRE, 
i.e., for Q=const, self-oscillation conditions will arise 
sawtooth oscillations of the 9 signal with respect to Q, 
where 99x ~ "min ~ °1* 52 


These self-oscillations, a consequence of hysteresis 
in actual control systems, introduce interesting peculiar- 
ities in the device operation. 

Thus, in the presence of self-oscillations, control is 
possible even if inequality (8) is not satisfied, i.e., if 
lp am <|q|. This is obvious from the graph in Fig. 4, 
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where two scanning cycles for the favorable (T, interval) 
and the unfavorable (T; interval) variants of the Q func- 
tion variation (solid line) are shown. As a result of the 
action of both variants, the peformance function is re- 
duced by h due to the fact that, in adjustment, the mu- 
tual rate of change in © and 7 is less than in the unad- 
justed state. This leads to a change in the duration of 
the above regimes: T,> T. 

Thus, as a result of hysteresis for a large |q| , the n 
signal follows the Q quantity under self-oscillation con- 
ditions, while the duration of favorable variants of Q 
variation (O< 0) is somewhat longer than that of unfavor- 
able variants (Q> 0) which ultimately leads to a sys- 
tematic reduction of the performance function. There- 
fore, in certain cases, the deliberate introduction of hys- 
teresis into the control system is entirely possible. 

Another peculiarity due to hysteresis is the possibili- 
ty of temporary derangement of the system to be con- 
trolled. This is due to the fact that a certain time is re- 
quired for the “trial” of the random variant offered by a 
GRE, If this variant is unfavorable (Fig. 4), the Q func- 
tion will be able to change during time T;. If several 
such variants (b# B*) are supplied, which is possible al- 
though not probable, the system to be controlled will be 
deranged to a certain extent. 

The possibility of such derangement renders the de- 
scribed control system selective with respect to different 
minima of the Q function if such are available. Actual- 
ly, after Q attained the intermediate minimum B (Fig. 








5), the system will continue its motion in the neighbor- 
hood of the point B until the system, in its random mo- 
tion due to the possibility of derangement, reaches the 
point C, From this point, the system can either return in 
the opposite direction or, due to the minimizing action 
of the described device, reach the neighborhood of the 
point A, In the first case, the whole process will be re- 
peated from the beginning until the tracing point reaches 
point A —the smallest minimum of Q. Naturally, motion 











Fig. 5. 


in the opposite direction, from A to B, is also possible; 
however, derangement variants which occur considerably 
less often are required in this case and, consequently, a 
longer time for overcoming the h,> hy barrier will be re- 
quired, Consequently, the probability of a transition 
from A to B is smaller than the probability of transition 
from B to A, Thus, the system, being able to dwell in 
any region, will remain for a conside-ably longer time 
in the neighborhood of its minimum value Q than in 
any other region. 


5. Experimental Verification 





The possibility of optimization by the random scan- 
ning method according to the proposed scheme was 
checked experimentally under laboratory conditions in 
application to simple self-adjusting systems. 

We considered the problem of the self- balancing of 
two independent small beam-weights which formed a 
system having two degrees of freedom. The beams were 
maintained in the horizontal position by means of leaf 
springs whose tension determined the deviations x, and 
X» from the equilibrium states of these beams. The ad- 
justment of such a system consisted in finding the equi- 
librium conditions for which the deviations x, and x, of 
these beams were simultaneously equal to zero; x4=x=0. 
The adjustment of this system was performed by means 
of mobile weight trolleys, which were driven by two 
electric motors, respectively. The performance function 
Q of the system was realized in two variants; Q= |xg+xJ 
and Q=max Ix,| (i=1, 2). 

Thus, the problem was reduced to the minimization 
of the Q function, which attained its minimum (Q"= 0) 
only when the system was adjusted, i.e., at the moment 
at which both beams were in equilibrium. 


The comparison signal generator consisted of a 
voltage divider ~a rheostat which was electromechan- 
ically driven. In the first case, the selection unit was 
based on a three-coil polarized RP-5 relay, where |x,| 
and |x,| were fed to two windings, and the comparison 
signal n was fed to the third winding. The second vari- 
ant was based on two two-coil polarized relays, The 
random event generator consisted of an ordinary PS-64 
scaling device, to the input of which noise from a tube 
was supplied; when the noise was disconnected, the 
trigger cells were fixed in a random position which de- 
termined the motion of the weight trolleys that balanced 
the beams. 

The performed experiments showed that: 

1, Adjustment, i.e., minimization of the Q func- 
tion in both variants, is possible through use of the de- 
scribed method, 

2. The adjustment itself is independent of the rate 
q of signal variation 1 in a wide range; for different val- 
ues of q, only the average adjustment time changed. 

3. There exists an optimum value q=q° for which 
the average time for the system adjustment is at the 

~minimum, all other conditions being equal, i.e., there 
exists a q for which the described adjustment system is 
optimized with respect to quickness of response. 

Another experiment was prepared for verifying the 
possibility of adjustment by the random scanning me- 
thod in the presence of a systematic disturbance factor. 
For this purpose we used a beam with two trolleys, one 
of which moved at a constant speed, thereby introduc- 
ing the disturbance, while the other was included in the 
above-described control system for the purpose of ad- 
justing (balancing) this single-parameter system. 

The performance function was realized in the nat- 
ural way: 


where x is the beam deviation from the equilibrium po- 
sition, Such a device is essentially the simplest statis- 
tical follow-up system operating on the basis of the ran- 
dom scanning method, 

Experiments showed that the performance of such a 
self-adjusting system is largely determined by the pro- 
cess speed and, consequently, by the selection of the q 
parameter. 


SUMMARY 

The proposed extremal control system, based on 
random scanning, can be successfully applied in auto- 
matic optimization of multiparameter systems to be 
controlled for any performance function, including dis- 
continuous functions. The advantages of this system 
are particularly pronounced in the case of control with 
respect to a large number of independent or dependent 
parameters of the system to be controlled. Computers, 
which are necessary in the case of predetermined scan- 
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This article presents a study of the second B, V. Bulgakov problem where the nonlinearity of 
one of the measuring units is taken into account and a method is described for solving the 
stability problem, which makes it possible, on the one hand, to find acceptable regulator 
parameters and, on the other hand, to determine the controlled variables range, containing 
the coordinate origin, for which the system stability is guaranteed. 


1. Formulation of the Problem 





In their paper, A. M. Letov and A, I. Lure [1] for- 
mulated one of the most important problems in automat- 
ic control; the determination of stability criteria for 
controlled systems described by equations of the follow- 
ing form: 


Me = S) Oma Ma + mE, 
«=i (1) 
5=/(5), 5= >) Pa (Ma) — (6), - 


where p,,(1,,) and 9(€) are functions of the (A) class. 


This paper also treats the effect of various parameters 
characterizing nonlinear dependences on stability of sys- 
tems to be controlled. The problem stated is of great 
practical importance, Actually, let us consider an air- 
craft controlled with respect to yawing by a nonlinear 
autopilot having rigid feedback. This autopilot has a 
damping gyroscope which can be described by the follow- 
ing equations: 


72x 4+ Ua + Ka = (6), 
3 = Mp+ N%, 
where & is the deviation angle of the damping gyroscope 


and y is the aircraft yaw angle. The %(9) function (see 
Fig. 1) is of the following form: 
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kp for > Go, 


$0(5)={ ko for |6]< do, (2) 
— kB for O< -- Gp. 
The T and U, coefficients are usually smal! and, in 


the case where the %(9) function is linear, they can be 
neglected [2]. Then, the equations of the aircraft motion 
with respect to the yaw angle can be approximately rep- 
resented in the following form: 


T+ Up + Kp+p=0, 
b=f' (0), c= ap+ Ba— 5 B, (3) 
a= @" (5), 6= Mp + N%. 
Here, ¥ is the yaw angle, ft is the rudder deviation 
angle, g*(9)=(1/K,) %, (2) and f*(0) is the servomotor 


characteristic (see Fig. 2), which has the following pro- 
perties; 
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The problem of the stability of system (3) in the 
case where the ¢°(G) function is linear has been thor- 
oughly investigated in [2]. The following question now 
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arises; will the stability of system (3) deteriorate if the 
¢*(@) function is nonlinear or will it be possible to find 
such autopilot parameters for which this nonlinearity will 
not impair the stability of system (3) in comparison with 
the case wherein the gyroscope characteristic is linear. 
The purpose of the present paper is to determine the ef- 
fect of the nonlinearity of g°(G) on the aircraft flight 
stability and to find the stability criteria for the system 
(3). It is obvious that system (3) is a particular case of 
system (1). 


2. General Form of the 
Stability Conditions 








In investigeting the stability of system (3), we shall 
use the second (direct) Lyapunov method. The Lyapunov 
V function is usually plotted in quadratic form, and the 
coefficients in the expression for V do not depend on the 
system coordinates. However, as is known, the latter 
condition is not compulsory. It is sometimes difficult to 
plot the Lyapunov function in such a way that the coeffi- 
cients of the V form are independent of the system coor- 
dinates. Actually, if some coefficients of the V form 
depend on the system coordinates, the problem consists 
in finding a continuous continuum of these coordinates, 
which contain the coordinate origin and which lie within 
a certain interval (in the control range), and in selecting 
the remaining coefficients in such a way that the V func- 
tion is actually a Lyapunov function for the given sys- 
tem. This method has been chosen here for solving the 
above problem.* 

Before we start plotting any diagrams, let us reduce 
the initial equations (3) to the normal and dimensionless 
forms, 

By introducing the following notation: 


dip 





m= ¥, m=TT, ¢= Tr, b,, = — K, 
U 
ba=— FF) 
N°K 
u=—i1, p=a, u=-—- jr 
(5) 
N2U\ 1 
2=(M—ar) FF: 
N?2 


p(s) = Eg*(s), 


i ie 
—e 8 r=, 


/(s) = Tf (9), 


= 


we obtain: 
mh = Ne, 
Na = Beit +- Doo e + Usp, 
p= 7 (9), (6) 


o= Pm + P(s)— 7H, 
3 = QM + 92% — TH. 


It is obvious from Fig. 2 that, in the control range 
|a|=<|o'|, the entire f(o) function is located between 
two straight lines y=h,o and y=h,o, where h}>h3. In 
this interval, h(o)o can be substituted for f(o), where 
h(9)=(1/9)f(9), 


hy = Th, >h(s) > Th; = ho, 
h (0) = lim + f (3) = hy = Thy. 
o-0 


Similarly, in the |G|= O° range, ¢(G) can be re- 
placed by k(G) 6, where 


k (6) = k= Egt>k (s) > a = ke. 
Thus, system (6) can be written in a different form: 
wh = Ne» 
Ns = bet + Sea%2 + Wop, 
p=h(s)[(m+akO)m+ (0 


+ gak (5) ns — (r + rk (5) p]- 


We assume that system (6) remains stable after the 
automatic device is switched off, i.e., that ail roots of 
the r*—beed —bg, = 0 equation have negative real parts. 
According to the well-known Lyapunov theorem, it is 
always possible to find a positive function 


F= § (Bunt + 2B + Bini) 


whose sign is determined and which satisfies the equa- 
tion 


OF OF 
om N2 + om (be:M1 4+- bo2N2) = (8) 


=W = — (Ayn? + 24am M2 +. Aoan}), 


where W is an assigned negative function whose sign is 
determined, i.e., the Ay, Ayp, and Age coefficients, 
which are not yet known, satisfy the conditions 


Ai, >0, Ay;4g — A? > 0. (9) 


From Eq. (8), the expressions for the By, By2, and 
Bez coefficients can readily be found in terms of the Ay, 
Aya, and Ag, coefficients: 


| 





By = wad [(bex — be.) Ay + 2ba1b22A 12 -. b2, Age|) 


(10) 


A At; — 521 Acs 
Ba=—7 2 





* The general idea of this investigation method has 
been presented in a paper by V. I. Zubov [3]. 








Let us now consider the following positive function 
whose sign is determined: 


i 
V=F(m, %)+>P- 


According to Eqs. (7), we find: : 
—V=—W— U2B Mp —UsBoonop - 
— h(3) [pr + ak (5)] ae — 
— qh (s) k (6) now + h(o) lr + rk (s)| y?. 


According to the Sylvester criterion, V will be a 
negative function with determined sign if the following 
inequality is satisfied: 


Ay, Ais P, 
D — Ais Ags P, > 0, (11) 
P, P, h(s)[r+7k(a)] 


Py = [wpBiy + h (0) (pi + rk ())], 


P, = + [ugBes +h (0) k (3) |. 

According to the well-known Lyapunov theorem, 
system (6) will be asymptotically stable if inequalities 
(9) and (11) are satisfied. Inequality (11) must be satis- 
fied for continuous continua of the 0 and @ coordinates 
which fill the |o|<0* and |G|=G* regions, It is obvi- 
ous from Figs. 1 and 2 that the o°and S* values corres- 
pond to the lower branches of the straight lines y=h,o 
and y=k,@ of the nonlinear functions, Therefore, we 
shall henceforth consider only domains where h and k 
change. The extent of these domains will determine 
the character of the system stability. This extent obvi- 
ously depends on the autopilot parameters. 


where 


3. Investigation of Stability Conditions 





We shall now pass to the investigatior-of stability 
conditions for system (6), In the particular given case, 
we assume that Aj.=0 and Ay;=Ag.> 0 for the sake of 
simplicity. Inequality (9) will then always be satisfied. 
By expanding the determinant D, we obtain: 


Aj,h (r + rk) — 4s {usBy. + h(py + qik)|? - 


(12) 
- om [uBos “++ qahk]* > 0. 
We shall denote: 
ghk=2z, ph+qhk=y. (13) 


By substituting (13) in (12), we obtain: 


( + wgBoo)® + (y + UgBy2)* < 4Aggh(r + rk). (14) 


We can readily find the expression 


he (7. 7h) ment geet ( 2 as 288 be oe 
(r+ rk) = ey + (2 — 2) ea Ay + Bs, 


where 


y payee pam rn 
Pi’ P1942 , 


Inequality (14) can be written in the following 
form: 
[zx + (U2Bo2 — 2AqoB)P + 
+ ly + (ueBis — 2Aa2A))* < 443, (A® + B*) — 


— 4A,.Bu,Bo, — 4A, AUB >. 
(15) 


Let us substitute the expression (10) for By, and Bay 
in (15) and represent it in the form; 
(X —a)*+(¥—y)? <4R%, (16) 


where 


ee Drags 9 = 2p — w(t — bn) 
X= ’ Y= , Zo = 2B berbes , 


Y= 2A+ 52, R= A+ B+, 
en on a ee 
Ce us| barbea ‘a: 


It is obvious that X and Y are not constant quantities 
in inequality (16), 

It can readily be shown that, in order that inequality 
(16) be satisfied, it is sufficient that the following in- 
equalities hold (see Fig. 3): 


(%9—V2R)<X <(%+V2R), 
(Yyo—-V2R)<Y¥ <(yot+V2R), 
or 


(io —-V 2 R) Age < qahk < (tq + V2 R) Ags, (17a) 
y 

















Fig. 3. 








(Yo — V2 R) Ars < pik + guhk < (yy + V 2K) Age. 
(17b) 
In accordance with inequalities (17a) and (17b), we 
can find an allowable domain where h and k change and 
determine the effect of nonlinearity on the stability of 
system (6). Before this, according to the initial transfor- 
mations of the coefficients, we find the expressions for 
R?, Xo, and yp: 


pe — Tt (Nial + NeK? 
= Gr (MT? —NUy 





+ T? (N*al + N*K) al+ K 














alKU (MT? — N2U) aeK ’ 
_ 2T (Nal + N*K) , T(1+K) wt Bra ag § 
% = —(MT?=— Nw) ~~ Ko * = ai +x: 
Denoting 
_ T (Neal + N2K) nthe: _(+K)T 
1=—yn—wy f=: S=-T—- 


(18) 

Equations (18) include the autopilot coefficients. In 

the given case, it is obvious that §=0 and s=0, and the 
expressions for R*, xg, and yp assume the following form: 


RE =: 2 (Bn? + Stn + H+ 8), 
(19) 
: m= (2En +S), yo == (1 +28). 


Let us find the conditions under which the inequality 

(12) or (16) is satisfied. According to inequality (12), it 
is necessary that the h=m > 0 relation holds where m is 

a certain finite number. Due to property (4) of the f(o) 

function, the latter inequality can always be satisfied. 
The R* > 0 condition, which is necessary for inequality 
(16) to be positive, follows from inequality (16). Let us 
plot the R?=0 graphs in the (n, €) plane and separate the 
domain where R*> 0 (Fig. 4). 
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H It is obvious that the R’> 0 condition is not sufficient 
H for satisfying inequality (16) and that inequalities (17) 

io must be considered in detail. Two cases, determined by 
i the sign in front of q, ~which coincides with the sign in 

i front of n— are possible here. 


v 
oe 
i 





1. In the case where q2=(1/T*)(MT?—N2U) = 0, in- 
equality (17a) yields the condition: 

Zo os V2 R > hyk, 

Ly — V2 R heka 





= a8 >i, (20) 


where the quantities «=(h,/h,) and 8 =(k;/k,) character- 
ize the domain where h and k change. 


Let us write inequality (20) in the form 


PLy — V2R <0, (21) 
where p = on <1. 


This inequality determines certain stability regions 
in the (9, &) plane of the autopilot parameters. The 
equation for the boundary of this region is pxy— /2R=0. 
After substituting expression (19), this equation assumes 
the following form: 


(4p* — 2) En? + 
(22) 
+ (4p? — 2) SEn — 26° — 26 + p*S?=0. ° 


The boundaries of the stability region for different 
values of p are shown in the right (, €) half-plane in 
Fig. 5. 
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It is obvious that inequality (17a) can determine 
only the domain of the h and k quantities. Their maxi- 
mum values are determined by inequality (17b). 

Actually, inequality (17b) yields the conditions 


[o —V2R)—#(@—V2R) |< 


<h<[m+V2R)—2(,+V2R)|=, 


Pi(w—V2R) wy, 
a2 —V2R)—nlo—V2R) ~ 


< 





Pr (Zo +V2R) alia 
qa (Yo + V2 RB) — 41 (&> + V2 A) 











The values of Ag, and k have not yet been deter- 
mined; suitable values of these quantities must be se- 
lected in order that the last two inequalities are satis- 
fied. 

2. In the case where qd» < 0, X9(m, €) in the left 
(n, €) half-plane has different signs for different values 
of n and €, as is shown in Fig. 4. We have two cases, 
respectively: 

for X9> 0, inequality (17a) yields the condition 


r—V2Rh<0, (23) 
for X9< 0, inequality (17a) yields the condition 
—prz,—-V2R< 0, (24) 
where p = rt < 1. 


It can readily be shown that the equation for the 
boundary of the stability region, determined by inequal- 
ities (24), is of the same form as equation (22). The 
boundaries of the stability region for different values of 
p are shown in the left (m7, &) half-plane in Fig. 5. 

As in the preceding case, the value of k will be 
similarly determined by the inequality 





(v9 + V2R) pr . 

0 —e 

q3(Yyo—-V2R)—qle+V2R) ~~ 
< (x9 -- V2 R) rs 





~ ga(vo + V2 R)—a1 (te —V2 A) 


It is obvious from Fig. 5 that the. stability region 
with respect to the 9 and € parameters is restricted as 
the value of p increases, i.e., as the domain wherein 
h and k change becomes larger. Therefore, generally 
speaking, the introduction of nonlinearity ¢(6), which 
characterizes the gyroscope nonlinear property, reduces 
the stability of system (6). However, if, for a certain 
nonlinear function ¢(G), corresponding, for instance, to 
P =P, 9= py =1, we choose the q and € parameters to 
be in the p =p, region, the nonlinearity will not affect 
the system stability. In this, in the region corresponding 
to p = 1, we can make the domain of k variations for 
which the system stability is guaranteed as large as we 
wish. 

In conclusion, the author extends his thanks to A. 
M. Letov for his valuable advice and remarks. 
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An exact method for determining periodic regimes in arbitrary sectionally linear systems is 
considered. The solution of the problem is reduced to the determination of forced oscilla- 
tions of a linear system with periodically changing parameters, These forced oscillations 
are found in an exact manner by means of a new approach, which is based on the solution 
of an auxiliary Fredholm integral equation of the second kind. The thus obtained periodic 
solution depends on the unknown times of motion along sections of nonlinear characteristics. 
By using this solution and the conditions for the transition from one section of the nonlinear 
characteristic to another, it is possible to determine the transcendental equations of periods. 





1, Formulation of the Problem 

By a sectionally linear control system we shall des- 
ignate a system which, for certain given time intervals, 
is described by different systems of linear differential 
equations with constant or variable coefficients. In the first 
case, such a system will be called a steady-state system, and 
in the second case, it will be called an unsteady -state system. 

The simplest example of sectionally linear systems 
is those which differ from the linear systems by the pres- 
ence of nonlinear elements having sectionally linear 
characteristics, Various spontaneously changing cyber- 
negic systems can be considered as the more complex of 
this class [1]. Until now, only steady-state sectionally 

linear systems were generally considered in the literature [2- 
5). This problem was first studied in [2] in which a me- 
thod was presented for the determination of periodic re- 
gimes by fitting the linear-system solutions corresponding 
to different time intervals, In so doing, the solution of 
the problem was always brought to the point at which a 
certain transcendental algebraic system of equations was 
composed with respect to unknown times of motion for 
different sections of the nonlinear characteristic ("equa- 
tions of periods”). 

The practical application of the fitting method in the 
form proposed in [2] often involves difficulties due to the 
following reasons; 

1. The period equations depend on the roots of char- 
acteristic equations of the corresponding linear systems 
which, if the considered system of differentia! equations 
is of a high order, precludes the possibility of expressing 
the obtained results in terms of initial parameters, 
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2, The amount of calculation involved in composi- 
tion of the period equations increases markedly with an 
increase in the order of the initial system of differential 
equations, 

3. This method does not make it possible to take in- 
to account the filtering properties of the linear portion of 
the system and to determine only a finite number of har- 
monics in the Fourier expansion of the periodic regime to 
be found, 

A modified fitting method for systems with a single 
nonlinearity, making it posssible to write the equations of 
periods in terms of the initial system coefficients, was 
proposed in [4]. In [3], a method different from those de- 
scribed in [2) and [4] was used in composing the equations 
of periods by fitting for the case where the system con- 
tained only one nonlinearity, the characteristic of which 
consisted of straight lines parallel to two given lines, 

The method used in [3] was based on a special affine 
transformation of the nonlinear characteristic and could 
not be extended to nonlinear characteristics of a more 
general type. 

A linear transformation making it possible to reduce 
the amount of calculation involved in the application of 
the above methods was proposed in [5]. 


* The basic results of the present paper were reported at 
the Seminar on Nonlinear Problems in Automatic Con- 
trol Theory, held at the Institute of Automation and Re- 
mote Control, Academy of Sciences, USSR, on December 
2, 1959. 








In the present paper, the solution of this problem is 
reduced to the consideration of the Fredholm integral 
equation of the second kind, providing the possibility of 
composing the equations of periods for an arbitrary sec- 
tionally linear system. An approximate method, making 
it possible to use only a finite number of harmonics of 
the periodic regime to be found, is also described. In 
this, all results are expressed directly in terms of the co- 
efficients of the initial system of differential equations. 


2. Transmission of Periodic Signals 
through the Linear Portion of the System 








Consider the differential equations 


ty = >) dyata + F,(t) (k= 1.2... 


a=] 


»n), (1) 


where b, are constants, and F (t) are periodic functions 
of the period T, which we shall consider to be sectional- 
ly continuous, It is necessary to find the periodic solu- 
tions of the system (1). This problem has been studied in 
detail in [6] and [7]. Another solution form, necessary 
for further derivations, will be given below. 
Let 

D (d) = | Oia — Myc | (2) 
be the characteristic determinant of the system linear 
portion,t and let A), be the cofactor of the element of 
determinant (2) whith is in the kth row and in the ath 
column, Subsequently, it is sufficient to assume that the 
equation 


D(h) = 0, (3) 


the roots of which are A = = (p=1,2,. 
have a root equal to 


-+, 1), does not 


w= es, (4) 


where gs is zero or any integer, 


Under the above assumptions, the periodic solution 
of system (1) is uniquely defined by the equation 


a (t)= 


of 
= xh Of, (t — 1) Fx (x) dt (s = 1, 2,...,n), (5) 
k=10 
where 
1s is (Ym) 53 


Functions (6) can be represented in finite form. If 
Hq and for p #a, we have 


Ars (Ao) ee 
— Day) ee? —1 (1) 


For y=0 and y=T, the sum of the series (6) is equal 
ji, (0) = 





OF, (y) = Dis (y) = 


to 


O12 (T) = 


1 = Ags - eth —— ‘e 


> 0O<y<T. 
: for 
Moreover, 
Dx; (0) - Dxs (T) = 8x5. (9) 


Consequently, for k#s, the a) (tr) function is 
continuous with the square given by 


O<t<T, 0 tT, (10) 


and the ar (t—T) functions have a discontinuity on the 
t=T diagonal. 

In using expression (7) in (5), it must be taken into 
account that, for 0<t=<T, the argument t—r of the in- 
tegrand in (5) varies within the -T<t—rt <T interval 
and, in correspondence with (7), 





OL(t—nD= 
x Ans (Ap) eet , 
’ = 7 for 0 < t —— Tt & 7 ’ 
2 D(A) eet 4 


n tT 
s} Aes (Ae) eT pe!) 
D(a) _'o* 3 








fo O<t—t<T. 


(11) 
If Eq. (3) has the roots A,,..., An (the multiplicities 
of which are v4, vg,...,¥,,), we shall have the following 
equation instead of (7): 


\ e=l1 e 











Wi; (y) = 
‘m m at 7 ee 
_ ke (—1)’ a e” 
=> >) “ey 1)! (d4,)"* ( rot and )oo < y < r), 
v=1 p=] (12) 
where the cfs quantities are coefficients in resolution 
into partial fractions: 
> S Ym m e. 
“2 at = ay (13) 


Also of interest are cases of F (t) functions with the 
period T as well as functions 
F(t) = — F, (t + T). (14) 


Under the above assumptions, in case (14), the peri- 
odic solution of system (1) will be written thus; 


a fT 
(= Y\FEt— Yar, as 
k=10 
where 
Piss (y) = —F 2 —_* é (16) 


Té,q is the Kronecker symbol, 





ee es ceeceed 


If Eq. (3) has simple roots, we also have 





Yis(y) = 
n A bev 
=2 Bray cu for O0<y<+, 
e * +4 
YE) = —¥h(y+)- an 


The case wherein Eq. (3) has multiple roots is con- 
sidered in the same manner as above. 
The derivation of equations contained in this para- 


graph is given in Appendix I. 


3. Self-Oscillations of a Steady-State 
System Having a Single Nonlinearity 








Consider the system of equations 


Ze = D) beatae + hef (3) 


(4 =1,2,...,n), 


6o= Py J s%e, (28) 


where b orn she tags 9 1), 90 are con- 
stants, and 


f(0) is a sectionally linear function. 
We shall introduce the notation 


Hy (0) = >} heAce (\), (19) 
s=1 


‘M (d) = + D jee) = — Dd) Dd hainAnr (0). 


k=1 s=1 k=—1 
(20) 


We shall asume that the polynomials (20) and (3) do 
not have common roots, Then without restricting the 
generality of treatment, we can consider that, for any 
number T, Eq. (3) does not have roots of the form (4), 
which can always be secured by introducing a new non- 
linear function: 


f(s) = f (6) —¢s, 


where c is a certain given number. 
We shall assume that system (18) has a periodic so- 
lution of a certain period T: 


a, (t) = 22 (t+ T) (s = 1, 2,..., a), 


a" (t) = 0" (¢+ 7). (21) 


f(t) =7(0" (t)) (22) 





will also be a periodic function of the T period. 
Consider the linear system 


Tp = z Dyula + Myf (t). (23) 


k=1 


According to the above assumptions, we have a non- 
resonance case, In this case, the unique periodic solution 
of system (23) will obviously coincide with the solution 
(21) of system (18). 

Since in the given case 


Fy (t) = yf (t), (24) 
then, according to (5), 





T 
ar()=\OFt—/(e)dr, (25) 
0 
where 
Try — — 1 5 He(¥m) inv 


(26) 


and, moreover, 


T 
of (t) = \ @? (t — x) f (x) dt, (27) 


where 
@7 (y) = 





n P ae Yon) m 
DAM = az DB » Bey’ eat 


(28) 


Let us now consider the function (22). We shall as- 
sume that the periodic solution type (i.e., the sequence of 
passage over different portions of the nonlinear character- 
istic) is given. Then, we have 


f(t) = 007 (t) +0, for tis<t<h 


(i= 4, 2,..., m), (29) 
where a, and b are assigned numbers, some of which can 
be equal to each other; m is the over-all number of tran- 


sitions from one portion ‘of the nonlinear characteristic 
to another during a single period; t,= 0; 


Oty <lgxc.ccctu al 


(30) 
are unknown times which must still be determined, 











Let us introduce the periodic step functions 
A(t)=A(t+T), A(t)=q@ for H<t<ch, 


B(t)=h for t-y<t<k. 
(31) 


B(t)=B(t+ T), 
Then, (29) can be rewritten thus; 


{(t)= A(oT ) + BY. (32) 


By substituting this expression in (27), we obtain 
the integral equation 


of (t) = 
T T 
*\ ©? (t — t) A(t) o? (t) de + \ ©? (t — 1) B(x) dr, 
whose kernel 
R, (t, t) = 07 (t—1) A(t) 6s 
and the free term 
 y 
M (t)=\ 07 (t — 1) B(x) dt (35) 
0 


are known functions of t,. 

We shall limit our consideration to values of t, for 
which the denominator of the resolvent of kernel (34) 
does not become equal to zero, Then, in correspondence 
with the general theory [8], we can write 

T 


oT (t, th, t2,..., T)=M(t)+ \ R(t, t) M (x)dr, 
. (36) 
where R(t,T) is the resolvent of kernel (34), 
If 


Rit, Di< a, 


the resolvent can be obtained by the method of succes- 
sive approximations. 

In the general case [8], the resolvent is given by the 
equation 








= d. (t, 
> (— 1)" + = 
Rit, tT) =7= x : 


while 


d,(t, D=Ri(t, 1), d =4, (38) 


and all remaining qd. (t, T) and qd, can be successively 
found from the recurrent relations 
T 


duets = \ de(e, t)d, (39) 


v0 


T 
dy (t, t) = Ry (t, t)da — n\ R, (t, y)dn—s(y, t) dy. 
’ (40) 


Considering (31), (39) and (40) can be rewritten in 
the following form: 


dy, (t, t) = dnai®" (t — z)- 


tj 


nda) OO — dvi, dy, hact<ty 


i=l tj-y 


dna = d,@" (0) 5} a; (t; — ti-4)- (41) 
i=1 


tj 


aa 
-n p's a;\ dt { 7 (t-—y)d,_,(y, t) dy. 
cag 0 ae 


It can be shown that 
d,, (t, t) = A(t) ba(t, T), (42) 


where b, (t, T) are continuous functions for t rT, 
By using (42), Eq. (36) can be simplified and reduc- 
ed to the form. 


of (t, A oe 
-- es 
—\L6, DB@ae= > bi | Le, dat, 49) 


i==[ ij nal 


where RU, 0) 


A(t) 





Lit, t\) = 
(44) 
is a regular function [9]. 


In practical calculations, one must always take the 
finite number of terms in the numerator and the denom- 
inator of (37). We can derive equations which make it 


possible to estimate the error resulting from this proce- 
dure. 


It is known [8] that 
Idn| <n? (MTY", 
; nti 
|dn(t, t)|\< ae (a+ 1) * (MTy***, 
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where M is the maximum value of the modulus of func- 
tion (34), which is found from (31), (7), and (20) in each 
concrete case, 

Let ng be an integer which satisfies the condition 


n, > (MTe)*. (46) 


Then by using (45), we can derive the error esti- 
mates; 























co 
d, " (MTe)™ 4 
2 at|S 2 atl Sie Vana 
— | 4 St Des (MTe)™*1 ‘ 
ish V (mq + 1)™+! Yng+1—MTe 








(47) 


Let us now assume that the expression for of (t, ty 
-++» T) has been found, Then, in order to determine the 
unknown times t,, it is necessary to use conditions which 
determine the transition from one portion of the linear 
characteristic to another (“switching conditions") [10). 

For instance, consider the simplest case wherein the 
transition takes place when the o coordinate attains a 
certain given value. 

Then, the following condition must be satisfied: 


of (ti) = ass? (48) 


where % ass are Certain assigned numbers, 
By using (43), we obtain 


OF (ti, ty)... , tm—a, T) = Sigee (49) 


Expressions (49) represent the equations of periods 
which were to be determined, For more complex switch- 
ing conditions [10], the equations of periods are composed 
in a similar manner. 

If a computer is used for the solution of system (49), 
the considerations in [4] remain valid. In this, it is ne- 
cessary to take into account that the periodic solutions 
of system (18) do not actually correspond to all systems 
RRP Big Pe pyemenla It is nec- 
t, for i< j, and that o attains its value 
oj a eo for t= t4, Moreover, for physical rea- 
sons, there may be some requirements imposed on the 
derivative of o, etc, [10]. 


Let us briefly consider the case of a symmetric pe- 
riodic regime. 

Let us assume that the number of transitions is 
m = 2n, and that 


ty — Cys = than —tein-, kK<n. (50) 





The nonlinear characteristic must be such that 


Onte=4,, Onix=—dy, ken. (51) 


We shall seek the odd periodic regime character- 
ized by the conditions 
T T T T 
Ls ()=—ai (1 +3) » 6 (t= —o7(t +). 


If (50) and (51) are satisfied, (52) 


A(t)= A(t+->), B(t)=—B(t++), 


(53) 
and, consequently, by taking (52) into account, 


f()=—7(t +). 


By using Eq. (15) and following the above reason- 
ing, we arrive at the integral equation 


of (t) = \ YT (t —1t) A(t) a" (t) de + 


+ | Fe—v) B(x) dr. 


Ce—rw!y Cee 


(54) 


The solution of this equation is written analogously 
to (43): 


T 


oT (t, thy ++ ybm—ay —z-) = | Lilt, 2) B(e) de. 


(55) 


nh 


Thus, in correspondence with (50), the number of 
undetermined times which have to be found is reduced 
by half. The equations of periods are composed as in the 
general case, 

The derivation of equations in this paragraph is giv- 
en in Appendix IL. 


4. Forced Oscillations of a Steady-State 
System Having a Single Nonlinearity 








Consider a system of equations differing from (18) 
by the presence of a certain period T on the right-hand 
sides of the periodic time functions; 


= 3) Prete +h i), 
Ty 2 kata - uf (3) + felt) (ek = 1, 2.60% n). 


6= Dieter Te (t) cage Te (t +T), 
(56) 











We shall seek periodic solutions of the system (56) 
having the period T. . The solution type is considered to 
be given beforehand, In contrast to the self-oscillation 
case, the initial instant of time corresponding to the pe- 
riodic regime will be unknown here, We shall assume 
that the periodic regime (21) is determined. Then, by 
consideration of, instead of (56), the linear system (1), 
wherein 


Fy (t) = hyf (67 (t)) + fe), (57) 
and by repeating the reasoning in the previous paragraph, 
we arrive at the integral equation 
T 
oT (t) = \ OT (t — t) A(x) o? (x) de + 


T 
+\O7¢—)B) d+ Mi), (58) 
where 
ee. & T 
M,(t)= Y D ie \ On —Dh(e)de 
s=lLk=1 0 (59) 


is a known function. 
By applying the superposition principle and the re- 
sults of the preceding paragraph, we find 


T 


tana) = \ L(t, t)B(x)dv + My, (t) + 


OT (f, bey. «0; 


T 
+ \ L(t, t) A(t) My (x)dt. (60) 
i) 

By use of the switching conditions we find the equa- 
tions of periods with respect to m unknowns to, ty, ..., 
tm-1° 
It is known that the action of disturbances within 
the period T on a nonlinear system under certain condi- 
tions can cause the appearance of periodic motions with- 
in the period pT, where p is an integer (“subharmonic 
oscillations”), 

The described method can readily be applied to the 
case of subharmonic oscillations. Actually, the f,(t), 
A(t), and B(t) functions can be considered to have the 
pT period. Therefore, the following integral equation is 


valid for the 0°’ subharmonic oscillations: 


pT 
oPT (t) = \ @PT (t — t) A(x) o?T (x) de + 


pT 


+ Orr(¢— a) Bar + Mi, 


where 


The solution of Eq. (61) is written similarly to that 
of (60). By using this solution and the switching condi- 
tions, we find equations of periods for the subharmonic 
regime. 


evmv 


(62) 


5. Approximate Determination of the 
Periodic Regime in the Form of a 
Trigonometric Series 











Let us first consider the case of self-oscillations. In 
all practical calculations, we can use a finite number of 
terms in series (28). In this, the integral equation (33) is 
transformed into an equation having a degenerate kernel 
[8] capable of exact solution, 


Actually, let 


T 
T(t) = \or (t — 1) A(x)o" (x) dr + 


= 
+ \Or¢e — 1) B (x) dr, 


0 
(63) 
where 


ae, ext) 





©7(t — 1) ae ae 4D (64) 


From (64) we can conclude that the solution of Eq. 
(63) is of the following form: 


n 


a (t)= Dy se” 


i=—n 


(65) 


where s, are constants to be determined. By substituting 
(65) in (63) and equating the coefficients of equal har- 
monics, we obtain 


(i) M (%) (66) 
s= at Ds Aj_p + Dewy 
where 
a, = _ b a,(e~** fs ek ‘e—1), 
s=1 
=—— >) b, (6% '#— oA “o—1 
-_ 2 . : } (67) 











are the coefficients of the Fourier function (31). 
After solving the linear system (66), we find 


1 7 M (p,) 
t)=z 2 DOr) ber, (68) 
k= 





Si (t;, te, “#9 


where A is the determinant of system (66): 


M M(x) 


® =| ae — D(%,) (,) 


(69) 





and A ; are the corresponding corecton. 
By substituting (68) in (65), we find O° (t, ty, te, 
, T) which, together with the switching conditions, 
yields the system of periodic equations which was to be 
found. In particular, for the case (48), we have 


n 


2! Si (tas tay -- 


k=—n 


. T)e*t= 


= Siass 


(i= 1, 2,..., m). 
(70) 


In the case of the odd symmetric regime, the num- 
ber of time unknowns will correspondingly decrease, 
while equations (66) will contain only quantities s; with 
odd subscripts. 

In the case of forced oscillations, the approximate 
integral equation obtained from (58) will be of the fol- 
lowing form: 

T 


3 (t) = ( G @—1) AWS (nae + 


0 


so 


+) O* (t — 1) B(x) dt + M, (0). ig 


0 


Its solution will be 


ST (t) = M,(t)-+ >) se%', (72) 


i=—n 


where the values of §; are determined by the relations 


si 1 <= ;M(%,) 
% = — m+ s 2 (Dn) by + mx) Bri (73) 


and m, are Fourier coefficients in the expansion 


M,(t)= > me*'. (74) 


i=——oo 


By using (73), (72), and the switching conditions, 
we obtain the equations of periods. 
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6. The General Case of a 
Sectionally Linear System 








Let there be a system which, for t 


_1< t< tp is de- 
scribed by the equations 
n 


T= D) Bate + fit) 


a=] 


(75) 


where Ba and fi, are certain given time functions, 
We shall introduce the periodic functions 


Bya(t) = Bya(t+T), Bra (t) = Bial(t) 
for Gtiis»<t<h, 
ox (t)=%¢+T7), oQ=AQ (% 
for ho, <tc. 


Then, if the periodic regime (21) has been attained 
in system (75), we have 


= >) Bye (t)ar +o, (t) (e=1, 2... 


a=] 


+» n).(77) 


After we separate the constant terms from the Four- 
ier expansion of the B,. of) function, (77) will be written 
thus; 


f= Doel + 3D Bult) + ol 
Tt es (78) 
(4 =1, 2,..., m), 
where 
danse . P 

b.e=a YD (2 A(t) dt, Bya (t) = Bya (t) —Bya- 

i=, i-y 
(79) 

Assuming that 


Fy, (t) = >) Byara(t) + x(t) (k= 1, 2... 0) 
a=1 (80) 


and by using (5), we arrive at the Fredholm system of 
integral equations of the second kind; 


a (t)= >) \ Rel, tax (t)dt + P,(), (81) 


k=10 











where 


Ris (t, t) = s O7.(t—t) Bme(t), (82) 
is od 5 {oF (¢ — t) @, (t) dt. 
10 (83) 


Here, in contrast to cases previously considered, the 
®,, functions depend on the unknown times t; which, 
however, does not involve any difficulties in principle. 

As is known [8 and 9], the system of equations (81) 
can be directly reduced to a single equation with argu- 
ments t and T, which are defined within the square given 
by 0StSnT and 0=7T=nT. Therefore, the method for 
solving system (77) does not differ from those described 
above. In particular, as in Section 5, system (77) can 
always be reduced to the degenerate equation. 

After finding the solution 

a, = 22 (t, ty,..., T) 
by use of the switching conditions, which can be most 
diverse in dependence on the design and designation of 
the controlled system, we arrive at the equations of peri- 
ods. In particular, if we consider the system 


Tp = >) deata + Mya (t) 6 + b(t) 


a=l 


for tiey < t < ti, (84) 


instead of (56), we find that all the results obtained in 


Sections 3,4, and 5 can be almost entirely applied to 
this system with the exception of the fact that instead of 


function (31), the following functions must be considered: 


A(t)=A(t+T7), A(t) =ai(t) 
for tiy<t<ch, 
B(t)=Bt+T), BO=h 
for uu, t< &. 


(85) 


APPENDIX I 


Derivation of Equations in Section 2 
Let the Fourier expansion of the F(t) function be 
of the following form: 





F.()= > fpr er’ 


r=—oo 


(86) 


We shall seek the periodic solution of system (1) in 
the form of a Fourier series: 


co 
z! (t)= > Xo, Ove 


k=—co 


(87) 


By substituting (86) and (87) in (1) and by equating 
the coefficients of equal harmonics, we obtain 


sini Yee Oo fam Sis (Ym) 


(88) 
k=1 Dm) 
Considering that 
S| 
Sym = ge \ F, (0) em" dt, 
: (89) 


by substitution of (89)and (88) in (87) and by changing 
the order of summation and integration, we obtain (5). 

Let us further consider the Fourier expansion of the 
function of the period T, which coincides with the isfy) 
function for 0< y < T. We have 











Tn A (A,) Apt 
1 ks \"e e —dnt 
r\ 3G) WT _,e "at 
lt (90) 
_ 1 wy Ane Ge) 1 1 Ans (p) | 
r 2 Dh.) e—% TT PY) 


which proves Eq. (7). 
Equation (12) can be proved in a similar manner if 
we take into account that 
A 
bs 


d*— 1 
any—* 





{ = 


—i,)” (v—1)! : Bai Jee. (91) 


From (7), we further have 


s Ans (A,) 


®,, (T)—®,, © = Bay 
e 


(92) 
e=1 


whence, by taking into account the well-known equa- 
tions from [11], we obtain (9). 

Equation (8) can be obtained by using the following 
well-known property of Fourier series: 


®,, (T) + ®,, O 


Of, (0) = OF, (7) = 5 - (93) 





Equations (15)-(17) are proved in a similar manner 
on the basis of the fact that, in the case (14), the Fourier 
expansions of the F;(t) functions will not contain even 
harmonics. 


APPENDIX II 


Proof of Equations in Section 3 
The recurrent kernels for Eq. (33) will be calculated 
by using the equation 
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T 
R,(t, D= VR (t, y) R,_,(y. T) dy- (94) 
0 


By induction, it can be shown that 
R, (t, tT) = A(t) L, tt, T), (95) 
where 


L, = ©7 (¢—1),..., 


n Fi 
L,= D% | oF t—v Lu Day 
t= yy 


(96) 


are continuous functions for n> 1. 

Considering that d,(t, rT ) is a linear combination of 
Ry, « « » » Rn» we obtain (42), whence, by taking into ac- 
count the absolute convergence of the numerator in (37), 
we obtain (44). Further, by multiplying Eq.(33) by A(t), 
we find that the function 


s(t) = 67 (t) A(t) (97) 


satisfies the equation 


T 
8 w= A(t) @7 (t — 1) s(t) dt + 
0 


(98) 
A(t) ©? (t — t) B(t) dt, 


+ 


owt 


whose free term is a function which can be represented 
by the kerne) 


Ni (t, t) = A(t) ©7 (t — 1). 
(99) 


The solution of Eq. (98) can be represented in the 
following form: 


T 
e()={ N(t, ) B(x) at, (100) 
0 


where N(t,T ) is the resolvent of kernel (99). 
By direct calculation, it can be shown that the re- 
current kernels Ny, are equal to 


N, (t, t) = A(t) L, (t, 1) (101) 
and therefore 
N (t, t) = A(t)L(t, 1). 
(102) 
By substituting (102) and (97) in (100), we obtain@3 


It now remains to prove the estimates (47). From the 
first of the equations (45), we have 
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co co k 3 M 
‘m=1 2) lS Da = ((1008) 
k=N% k= Te 


Since, according to Stirling's formula [12)}, 


Tein (104) 
then, \*) 





< MTe\k &y / MTe \k 
ms (Ve) <2 (Va) 08 


k=N, 


If we assume that ng satisfies condition (46), then, 
by summing the geometric progression on the right-hand 
side of (105), we obtain the first of the equations (47). 
The residue of the numerator in (37) is estimated in a 
similar manner. 
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It is shown that the describing-function method, when applied totwochannel systems with 
modulation and a nonlinear element in the ac circuit, becomes an exact method for determin- 
ing self-oscillations for an unbounded increase in the carrier frequency. An approximate 
criterion for theself-oscillation stability in such systems with a relay component is formula- 
ted, The transfer function of a system which is linearized by self-oscillations with respect 


to slow control inputs is determined, 


In ac follow-up systems, the infralow-frequency er- 
ror signal is transformed into a signal with the carrier 
frequency 2, after which it is amplified and then again 
transformed into an infralow-frequency signal. Alternat- 
ing-current follow-up systems where the infralow-fre- 
quency portion of the circuit has two channels, in partic- 
ular, two identical channels and antisymmetrical cross 
couplings [1], are used to a considerable extent. Some 
automatic ac bridges [2], some follow-up hydraulic plat- 
forms, etc,, belong to these systems. If the ac factors 
and the two-channel infralow-frequency part of the cir- 
cuit are linear, such systems can be conveniently inves- 
tigated by the method of complex transfer functions [1) 
with the application of ordinary frequency criteria for 
stability and performance, However, in a number of 
cases, two-channel relay systems with a relay element 
in the ac signal circuit can be recommended, A similar 
situation occurs in the case of considerable and unavoid- 
able changing of the ac channel amplification factor 


basic operating conditions are the conditions of “spatial” 
self-oscillations of rather high frequency and small am- 
plitude. With respect to control inputs, which change 
slowly, the system behaves as a linear system due to vi- 
brational linearization [5 and 6] by self-oscillations, 

The purpose of the present article is the presenta- 
tion of a relatively general method for the investigation 
of similar two-channel systems with a relay element in 
the ac circuit,as well as the analysis of certain common 
properties of these systems. 

The method under consideration is suitable for the 
solution of actual technical problems. 


1. Structural Diagram and the System 





Transfer Functions 





The quantities x;jn and xgip are supplied to the in- 
put of a two-channel follow-up system (Fig. 1). The out- 
put quantities x; 5,, and xgout are subtracted from these 
quantities, and the following error signals are formed: 




























































































and in a number of other cases, In such systems, the & =—Zin—*out & = Zin Trout 
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Fig. 1. General diagram of a two-channel follow-up system with a relay element in 


the ac circuit, 
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The error signals are fed to modulator 1, where they 
are transformed into a signal with the carrier frequency 
Qs 


u = 8, sin Qt + e, cos Nt = e sin (Nt + Q). 


This signal is transmitted through the first linear 
portion of the ac circuit, which has the transfer function® 
W,(D), and then through the relay element and the 
second linear portion, whose transfer function is W3(D). 

The output signal z of the single-channel portion of 
the ac circuit i fed to the demodulator 3, where it is 
multiplied by sin Qt and cost. The z sin Mt and 
z cos Mt quantities enter the two-channe] infralow- 
frequency part of circuit 4. 

We introduce the complex quantities 


Zin = Zin + 2 in 2 out = 71 our + 1% out 
z = zsin Ot + jzcos Qt = jze™, () 
and we write 


@ = zin— Zour Z out= WwW (D) Z. (2) 


Taking into account the possible presence of anti- 
symmetrical cross couplings in the two-channel portion 
of the circuit, we shall consider W (D) as a complex 
function in the general case. 

The equations of the instantaneous values of vari- 
ables in the circuit single-channel portion are of the 
following form: 


u=esin(Qt+q), v= W,(D)u, 


(3) 
y= (2), z= W,(D)y, 
where f(v) is the nonlinear factor characteristic, 

The single-channel portion of the circuit consists of 
an ac circuit, and the variables in this circuit can be 
conventionally represented by complex amplitudes. 

We shall assume that the variable x = X, sin Qt + 
+ X, cos Qt has the complex “amplitude” X, = Xy + iXe, 
which is a function of time in the general case. 

In the Appendix it is shown that the complete sys- 
tem of equations in terms of complex variables is of the 
following form: 


& = tin — W (D)z, 
V = W, (D + jQhe, 
Y, sin Q¢+Y, cos Qt = Alf (V; sin Qt + Vz cos Q2)], 





Here, the A[f] symbol denotes the first harmonic of 
the almost periodic function f(V, sin Qt+V_cos Nt), 
Thus, ¥ = Y, + JY, represents the complex “amplitude” 
of the carrier frequency first harmonic at the nonlinear 
element output, to the input of which a carrier frequen- 
cy signal with the complex "amplitude" V = Vy + JV, ts 
supplied. 

The frequency spectrum of the z§ component begins 
in the vicinity of the carrier frequency Q and then ex- 
tends toward higher frequencies. 

Equations (4) correspond to the structural diagram 
with complex transfer functions and complex variables 
given in Fig. 2. The following fact should be noted here, 
The transfer functions of the single-channel linear por- 


tions of the circuit for the envelope in the considered 
system with a two-channel infralow-frequency portion 


are obtained by the simple substitution of D+ jQ for the 
instantaneous values of the D operator in the transfer 
functions, where 2 is the carrier frequency (the factor 4 
in front of W, is generated by the demodulator—synchro- 
nous detector). 

Such a simple rule for the formation of the transfer 
function for the envelope is valid only for sys- 
tems with identical channels in the infralow- frequency 
two-channel portion of the circuit. Only in such systems 
can the sinusoidal as well as the cosinusoidal output sig- 
nal components in the ac channel be used to the same 
extent. As is known [3], in ordinary ac follow-up sys- 
tems with a single-channel infralow- frequency portion, 
the formation of the transfer function for the envelope 
is somewhat more complicated, 


2. Steady-State Self-Oscillations 

As was mentioned above, self-oscillations represent 
the basic operating conditions of relay systems of the 
class under consideration. 

In order to determine steady-state self- oscillations, 
we shall assume that x; = 0, and we shall seek the so- 
lution of system (4) in the following form: 





Pitas shy tiie date 
6 = e,ea' , V=V ie" , 


(5) 


= — a4 ae 
= Ze * > Y= Y,e'%a' , 


where €, Vo, Zp, and Yq = const, and w, is the self- 
oscillation frequency. 















































z= : W,(D + jQ)¥ + 25. m_ Ss D=> everywhere, 
4 
fr; w,(D+j2) YeWw,(D +2) (0) }--+15 
out 
Fig. 2, 
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Due to the presence of the 35 term, this is an ap- 
proximate solution, and its form entirely corresponds to 
the describing -function method [4]. 

However, in application to the considered systems 
with a relay element in the ac circuit, the describing- 
function method has a special feature. This feature, 
which renders the describing-function method an exact 
method when applied to such systems, consists in the 
following. Here, instead of neglecting harmonics with 
frequencies which are multiples of the self-oscillation 
frequency, we neglect harmonics with frequencies which 
are multiples of the carrier frequency Q or of the com- 
bined frequency Q + w.. 





Actually, for Xin = 0, by substituting expressions (5) 
in equations (4), we find 


_& = — W (jaa) Zo, 
Vo = (j@a + jQ) &>, 
Yo= W Vo, (6) 


- 1 , \ & —* —jW, 
29 = > Ws (j@a + jQ) V9 + Zé ane 


Here, Wp = + is the amplitude response of the 


nonlinear element, which is determined by the relation 


Y, sin Qt + Y, cos Qt = Yy cos (@at + @,) sin Qé 4- 
+ Yosin (@at + g,) cos Qt = Yo sin (@gt + Qt + @,) = 
= Aff (V, sin Qt + Vz cos Q1)] = A {/ [Vy sin (@et + Qt + @))}, 


Y, = Y,c, 7, =V,e%. 


For a relay element without an insensitivity zone 
and hysteresis, which will mainly be considered below, 


Ww _ 4a 1 


as rf Ve F (7) 


where 
a = const. 


The %%e"J“a' component, the neglecting of 
which renders the describing-function method an ap- 
proximate method, contains higher harmonics, From 
Eqs. (26) and (27) in the Appendix, it follows that, in the 
given case, the lowest frequency of these harmonics is 
equal either to Q + wg (the nonlinear element response 
is nonsymmetrical, and there is a constant component at 
the ac channel output) or 2(Q2 + w), while, in ordinary 
cases where the describing-function method is used, the 
lowest frequency is Qu. 

In the majority of practical cases, the carrier fre- 
quency 2 is considerably higher than the possible self- 
oscillation frequency. Thus, for angle-tracking systems 
with conical scanning, Q is equal to (607 - 2007) sec” 
For ac systems with standard frequencies, 2 is equal to 
(1007 - 8007) sec™*, Therefore, in the majority of 
cases, the Ze~J“at component is almost completely 
suppressed in the two-channel portion of the circuit, 
which always acts as a filter for low frequencies, and 
the neglect of this component yields only slight inac- 
curacies, 

Neglecting this component, we write 


&y = —W (joa) %, Vy = Wy lj (Q+ oa) &, 
Y, a WY, Zo om Wali (Q + @a) Y,. ws 
Hence, 


V.= — W (joa) Wi [j(Q + 





+ @)] W,[j(Q + @2)] WV, 


or 
1+ W (joa) Wy = 0. (9) 


Here, the quantity 


WwW ' (J/@q) = 
4W (joa) Ws 1) (2+ oe) Wali (2 | w2)) 


can be called the over-all transfer locus of the system 
linear portion, and the quantity 


W, (D) = +-W (D) W,(D + jQ)W,(D 4-2) 


can be called the complex transfer function of the iinear 
portion. 
Equation (9) and the equations equivalent to it, 


W=—Wa', Wi'=—M, 
where Wj]! and W,! are the inverse loci of the linear 
and the nonlinear portions, respectively, whose shapes 
coincide with the relations of the describing function 
method in L. S. Gol’dfarb's form [4]. The difference 
consists in the expression for W; and the transformation 
of the relation into an exact equation for Q-> ©. 

It should be emphasized that, in the considered two- 
channel systems, where one portion of the circuit oper- 
ates with alternating current, the transfer locus, 


W) (jo) = —- W (jo)W, (jo + jQ)X 
x W, (jo + jQ) 
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is generally not symmetric with respect to the real axis 
(i.e., Wy(jw) # Wy(- jw)) even if there are no antisym- 
metrical cross couplings between identical elements in 
the circuit two-channel portion. Therefore, the transfer 
locus must be considered in the range where w, changes 
from -® to ®, As usual, the values of wa at the points 
of intersection between W,(jwg) and W;,! (or between 

Wj *(jw,) and Wy) correspond to the frequencies of pos- 
sible self- oscillations, and the Wn = Wn(V,) values at 
these points determine the amplitudes V, of possible self- 
oscillations (Fig. 3). 











Fig. 3, 


3. Stability of Self-Oscillations 
In order to find the conditions for the stability of 
self- oscillations we shall represent them in the form 





& = (e, + Ae) ee, V = (V+ AV) ec", on 

Y =(¥,+ AY)e**', 2=(z,+ Az)z0', 
where, as above, the quantities with the "0" subscripts 
correspond to the amplitudes of steady-state oscillations, 
and the second terms on the right-hand side represent 
deviations. Thereafter, we shall consider the deviations 
to be small, i.e., we shall consider stability in Lyapunovs 
sense, 
In order to write the equations in terms of devia- 
tions, it is necessary to establish a relation between AY 
and AV at the relay element output and input. 

For the relay element under consideration, 


V. , 4 V 
hezep hagey 


V =) Vi4- V3. 


Hence, with an accuracy to small higher orders, we 
obtain 


914 





4 AV, V , 
AY, = = 8 kee = (V,AV, + V.AV,) , 
0 


" 4 AV. V , ‘ ; 
ate= a (Te WidMs + VAV2) |: 
However, ; 

VY, => V, cos (al ws 9), V, = - sin (Wal ‘ p). 


Therefore, 


AY = AY, - jAY. oe 
ha AV eial +) { 


Ts” 


|- Sin (@at -+ @) avai} = 7 | av — 





COS (Wat -+ @) AV, 


2i(wgt-+) 2i(wgt++) 
1 , -—! 

e—1t ae, jp 2 Av, | =. 
(11) 


= [AV — av ent) 


Here, AV® is the complex conjugate of AV. 
By substituting expressions (10) in (4) and consider- 
ing (11), for Xj, = 0, we find 


Ae = —W (D + jaz) Az, 
AV = W, (D + j@a + jQ) Ae, 
Be = + W,(D + joo + jQ) BY + ze", (12) 


AY = = [AV P= AV ea! +) . 


Thus, in the general case, the equations for the self- 
oscillation deviations contain a periodically changing 
coefficient e*/(Wat+), If the self-oscillations vary 
slowly, and this is the most typical case, AV and AV* 
are functions which vary slowly in comparison with 
eti(Hat+?) and the effect of the second term on the 
right-hand side of the expression for AY can be neg- 
lected. The 25 e}“# term, which, according to what 
has been said above, contains higher frequencies, can be 
even more readily neglected, 

In connection with this, the approximate equation 
of deviations of self-oscillations, which is obtained from 
Eq. (12), is of the following form: 


E = oa Wy (D -b jw.) | AV = 0, (13) 
where 
W (Dd -|- J@a) = a W (D -- j@a) 


Wy (D <- ]@a + jQ) Ws. (D -!- ]@q +- fQ) 


is the over-all transfer function of the circuit linear 
portion for the envelope of the self-oscillation frequen- 


cy Wa. 











If we know the equation of deviations, we can use 
any criterion applicable to equations with complex co- 
efficients for determining the stability conditions, By 
applying the Nyquist criterion, we find: if the transfer 
function W;(D) of the linear portion of an open-loop sys- 
tem does not have poles in the right half-plane, those 
and only those self-oscillations are stable for which the 
points (1/2a)Vo, j0 (Vg is the self-oscillation amplitude) 
are not encompassed by the transfer locus W;(jw) when 
w changes from -@ to , 

By comparing this condition with Eq. (9) for steady- 
state self-oscillations and the corresponding graphic plot, 
we readily arrive at a different formulation of the self- 
oscillation stability criterion for the relay system under 
consideration: if (-1, .j0) is the intersection point between 
W) and W_* that corresponds to the given self-oscilla- 
tions, then, for a W,(D) which does not have poles in the 
right half-plane and if w varies from -® to ®, these self- 
oscillations will be stable only if the transfer locus 
W;(jw) does not encompass the(—21, j0) point, whose 
distance from the imaginary axis is twice as large in 
comparison with the first point. Thus, the self-oscilla- 
tions with the wa, and wa, frequencies in Fig. 4 are un- 
stable, and self-oscillations with the wa, and wa, fre- 
quencies are stable. Hence, it follows that self-oscilla- 
tions with the largest amplitude are always stable. 

It must be emphasized that the above stability cri- 
terion for self- oscillations corresponds to the approxi- 
mate equation (13) and is valid insofar as this equation 
is valid, 

Another fact should be noted, As can be seen from 
the above, stable self-oscillations with a positive as well 
as negative frequency w, can exist in the systems under 
consideration. The sign of wg in two-channel systems 
corresponds to the sense of rotation of the tracing point 
or to the direction of the actual three-dimensional coni- 
cal motion in two-channel follow-up systems with me- 
chanical output quantities. The frequency w, can have 











Fig. 4. 


the same sign as the carrier frequency 2,or it can be of 
opposite sign. In angle-tracking systems with conical 
scanning, the first case corresponds to the adding of the 
scanning and self-oscillation angular velocities, and the 
second case corresponds to the subtraction of these ve- 
locities, In any system of the class under consideration, 
the first case corresponds to the operation of the circuit 
single-channel portion at the combined frequency, and 
the second case corresponds to operation at the frequen- 
cy difference. 

The realization of one or the other direction of self- 
oscillations depends on the initial phase shift in the cir- 
cuit (the phase shift for w, = 0) and the shape of W)(jw). 

The change in the initial phase shift, produced, for 
instance, by changing the carrier frequency phase in the 
modulator, means that the entire transfer locus W)(jw) is 
rotated in the corresponding direction. In such rotation, 
the self-oscillations in one of the two possible directions 
can always be rendered unstable (see Fig. 4). 


4, Reproduction of Control Inputs 

With regard to the reproduction of control inputs, 
which change slowly in comparison with self-oscilla- 
tions, the behavior of the two-channel systems under 
consideration is analogous to other vibrationally linear- 
ized nonlinear systems [5, 6]. For small deviations, 
the dynamic properties of such self-oscillation systems 
are close to the characteristics of linear systems. 

In order to find the transfer function of a self-oscil- 
lation. system with respect to the slowly changing con- 
trol input Xjp, we shall write the variables in the follow- 
ing form: 





e = e,e*a' + Ae, V =V,e2' + AV, 
Y =Yyeo! + AY, 7 = Ze"! + Az, 
where the first term on the right-hand side in each ex- 
pression corresponds to self-oscillations, and the second 
term, according to conditions, changes slowly in com- 
parison with the first. 
By substituting these expressions in equations (4), we 
find 
aie! + Ke = — W (jug) ee! — W (D) Be, 
Veta! + AV = W, (jQ + ja) bea" +W,(D+ 
+ 0) Ke, (14) 
ite! + Az = SW, (j2 + joa) Yee"! + Wy (D+ 
+ jQ) AY + z,; Y,sin(Qt + @at + 9) + 
+ AY; sin Qt + AY, cos Qt = A {/ [V, sin{Qe + 
+ Gat + @) + AV, sin Qt + AV, cosQz]}}. (15) 
In the considered case of a relay element with the 
ideal characteristic f(x) = asign(x), the first harmonics 
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in the expansion of the f function are expressed by the 
well-known equations [7] 


7(V, sin (Qt + wat + —) 4- AV, sin Qt + 016) 
-- AV, cos Qt] = asign [V, sin(Qt + o t+ @) + 
+ AV sin (Qt + 9,)} =~ NW sin (ax + Wat + @) + 


+ 22 Msin (Qt +91) 4 gg = PE, Se ee 


where 
AV; 
Pi = arctg WW,’ 


Nab a()a1—t (+. 


moor Lat (8t)- (aot) « (42) + 


=14+4-(F)+... 


and E and K are complete elliptic integrals of the first 
and the second kind, respectively. 

Here, the quantities with the "0° subscript-the self- 
oscillation amplitudes~are always constant, and the in- 
crements AV, A€, etc, change slowly or are constant. 

From (15) and (16), it follows: 


4a »Vo 


(17) 
where 


AY= ** M Av. (18) 


By equating the coefficients in front of e/“at and 
the slowly changing terms in equations (14) and by 
eliminating the intermediate variables by means of ex- 





pressions (17) and (18), we find. 
1+ W, (jo.)W, =0, (19) 
2a , 
* wer Wi(D) | 
Tout = — Zin 
1+ Ww. W} (D) (20) 
Here, 
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W; = +-W (D)W,(D + jQ) MW, (D + jQ), 


and W;(jw,) has the previous value; 


Wi (jee) = + W (jaa) Wy (j2 + joe) Ws (jQ + jag). 


In the general case, for large = deviations, the 
0 
N and M quantities basically depend on these deviations, 


In this, Eqs. (19) and (20) constitute a system of in- 
terdependent nonlinear equations, and the self-oscillation 
amplitude depends on the control input. 

' However, the basic system operating conditions are 
those conditions for which the deviations, which are 
caused by control inputs, are considerably smaller than 
the self-oscillation amplitude, i.e., AV/V»_ <1. Under 
these conditions, it can be considered that N = 1 and 
M=1, 

AV 

Actually, even if Ve = 0.5, N = 0.95 and M = 1,02, 

Thus, in the os“ = 0.5 interval, the N=1 and M=1 


0 
equalities are satisfied very accurately. If these equali- 


ties are satisfied, Eq. (19) is transformed into the equa- 
tion for steady-state self-oscillations, 


1+ W) (joa) W, = 0, (21) 
which has been considered above, and Eq. (20) corre- 
sponds to the — function of the linear follow-up 

1 
® (D) = sae ae 
i+ 








, Zout =©® (D) Zin 
(22) 


we. (D) 





which, in the open-loop state, has the transfer function, 
2a 
nV, W, (D) = 





= F+W (D)W,(D + jQ)W,(D + jQ). 


It should be noted that the characteristic equation 
which corresponds to the transfer function (22) for slowly 
changing inputs coincides accurately to the jw, term 
for D with the characteristic equation (13) for the estab- 
lishment of self-oscillations, 

This means that, if the self- oscillations are stable, 
a follow-up system with self-oscillations is similar to a 
stable linear system with regard to the reproduction of 
slowly changing inputs, 

The self-oscillation relay system under considera- 
tion has the property that its reaction to slow disturb- 
ances is independent of the absolute magnitude of the 


open-loop gain. 








t 


—= ./45 ee 











Actually, it follows from Eq. (9) that, if the modulus 
of the open-loop gain changes slowly, the self-oscilla- 
tion amplitude V, increases in proportion to the gain: 


w;'= Vo = —W, (jas). 


2 
The = W,(D) quantity, which enters Eq. (22), re- 
0 


mains constant, which determines the invariability of 
the closed-loop system transfer function. 

Here, it should be added that, if only the gain of the 
first factor W, of the circuit single-channel portion 
changes, which occurs often in practice, then not only is 
the over-all transfer function constant, but also the am- 
plitude €9 of self-oscillations at the point of error meas- 
urement is constant (V, changes), 

In conclusion, it should be mentioned that this 
method and the considered facts can be relatively simply 
applied to the case of a relay element with hysteresis 
and an insensitivity zone. 


APPENDIX 

It is known [8], that the transfer function for com- 
plex amplitudes (for the carrier frequency 2 envelopes) 
is obtained from the transfer function for instantaneous 
values by substituting D+jQ for the D operator. 

Thus, 


V=W,(D+;2)0, U=e. (23) 
The quantity, 
y = f (v) = (V1 sin Qt + V2 cos Qt), 


where V = V; + jVg, can always be represented in the form 
y =f (Vi sin Qt + Vz cos Qi) = 
= Y°+ Y,; sin Qt + Y2 cos Qt +- 


CO 
+ >} (¥( sin ine + Y{% cos iM), 
i=2 
where Y°, Y;, Yo, Y;{i), and ¥,(i) are functions of time in 
the general case, 
In particular, if V,, V.=const and the relay element 
has a z-shaped characteristic (f(v)=asignv), then 


¥°=0, Y=/VY¥i+¥i = a, 


at = 2 (there is no phase shift in the relay element), and 
1 





YP =VP =o, YO = VPs OP = fa, 
. . (24) 


. 
* ©. @€ © 2.0°8900¢0 6 € © £6 1 620.226 ee 


Thus, the relay element output can be character- 
ized by the complex amplitude of the first harmonic 
Y= Y,+ j¥, and by the complex amplitudes of higher har- 
monics Y(i) = y,(i) + jy,(i), 





In correspondence with this, the complex amplitude 
of the first harmonic of the factor W, output will be ex- 
pressed by the equation, 


Z = W,(D + jQ)¥, (25) 


and the complex amplitudes of higher harmonics of this 
quantity will be determined by the relations, 


Z® — W2(D + ijay Y™. 


The instantaneous value Z of the ac channel output 
is equal to 


z= Z,sin Qt + Z, cos Qt + Z, (26) 


where Z,+jZ, = Z, and Zg is the component correspond- 
ing to high harmonics oak the component of the zero 
carrier frequency (if Y°#0, also such a component is 
present), 

However, according to (1) and (26), 


z =: sje 7M (Z, sin Qt + Z,cos Qt 4- zie jor __ 


a: + |Z) -- > (Z; ae }Zs) e 2/at jz,eI™ 


124: 


2 oe 





(27) 


Here, 


-6 


2 ~ ry hee oo 8 + (Zi ee |Z) Caos 


is the component which, under steady-state conditions, 
has components with a frequency of the order of 9 and 
higher, This follows from the fact that Zg contains har- 
monics with double and higher frequencies. Thus, the 
closed systems of equations for the circuit under consid- 
eration is of the following form: 


@ = zin— W (D)z, 
V = Wi(D + jQ)e, (28) 
Yi sin Qi + Y2cos Qt = f (V; sin Qt+ V_cos Qt) 4 > . 
pee oy TF — 
t= —-We@ + fQ)F + Za 


where = denotes higher harmonics, 
Example. Consider the case where 


W1(D) = const, Ws (D) = e*?, W = Dat . 


This case corresponds to the inertialess first ele- 
ment in the ac channel (for instance, a broadband ampli- 
fier), to an element with a constant lag at the relay-ele- 
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ment output (in particular, it can reflect the actuation 
time of an actual relay), and to a two-channel final ele- 
ment in the form of inertial and integration elements 
which are connected in series (for instance, electric mo- 
tors). 

Considering that a phase shift of the reference sig- 
nal can occur in the modulator (or in the demodulator), 
we assume that 


W (D) == ke. 


where k; is generally a complex quantity. 
Thus, in the example under consideration, the over- 
all transfer function of the linear portion is* 


Wy (D) + W (D) W1(D +- /Q) Wa(D + jQ) = 


4 Rye Ot) — gO 
2 DAFTD) DTFTD) ' 





where 


k = Bs kie 


—jra 


For a relay element, expression (7) holds, and Eq.(9) 
for the self-oscillation frequencies and amplitudes as- 
sumes the following form: 


> —jt@, 


ke ge nVo 
j@, (1 + /T@,) a 








The corresponding graphic plot is shown in Fig. 5. 
The root locus portions corresponding to Twag>™ are not 
shown, since the approximation of dynamic properties 
by means of a constant-lag element ceases to be valid 
in the region of frequencies higher than /r for arelay 
and many other elements, 











Fig. 5. Root locus of Wj(jw) for the 
example under consideration. 
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In dependence on the ¢=arg k values, it is obvious 
that one or two self-oscillation regimes can exist in the 
system under consideration, which differ from each other 
by the magnitude and sign of wa, (the direction of coni- 
cal motion) and their amplitude. 

The self-oscillation equations in terms of real quan- 
tities are of the following form: 


To, + arcigTo,= + > +9, 
* 4ak 
lel i+ Te 





Vo 


If r =0, self-oscillations with finite frequency can 
occur only for ¢ #0, i.e., if some cross coupling is pre- 
sent. Under these conditions, for 1/2 >¢> 0, only self- 
oscillations with negative frequency can exist: 


and for-( 2)< ¥ < 0, only self-oscillations with positive 
frequency can occur; 


o,= 4. are te(¢ ++) . 

As follows from the above, the self-oscillations are 
stable in both cases, For T=0 (only integrating ele- 
ments are present in the system two-channel portion): 


4 eetate 9. 


@® = —___—_—____.. 
. tT 1 | @, | 


In this case, self-oscillations with positive as well 
as negative frequency can exist either for “= 0 or #0, 

However, it follows from the above self-oscillation 
stability criterion that, for ¢> 1/6, only self-oscillations 
with negative frequency are stable, and, for ¢<—(1/6), 
only self-oscillations with positive frequency are stable, 

In reproducing slowly changing control inputs, the 
self-oscillation system under consideration has, in ac- 
cordance with (22), the following transfer function: 











2ak e*D 
® (D) = mo . z 
Vo 


In particular, for T=0, because of the relation 
We" aial , the transfer function assumes the following 
form: 


i —tD : 
a lw. Je +39 |, |e 


2De™? + |W, | e'* 





® (D) = 


D + lug |e“? +# 














Here, w, depends only on the phase shift ” and the 
lag T: 
£x/2+9 

- ; 


Wq = 





Thus, the transfer function for slowly changing con- 
trol inputs is invariant with respect to changes in the co- 
efficient k modulus. 

This illustrates the general aspect of the determined 
self- adjustment property of the considered self-oscilla- 
tion systems, 
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A method for the calculation and experimental determination of the dynamic frequency- 
response curves for hydraulic couplings is presented, A brief description of the equipment 
used in experiments and the results of the experimental determination of characteristics 


of some couplings are provided. 


Hydraulic couplings in the shape of pipes filled with 
a liquid are used in automatic control systems (ACS). In 
designing ACS, the dynamic properties of such couplings 
must be taken into account. A method for the calcula- 
tion and experimental determination of the dynamic 
frequency-response curves for hydraulic couplings is 
given below. 

I. A. Charnyi [1] performed a theoretical investiga- 
tion of the effect of piping on the accuracy of the pres- 
sure gauge which registers pressure variations, The au- 
thor based his investigation on the N. E. Zhukovskii equa- 
tions of the hydraulic shock theory [2], which had an ad- 
ditional term whereby the pressure-head loss due to 
friction was taken into account, 

Following I. A. Charnyi [1], our theoretical investi- 
gation of the dynamic frequency response of hydraulic 
couplings will also be based on these equations: 


op __ (av Op _ Ov 
—F = (5+ 20), “se OO 
K=—“Q3, a= 
od 
i+ Fe, (2) 


where p and v are the pressure and velocity at the pipe 
section which is at the distance x from the origin, re- 
spectively, p is the liquid density, Ky is the actual elas- 
ticity modulus of the liquid, E is Young's modulus of the 
first kind for the pipe material, 59 is the pipe wall thick- 
ness, d is the pipe inside diameter, and v is the liquid 
kinematic viscosity coefficient. 

The dynamic characteristics will be determined for 
two cases: for piping without discharge and for piping 
with liquid discharge, 

For the first case, the boundary conditions will be 
written in the following form: 


z= 0, p= peciat, (3) 


z=l, v=0. (4) 
For the second case, considering that the discharge 


velocity is approximately proportional to pressure, i.e., 
v= Mp, we obtain the following boundary conditions: 


t=0, p= pei, (5) 
=l, v=Mp, (6) 


where 1 is the piping length, and M is the proportional- 
ity factor. 

We shall consider only forced oscillations, and we 
shall seek partial solutions of Eq. (1) of the following 
form: 

p= (Acosnz + Bsin nz) e, (7) 
v = (Ccosnz + Dsinnz) e. (8) 


The values of A, B, C, D, and n are found from 
Eqs. (1), (7), and (8) and boundary conditions (3)-(6). 
Let us consider the first case. Here, 4 — Po, B= 


= ; tw 
Potgnl, C= —— Po tg nl, D=—-—y Pw 


n= (@' —ia'), 








’ 


where o’ = /? wt + 40%? + w?* 
are 2 


? os / ee 
= ; F 


AK 








c is the sound velocity in the piping. 


If a is small in comparison with w, a'=0, andw*=w, 
By substituting the values of A, B, and n in (7), for 
x=], we obtain: 


Pi = (po cosnl + py tg nisin nl) ett, (9) 











ve 








After transforming (9), we find 


tol 











7 Pot 
— rt) n’! ein Oe gh 2! 
cos 7." 2. deceit hadi W 
or 
Pi = PoA (@) etlot—(o)] (10) 
where 
1 
_ (o) 2 l l 
a ae (11) 
cst: aoe S Sane 
ch P= sin’ — 


is the amplitude response, and 


Tre @’l 
@ (@) = arc tg (th— tg ) (12) 


c 


is the phase response, 


It follows from (11) and (12) that the amplitude re- 
sponse A (w) has not only one, but a number of resonance 
maxima for 


w'l 


, 


nm 
= > + mx, (13) 
where m=0,1,2,... 

Consequently, the resonance frequencies will be 
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ra c 
a, = (0.5 + m) 57. (14) 


With an increase in oscillation frequency, ch + in 


(11) will also increase, due to which the amplitude re- 
sponse will have its greatest value in the first resonance 
maximum, The amplitude response in the subsequent 
resonance maxima will gradually decrease, 

Figure 1 shows the amplitude response A (w) and the 
phase response ¥ (w) curves, which were calculated ac- 
cording to Eqs, (11) and (12) for the following case; 


l = 1000 cm» v = 1.0 cm*/sec, 89 = 
= 0.1 cm, p = 0,8 886", d= 


=0,4cm, E = 0,8-10*% kg/cm’, Ko = 14000 kg/cm?, 
It is obvious from Fig. 1 
that, for frequencies of up to 200 cps, the amplitude re- 
sponse has three resonance maxima: at frequencies of 
10, 73, and 134 cps, At the resonance maxima, the 
3 

phase-shift angle is equal to = 5 T, >, respectively. 

For the case with liquid discharge, by using bound- 
ary conditions (5) and (6), we obtain: A=pp»,n =i (w'-ia") 
and 


io 
rk Sin ne + M cos nz 





B= 
Pox ote 
ak cos nz — M sin nz 


By substituting the values of A, B, and n in (7), we 
obtain after transformations; 
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ota Pot 
a= ol  a’'l : 
cos ch — B (@P — a,Q) + i 
where MK 
P=: 
V o + 4a? + @ V o + 4a? — 
a= x ,_o. = om ° 


P =sin@! ch?! ,Q= cos! sin *'. 
c c c c 
If a is small in comparison with w, a,=1 andw,=0. 
Thus, in the given case, 


glut 
Pt = Por ig = Pod (w)etot—™!, (16) 
where the values of r and S are clear from expression@'} 

The amplitude response is determined by the equa- 
tion 























[sin sh << +B (w1.0 — a:P)| (15) 


i 
Verdes’ (17) 


and the phase response is given by the equation 


A(@) = 


¢ (@) = arc tg S. : (18) 


For low-viscosity liquids (for instance, kerosene), 
the value of a is small, In this case, the approximate 
equation for A(w) is of the following form: 








A(@) = : 


4-— sin? = (1 — Bp) (9) 
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It follows from (19) that, with an increase in 8 , the 
values of A(w) decrease. Figure 2 shows the A(w)=f(8) 
dependence for the first resonance maximum “ts = : 
whence it can be seen that the sharpest decrease in the 
amplitude response occurs for 8 =0-8, For 8> 8, no res- 
onance phenomena are observed with an increase in the 
oscillation amplitude. 

The dynamic characteristics of hydraulic lines were 
determined by experiments where the frequency method 
was used [3]. The pressure oscillations at the piping in- 
let were secured by means of a special stand, the dia- 
gram of which is shown in Fig. 3. The stand consists of 
the following components; a cylinder, 5, which is filled 
with the liquid, a piston, 4, a rod, 3, a weight, 2, and a 
centrifugal vibrator, 1, which is driven by an electric 
motor, The rod, the weight, and the vibrator are rigidly 
connected to each other. Their combined weight creates 
a certain initial pressure in the cylinder. As the vibrator 
eccentric masses rotate, a harmonic disturbing force is 
created, By acting on weight 2, this force causes it to 
vibrate, due to which harmonic pressure oscillations 
arise in the cylinder and, consequently, in the hydraulic 
line 9. Centrifugal LV-500 vibrators, which were made 
in the Bauman MVTU, were used for this stand. A 
2780-S “Orion” (Hungarian People's Rep.)two-beam pres- 
sure indicator, 6, with a piezoelectric quartz crystal was 
used for recording pressure oscillations at the hydraulic - 
line inlet and outlet, The data transmitters 7 and 8 of 
the pressure indicator have a high natural oscillation 
frequency, which attains 20 kc. Our experiments were 
performed in the frequency region from 0 to 160 cps, 
The average pressure was equal to 25-30 atm, and the 
total (double) amplitude varied from 6 to 20 atm. 

Figures 4 and 5 show the experimental dynamic fre- 
quency-response curves for loop copper hydraulic lines 
with a length of 1000 cm, an inside diameter of 4 mm, 
and an outside diameter of 6 mm, which were filled 


with transformer oil (Fig. 4) and kerosene (Fig. 5). The 
data in Figs. 4 and 5 indicate the presence of two res- 
onance maxima of the amplitude response. The phase- 
shift angle at these maxima is 7/2 and 3x/2 At the 
resonance maxima, the amplitude of pressure oscilla- 
tions increases 2,8- to 4.8-fold. 

For hydraulic lines filled with kerosene or trans- 
former oil, in correspondence with Eq. (11), the frequen- 
cy of the first resonance maximum for sin u/1/c* 1 
is approximately equal to 


Op = —_ (/ incm). 


The frequencies of the second and the third reso- 
nance maximum will be three and five times as large, 
respectively. Therefore, the theoretical amplitude re- 
sponses of the hydraulic lines under consideration must 
have resonance maxima for the 30, 90, and 150 cps fre- 
quencies, which is in good agreement with the experi- 
mental results, 

Thus, the dynamic frequency-response curves for 
hydraulic couplings were determined theoretically and 
experimentally,and their flow characteristics were ex- 
plained. 
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One possible approach to the analysis of matrix decoders made from semiconductor diodes and 


having a relay output is given. 


Decoders which use unicombinational matrices (C1 ) and which provide a defensive 
rejection (for any distortion of the code) are studied. 

Relations for a selection network which interrelates the parameters of the output relays, 
matrix, and power supply voltage are derived for both continuous and pulse powering. 


The industrial development of semiconductors, in 
particular, diodes with a large back resistance, perrnits 
one to build and use diode matrix circuits with relay 
output in telecontrol and telesignaling (TC- TS) appa- 
ratus, 

The different types of such decoders are distin- 
guished from one another by the object selection code 
(uni- or all-combinational matrices are used), by the 
means of generating the input (a voltage input or an in- 
put through relay contacts is used), by the switching 
method, and by the operating mode of the output relays. 

The decoders mentioned, in contrast to the widely 
used matrix decoders with various kinds of trigger-cir- 
cuit outputs, require rather large currents of prolonged 
duration since their input resistances are relatively small. 

Moreover, a significant number of parallel circuits 
through the forward and back resistances of diodes is 
formed in unicombinational matrix decoders, depending 
on the selected code, and this raises the requirements on 
the quality of the diodes, and increases the power con- 
sumed, 

The matrix parameters, and sometimes the input 
voltages, will be different, depending on the required 
operating conditions and the relay parameters chosen. 

Therefore, it is best to be able to evaluate this or 
that type of decoder beforehand, and to compute its pa- 
rameters, 

One possible approach to the investigation of matrix 
decoders with relay outputs is given in this article, and 
relations are derived for decoders having defense in their 
circuits (defensive rejection for any sort of code distor- 
tion), and which use unicombinational matrices (CP), 
for two cases—continuous and pulse (through a capacitor) 
powering of the output relays. 

The effect of a change in the decoder parameters 
and of certain external factors on the operation is also 
examined, 
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Relations for the Circuits of Unicombina- 





tional Matrices (Ci) 





Let us examine a unicombinational matrix separate- 
ly, before deriving the relations for the decoder. 

A C2 matrix is shown in Fig. 1; the generalized pa- 
rameters for a unicombinational C™ matrix are desig- 
nated on it.’ It follows from an examination of this cir- 
cuit that the total number of input circuits is equal to 
the total number of combinational elements n; the num- 
ber of output circuits is equal to the number of possible 
combinations of the assumed N = C™; the number of in- 
put buses connected to one output bus is equal to the 
number of elements in each code m. The number of con- 
nections of one input bus to output buses s can be found 
from the formula s = cm : The number of output cir- 
cuits which are not connected to the given input circuit 
is equal to N-s, 


An Investigation of Decoders with Contin- 





uous Powering of the Output Relays 
Let us turn to the study of the decoder operation. 
According to the equivalent diagram (Fig. 2,c) for 
the circuit under consideration, we can write 





u=i(Ri + Reg), (1) 


vhere R, is the input resistance of the given decoder cir- 
cuit, Req is the equivalent resistance of the rest of its 
circuits, and i is the current of the equivalent circuit. 

The different operating modes of the decoder are 
caused by the change in the equivalent resistance Rog, 
and current i; the resistance R,, and applied voltage u, 
remain constant, 

Let us explain how the change in the above-men- 
tioned quantities is determined, 

Having used the intermediate equivalent circuit 
(Fig. 2, b), we get 









































R RR, Ry The resistances Rj and Ree in turn, are determined 
Req = Rhy + RAR, + KR, Tr R, Ry, * (2) from the following formulas, expressed in matrix param- 
i+n7 th eters: 
a 
SR R Ry (R fd + Ray) 
Resistances R, and R, will be equal to R zim) | rd RR,+R fd +Roip, 
a a r Ry (Rig + Ray) ie 
R, - (? fd +- R,) we ° (3) td + kRy+ R fd + Ray, + s(s— ) 2(m) 
i 
Ry=(R fa + Re) ——. (4) -~ 5 1 4 : 
5 Rom) * R Ry [Req + By yy | 
according to the equivalent circuit of Fig; 2, a, under the rd + kRy+ R fd +R, (x) 
condition that all elements with the same designation (5) 
are identical, 
Resistance Rg in the circuit is equal to R= { R 4 Ry (Rtg + R, (ep? 
c N—s—i rd kKRo+ Rey + Ray |, 
R R 
R, = —f4 + ace) : 
4 (6) 
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where Reg and R, are the forward and back resistances 
of the diodes; R, is the resistance of the relay coil; 
Rem) is the value of the matrix resistance R for the in- 
put circuits pertaining to the given code; Reg) is the 
value of resistance R, for the defensive circuits of the 
matrix; k is the number of defensive circuits. 

Let us find how a change in the resistance R, effects 
the equivalent resistance , all other quantities re- 
maining constant. It follows from formulas (2)-(5) that 
the equivalent resistance will equal 


R= 
eq R a ° 
1+—-+—F mM 


when Ram =* , where 








Ry ei Rym (s — 1) 
_ Ry (Reg + Faw) ’ 
Sa Peg RROt+ Keg + Ray 

Rp Ry (N —s — 1)(n — m) 





—-— —_— 
R Ry (Rea + Ray) 
b ms Pp“ fd 2(k) 

Reg Wet hig + ER + Rea + Ra 





When = = Reck) = 0, the equivalent resistance 
will be designated by Req’ which is expressed by the fol- 
lowing formula; 








ts R 
ea ae te 5 
R eq= ee 
itp 2 4 
’ (8) 
R 
— — p 
si Rim R, (N —s — i1)(n — m) 
P p 
‘hig * Ry Req 


Rfq (N—s—1)+Req + ERT Rea 


A numerical example is given below (Appendix I) 
to evaluate the effect of the second and third terms in 
the denominator of Reg. 

We note that Rag is determined with the output re- 
lay of the given circuit operated, and Re, is determined 
when that relay is idle. 

By examining the expressions introduced above, we 
can draw the following conclusions, 

1) For a given choice of diodes and relays, the 
equivalent resistance and current in the circuit vary asa 
function of Re. The largest value of R,q occurs when 
Rg = ©. 

2) The back resistance of the diodes substantially 
effects the value of R,,,, decreasing it significantly when 
Req is small, and very little when R,g is large. 

The absolute value of resistances R, and R,, which 
are determined from formulas involving the back resist- 
ances of diodes, depends on the number of elements in 
the code, the total number of elements, the number of 





output circuits, and the number of connections between 
the buses, that is, on the chosen C™, 

3) If it is assumed that the back resistance is suf- 
ficiently large (Rpg * ©) and the matrix resistance 
Rem) = «©, the equivalent resistance will be equal tothe 
relay resistance, Rag = Rp- The resistance, Rem) = @, is 
characteristic for decoder inputs coming from relay con- 
tacts; in all other cases, including voltage input (if Rem) 
is not infinite), Req will be less than the value men- 
tioned. 

4) The equivalent resistance R&q will equal 








R ae p 
eq Rym . 
1+ 
Rtg 


for the same diode back resistance (Rpg = ©), and for 
4 1 
Ram) = Roxx) = 0( Ra =~ Rta =), that is, it will be 


m 


significantly less than Reg. 

We will derive the design formulas relating the de- 
coder parameters to the applied voltage. 

When the given voltage is applied, the decoder 
must ensure reliable operation of the output relay of the 
circuit which corresponds to the code, and must ensure 
that the rest of the output circuits do not operate. 

Moreover, the decoder, according to the condition 
imposed above, must ensure a defensive rejection for 
any Code distortion, that is, if there are more or less 
than the proper number of code elements m. 

The latter condition is equivalent to stating that 
the output relay must operate only when all defensive 
circuits are disconnected (Rem) = @), and must not op- 
erate in any other case, including the case in which 
only one defensive circuit is energized. 

Reliable operation of the relay will be ensured by a 
suitable choice of the matrix input resistance, R,, forthe 
condition, Rem) = @, and for a given relay pull-in re- 
serve factor. 


One can write 


i (Rit Req) = Rai’ + fogs Copltp=4, ©) 





for this case, where igp is the operating (pull-in) cur- 
rent of the relay; krs(op) is the pull-in reserve factor; 
i’ is the current flowing through R, in the mode indi- 
cated, equal to 


i’ = lop rs(op) + tpg [8 — 1) m + 
+ (N —s—i)(n—~m)}. 


Substituting the current i* in equation (9), and 
solving for Ry, we get 





R,-; 1 eee? is 
‘op *rs(op) ting (m(s—1) + (n—my(N—s— 1)] 








If the back resistance of the diodes is large, the 
second term in the denominator can be ignored, and the 
approximate formula 


u 


R,=z-— — — ft,. 
‘op Kes(op) . 





can be used, 

This formula, expressed in terms of the drop-out 
current and reset coefficient of the relay, will take the 
form 


where 


The reverse current which flows through the diodes 
can be found from the formula, 


PS ep els 8 op ‘rs(op) (y+ Mr) 
SS Ri m(s —i)+(2— m)(N —s— 1)}° 





By using the relations for the second operating state 
of the decoder (relay unoperated, where Rem) = 9, and 
Red = ©, R&q), we find the forward current of the diodes, 
This current, ifg¢,), depends on the number of defensive 
circuits formed, and can be determined from the ex- 
pression 
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where k=1,2,...,m,...,m is the number of de- 
fensive circuits; the remaining designations are the same 
as before. 

The largest current through the diodes will occur 
when there is one defensive circuit. It is useful to know 
this current in order to choose diodes correctly, so that 
they will not be overloaded during prolonged operation, 

The power required by the matrix decoder depends 
on its operating mode, and a large amount of power will 
be required if no code is present, 

This power can be determined from the formula 


. u2N 
P = ig 8 = ———, 


on the basis of the equivalent circuit (Fig. 2, c), under 
the condition that all the circuits are identical; Rég is 
the equivalent resistance in the absence of code, and 
itor is the total current required by the decoder. 


An Investigation of a Decoder with Pulse 
Powering of the Output Relays 

A matrix decoder, whose output relays are connect- 
ed through capacitors, is shown in Fig. 3. 

An examination of the equivalent circuit (Fig. 4, a) 
shows that the given circuit has parallel circuits of two 
types; a) circuits which contain inductances and capaci- 
tances connected to them through the back resistances 
of diodes, and b) circuits which are connected through 
the forward resistance of diodes, and contain the resist- 
ance Rg. 

An analysis of these circuits shows that the first 
type does not have much effect on the operation of the 
given circuit if the back resistance of the diodes is suf- 
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ficiently high, since the transient and steady-state cur- 
rent are basically determined by Ryg. 

Therefore, all these circuits, for practical purposes, 
can be replaced by an equivalent resistance expressed in 
terms of the back resistance of a diode. 

The second type of circuit contains the resistance 
Re, which, as in the case examined earlier, determines 
the operating mode of the relay, and may change from 
0 to @, 

The circuits mentioned are shown in Fig. 4, b, where 


Rrd 


R= (n —m)(N —s— 1) +-m(s — 1)° 





The computational equivalent circuit is shown in 
Fig. 4, c, where Req will assume different values which 
are a function of the operating mode. When R, = ©, the 
equivalent resistance will equal 





. R 
nen rd 
Req= (n — m)(N —s — 1)+ m(s —1) Cn 
and when R, = 0 
° R . 
R eg= —_. (11) 


We will derive relations for the given circuit, as- 
suming that all quantities entering them are constant. 
We can write the equation 


. 
u=(ie+ is) Rit icky +Ly seu, (12) 


from the computational equivalent circuit. 

In this equation, the current i, flowing through the 
given circuit must not be less than the relay operating 
current, when the relay is operated, 


t= lopkes(op) 




















but must be less than the relay drop-out current for all 
the remaining cases: 





fe ido 

= =< 
° — ¥ ts(do) 
where kys(ony > 1, and do > 1, 


Manipulating equation (12) further, we obtain a 
secon¢ order differential equation 


c 


dtu RyR+ Ri du, ‘ 
PT RL, “at Lc Ue 





= U, (13) 


where R = (rez + 1), Lp is the inductance of the 
eq 


relay coil; Rp, the relay resistance; Ry, the matrix input 
resistance; C, the capacitance of the isolating capacitors; 
and u, the applied voltage. 

The roots of the characteristic equation correspond- 
ing to equation (13) are 
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Integrating equation (13), we find the current 


through the relay and the voltage across the relay coils 
and capacitor: 








uR, { 
ne ™ Bit Req [1 ~ ta 1 (ase — aye%)| (15) 
gs uReg 
1 ~ (Ri + Reg) by (aa — a) (em! —e*#), (16) 
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Substituting the values of the equivalent resistances 
for the various modes in expressions (15)-(17), we find 
the currents and voltages corresponding to them in the 
circuit being studied, A numerical example is given in 
Appendix II. 




















CONCLUSIONS 


One can draw the following conclusions from the 
given analysis of the operation of the decoders men- 
tioned above, and also from experimental work on them 
done by the author. 

1, The choice of decoder matrix parameters (pa- 
rameters of diodes, relays, etc.), and of the supply volt- 
age for continuous and pulse powering of the latter,de- 
pends on the type of matrix (N = CP), that is, the total 
number of output circuits, combinational elements, code 
elements, number of defensive circuits, and number of 
connections between the input and output circuits, 

2, The operating mode of the decoder~a reliable 
selection of the output relay corresponding to the given 
code, and a defensive rejection if this code is distorted 
in any way~is basically determined by the matrix resist- 
ance Rg. The best results in the given two cases are ob- 
tained when this resistance changes from 0 to @, 

3. The back resistance of the matrix diodes sub- 
stantially affects the operation of the decoder output re- 
lays, and this effect increases as the total number of ma- 
trix output circuits grows. 

A decrease in the diode back resistance or an in- 
crease in the total number of output circuits leads to a 
decrease in the relay operating reserve factor when a 
code is present, and in the drop-out reserve factor when 
no code is present, if the voltage is fixed. In a number 
of cases, this circumstance makes it necessary to in- 
crease the supply voltage. 

4, The energy content of matrix decoders with re- 
lay output is relatively large, and is determined by the 
type of matrix and the output relay parameters, and also 
depends on the operating mode of the decoder. A large 
amount of power is required in the absence cf code, 
both for continuous and for pulse-powering. 

Therefore, it is best to connect this type of decoder 
to supplies of low power only during code reception, 
even if there are only an average number of output cir- 
cuits (15)-(20), 

Moreover, if this is done, one can lower the require- 
ments on the ‘permissible amplitude of diode forward 
current when the defensive circuits are formed, since 
the diodes can take a significant overload for a short 
time. 


APPENDIX I 
Numerical Example. Let us assume that a C2 ma- 
trix is used in the decoder, and that it is made from 
type D7Zh diodes, and that type RKN telephone relays 
are used in the output circuits, In this case, the basic 
data will be 
R, = 1000 ohm, Rp, = 1Gohm, 4 





= 5-10%ohm, = 2, 

N=10, n=5, m=2, s=4, N—s—1i=5, 
Raiyy= % Racmy = % (for Req )s 

Racy = opm) = O (for Red 





Computations give the following results; 





° Ry Ry 
Reqa= TF 0.011 + 0,08 = 1.04 =996 ohm, 


R" Ry Ry 4.97 ohm 
eq” 1+200400 ~ 01.0 -* . 





APPENDIX II 


As an example, let us find the current flowing 
through a relay for the case computed in Appendix I, as- 
suming that the relay coil inductance is L,, = 10 h, and 
the capacitance C = 16 pf; all the remaining quantities 
are the same as in the preceding example. 

The operating and drop-out currents of the relay 
are, respectively, ign = 20 ma, ido = 9 ma, Rg=500 ohms, 

When the relay is operated, we get 

Red 
Reg = (a—m)N —s— 1) + mis —1) ~ 


500 - 10? 





When the relay is idle, we get 


P Reg 10 


Req = m =F =ohm 


Let us find the resistance R for these two cases: 





R 

R' = —>— + 1 = 1,021 ohm, 
R’ 

Rom or + 1 = 10 ohm. 
eq 


The roots of the characteristic equation for the re- 
lay in the operated state are 


‘and for the relay idie- 
a, = — 34,6; a, = — 66. 


The current through the relay while it is operated 
is found from the expression 


‘ uReq 
"(Rit Reqly(a,— 2,) 


Since the second term in the parenthesis is signifi- 
cantly less than unity when t = op)» it can be 
neglected, and then we find that i~ 40 ma. 

Since the operating current of the relay is equal to 
iop = 20 ma, an operating current reserve factor, 

= 2, is ensured, 
“top idle relay current is found from the formula 
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Rr ~ - 
a eq (et Scat eat), 





~ (Ret Req ky (a, — @,) 


When t = tgo, we find that i* » 0,46 ma, which is 
significantly less than the operating current. 
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ON THE INTERFERENCE RESISTANCE OF PULSE-TIME 
AND PULSE-WIDTH TELEMETERING SUBJECTED 
TO STRONG FLUCTUATION NOISE 


N. V. Pozin 
Moscow 


Translated from Avtomatika i Telemekhanika, Vol. 21, No. 9, pp. 1320-1322, 


September, 1960 
Original article submitted January 30, 1960 


The possibility is discussed of raising the interference resistance of pulse-time and pulse-width 
telemetering by introducing a time selector into the receiver; the selector would eliminate the 
effect of short duration, but highly probable, errors formed by the action of fluctuation noise. 


Video pulse errors are formed in the video channel 
of a receiver by the action of strong fluctuation noise. 
In general, these errors are false pulses and false spaces 
of different duration.° 

The influence of strong fluctuation noise on the 
transmission of PWM and PTM telemetering is analyzed 
in reference [1]. It is shown that the average and mean- 
square errors are determined by the probabilities, P; and 
P.., of the unit errors of the space and pulse. The unit 
error is defined as the formation, by noise, of a false 
pulse or false space of duration At = 1/W = 1/2f,, where 
f, = W/2 is the cut-off frequency of the low-pass filter 
in the video channel, which is replaced by a band-pass 
filter if there is no secondary modulation before trans- 
mission, 

It is assumed in [1] that the receiver is sensitive to 
errors of duration At. 

Is this really necessary? It is true that an input 
filter with pass-band W creates the necessary conditions 
for passing pulses At wide. But the pass-band of the in- 
put filer is designed to be wide enough to pass working 
pulses with a sufficiently fast rise-time. If transmission 
stability requirements oblige one to pass pulses with a 
fast rise-time, this requirement does not extend to the 
width of those pulses, 
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Fig. 1, Pulses in the video channel of a receiver. 
a) before the time selector; b) after the time selec- 
tor (ty <t< te) 
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By taking special measures, such as introducing cir- 
cuits in the video channel of the telemetering receiver 
which delay time by t, and which make the receiver in- 
sensitive to short-duration errors, one can significantly 
increase the interference resistance of the receiver. 
This remark is intended to draw attention to this simple 
method, and to make the necessary additions to the 
analysis of interference resistance in reference [1],where 
it is assumed the’ main influence on the indications of 
the receiver is the extremely short, but highly probable, 
errors of unit duration, At. 

Ideally, such a device, which can be called a time 
selector (filter), or duration selector, is triggered by the 
leading edges of the incoming pulses. Pulses are gener- 
ated at the output which are time-shifted by t, and 
which duplicate the input channel, providing it does not 
contain pulses or spaces shorter than t. Pulses shorter 
than t do not trigger the device, and likewise, spaces 
shorter than t cannot restore it to the initial state. Thus, 
the time selector does not pass signals which are less in 
duration than some set value, t (Fig. 1). 

This type of apparatus can be made in different 
ways, The simplest is a relay which is slow to operate 
and release, and which is actuated by constant ampli- 
tude pulses, In general, if there is double limiting in 
the video channel (Fig. 2), a low-pass filter can be used 
as a time selector (for instance, an integrating RC cir- 
cuit), if the method of demodulation after this is such 


* Telemetering subjected to strong noise is ordinarily ac- 
complished by transmission in which the rise time of the 
working pulses is very short, The permissible width of 
the leading edges in this case is of the order of the inde- 
terminateness permissible in measuring time intervals, 
and distortions of the leading edges by noise can be neg- 
lected, 
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that the “piling” of the leading edges of the working 
pulses in the filter does not lower the measurement ac- 
curacy. If this is not so, a more complex time selector 
must be used . 

As is mentioned in [1], the duration of the unit er- 
ror is the correlation interval of the smooth noise at the 
filter output, and thus the formation of unit errors can be 
considered independent. If the probabilities of unit er- 
rors (that is, errors of duration At) in a space, Py and 
pulse, P,, are known, then the probabilities of errors 
which have a duration equal to the delay time, t, are 
found in the following manner: 











P(t), = Pyoand P(t) = Fe (1) 


It is certain, from the standpoint of interference re- 
sistance, that the greater the delay, 





t =kAt, (2) 





the better~that is, the less—the error due to noise. 

However, the delay cannot be too great; it is lim- 
ited by the following factors, 

First of all, an increase in the delay causes an in- 
crease in the minimum values of the pulses and spaces, 
whose duration must exceed the delay introduced by the 
time selector. As a result, the useful time deviation 
which carries information is decreased somewhat, and 
the telemetering accuracy is reduced, or the rapidity of 
response is lowered as the cycle duration is correspond- 
ingly increased, But the magnitude of these factors is 
not great, and it is seen that in practical devices, they 
do not limit the delay. 

The second most important factor limiting an in- 
crease in the delay is the instability of the device which 
introduces this delay. 

The greater the delay, the more difficult it is to 
meet the stability requirements, Hence, it is evident 
that a relativelv unstable, but small, delay introduces a 
large error in the indications of the telemetering re- 

























A telemetering receiver insensitive to short errors 
of duration At would be of great practical value. To a 
lesser degree, it is also desirable to eliminate the less 
probable errors of duration 2At, etc, If the duration, 
At, is comparable to the rise time of the pulse, that is, 
to the telemetering accuracy under conditions of strong 
interference, it is not difficult, in practice, to ensure a 
sufficiently stable delay, t = (2~3)At, that is, a coeffi- 
cient k = 2-3, 

In order to evaluate the interference resistance of 
PWM and PTM telemetering apparatus,incorporating a 
time selector with delay t(2) in the receiver, the proba- 
bilities P; and Pp in the design formulas derived in [1] 
must be replaced by PK+! and pPk+! respectively. 

Simplified expressions for the average (5 ay), and 
mean square (5 ms) reduced errors are given in the fol- 
lowing table;t 











Method of Without delay With delay 
modulation + 2 
Say Sms Say oms 
tsie & | h pk+t p+! (+1) 
2 2Nn 2 2Nn 
bres nPs nP, | npkt? npktt 
6 12N | 6(k+1) | 1L2N(k+1) 














Tin practice, the ideal time selector is based on a scaler 
circuit, triggered and blocked by the réspective leading 
edges of pulses, 

tWe must keep in mind that the formulas presented in 
the table do not take into account possible error at the 
edges of the pulses and spaces, The probability of these 
errors is low, but the error component due to “edge” er- 
rors will increase as the delay is increased. 
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time selector, F-band-pass filter, 


Fig. 2, Block-diagram of a telemetering receiver using PWM and a 


D-amplitude detector, L~limiter, 


TS-time selector, DM-demodulator-this transforms pulse duration 
into proportional amplitudes of voltage or current, 














SUMMARY 

The introduction of duration selection in the video 
channel of a receiver significantly increases the inter- 
ference resistance of PWM and PTM telemetering. Such 
a method is especially valuable in telemetering appara- 
tus working under conditions of strong interference, 
since there is no strict relation between the durations of 
the pulse, and of the leading edge: the pass-band of the 
filter determines the rise time, but the duration of the 
shortest pulses can vary over wide limits, 

The introduction of long delay, two to three times 
more than the duration of a unit error, At, is not diffi- 
cult to realize, and can be used practically in all time- 


type telemetering apparatuses, In connection with this, 
it is best to examine the data of reference [1] on the 
comparative interference resistance characteristics of 
telemetering devices which use different modulation 
methods, Such a comparative analysis will be the sub- 
ject of a special article. 
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A form of the differential equation for heating a thermistor by the current flowing through it is 
proposed; with it, one can obtain all the variables necessary for calculating the transient directly 


from a graph of the volt-ampere characteristic. 


A graphical method of computing the transients 
arising in a thermistor circuit (TC) excited by a step 
change in the parameters which characterize the oper- 
ating mode of the circuit is presented in [1], and is fur- 
ther developed in[2, 3]. This method has a great disad- 
vantage: the thermal time constant, r, whose magnitude 
greatly depends on temperature, enters the calculation, 
Assuming + to be constant results in significant error in 
caiculating the transients, which are accompanied by a 
temperature change of hundreds of degrees in the work- 
ing body of the TC. This disadvantage can be elimi- 
nated by a method suggested in [4], where a form of the 
differential equation for the heating of a TC is given in 
which the thermal time constant is eliminated from the 
calculation. However, the solution of the problem in 
this case remains quite complicated. 

Still another form of the differential equation is sug- 
gested below; it permits one to curtail the amount of 
computation significantly, and makes the calculation 
quite graphic. 

The following form of the heating differential equa- 
tion is recommended in [1] for convenience in calculat- 
ing the rate of temperature change in the working body 
of a TC at various temperatures: 


or GR. 2-9, 
dt ~ H(R+re” =< ° 





(1) 
where T is the instantaneous temperature of the TC 
working body, in °K;. R, the instantaneous resistance of 
the TC,ohms; H, the heat capacity of the TC, watt-sec/ 
/deg;r, the linear resistance in the circuit, ohms; 
Up, the voltage applied to the circuit,v; r, the thermal 
time constant, sec; and T», the ambient temperature, 
x. 

The power liberated in the TC at any time during 
the transient can be found from the formula, 
u2 
Tee RU 


where U" is the instantaneous voltage drop across the TC, 
and I’ is the instantaneous current, 
The time constant is defined by the expression 

t = H/b, where b is the dissipation coefficient. Taking 
this into account, we can rewrite expression (1) as 

aT U'I’—b(T —T,) 

ae if 

The product of the dissipation coefficient b and the 

temperature of overheating (T—T,), which appears in 
the above expression, is equal to the power dissipated by 
the TC into the surrounding medium under steady -state. 
conditions at a temperature of the working body T, i.e., 





b(T —T,) = UI, 


where U is the voltage drop across the TC in the steady- 
state, and I is the steady-state current through the TC. 
Consequently, 


aT U'tt'—UI 
 Seeee kunt (2) 

Turning to the volt-ampere characteristic on which 
the load line (Up-iy is drawn (Fig. 1,4), it is easy to see 
that all the quantities dependent on temperature which 
are needed to find the TC temperature change rate can 
be taken directly from the figure; these are the coordi- 
nates of points K* and K, lying on the resistance line 
which corresponds to a given TC temperature, T, One 
of these points belongs to the dynamic volt-ampere char- 
acteristic, and the other, to the static. We will call the 
dynamic volt-atnpere characteristic the geometric locus 
of points, whose coordinates are equal to the voltage 
drops across the TC and the current through it at any 
time during the transient. In the given case, the dynam- 
ic volt-ampere characteristic is the segment 1'-2, of the 
load line, The heat capacity can be found by the meth- 
od suggested in [4]. 

Let us illustrate the technique of computing a tran- 
sient by a concrete example, It is necessary to plot a 
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curve of the current increase in a circuit consisting of a 
type KMT-11 TC and a linear resistance, f, connected 
in series, The transient arises when r is changed from 

ty = 180 kohms to tr, = 12,65 kohms in one step, The TC 
is immersed in transformer oil at 20°C. The volt- 
age, U,, applied to the circuit is 220 v, The tempera- 
ture characteristic of the TC is described by the expres- 


4600 
sion # =:0,281e 7, and the heat capacity, H= 2.78 mw- 


sec/deg. The volt-ampere characteristic of the TC, 
with resistance lines corresponding to eleven different 
temperatures, is shown in Fig. 1, a 


We begin the calculation by determining the steady- 
state modes of the TC before the start and after the com- 


pletion of the transient. Load lines, U,~Iry, and U,—Ir,, 
must be drawn for this purpose. Point 1 characterizes 
the state of the TC before the transient starts (R=860 
kohms, T=309°K). After the transient is completed, 














the state characterized by point 2 is established (R=4.4 
kohms, T=474°K). 

The temperature change rates of the TC at differ- 
ent temperatures occuring during the transient are cal- 
culated from formula (2). For example, when T=408° 
K, the point K’ on the dynamic volt-ampere char- 
acteristic has the coordinates: U* = 134 vy, I'=6.2 ma, 
and point K on the static characteristic has coordinates 
U= 88 v, I= 4 ma, and consequently, 

. e 2 88-4 ny ities. 

The rates found by this method are plotted on the 
graph dT/dt = f(T), for the corresponding temperatures 
(Fig. 1, b) 

Then a curve of the temperature change as a func- 
tion of time is plotted(Fig. 1, c). To do this, the entire 
temperature change interval is divided into a series of 
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Fig.1. An example of the calculation of a transient, a) initial volt-ampere character- 
istic; b) curve of the dependence of the TC working body temperature change rate on 
temperature; c) temperature change curve; d) current change curve, 














segments, AT (16 parts, in the given example) within 
which the temperature is assumed constant, and the cor- 
responding time intervals are computed, Thus, when the 
rate is dT/dt = 9 deg/sec, a temperature increment 

AT = 11 deg (interval 1) occurs during the time At= 
11:9 = 1,22 sec. By adding all the time intervals 
At, we construct the curve T = f(t). 

The construction of the curve I = f(t) is clear from 
the drawing (Fig. 1, d). 

The curve of current change obtained by calcula- 
tion almost coincides with the experimental results (the 
dashed curve in Fig. 1, d). The experimental curve, like 
the initial data for the example under consideration, was 
taken from [5]. 

Transients arising as a result of a step change in the 
supply voltage or a change in heat-transfer conditions are 
calculated analogously. The volt-ampere characteristic 
needed for the calculation can be found analytically [6]. 
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THE SEVENTH SCIENTIFIC -TECHNICAL CONFERENCE OF 
YOUNG SCHOLARS OF THE INSTITUTE OF AUTOMATION 

AND REMOTE CONTROL OF THE ACADEMY OF SCIENCES 
OF THE USSR 


On March 14-16, 1960, was held the traditional 
(seventh) Scientific- Technical Conference of Young 
Scholars of the Institute of Automation and: Remote Con- 
trol (IARC) of the Academy of Sciences of the USSR on 
Problems of Automatic Control, The conference discussed 
the work of young scholars which was completed in 1959. 

More than 400 persons participated in the work of 
the conference, including about 200 representatives of 
various organizations in Moscow and the Moscow area 
(scientific research institutes, institutes of higher learn- 
ing, planning and design organizations, etc.). 

Seventy-five reports and communications were heard, 
encompassing a broad range of problems of automatic 
control, The materials presented evoked a lively ex- 
change of opinions. 

In opening the conference the director of the Insti- 
tute of Automation and Remote Control, Academician 
V. A. Trapeznikov, wished the young scholars success in 
their work and charged them to link their scientific work 
more closely with the practical, with the concrete needs 
of the national economy. At the plenary session Profes- 
sor M. A, Aizerman, D. T. S., presented a report on “The 
scientific problem of the theory of finite automata." 

At the final plenary session £, VY, Voloshina and 
E, V, Shtil'man presented a report on “Modeling the 
processes of instruction.” In the report they reviewed the 
charactetistic features peculiar to the processes of in- 
struction of animals and man, surveyed a number of auto- 
matic devices which fulfill certain functions of the hu- 
man brain,and showed analogous features between certain 
processes in automatic systems and processes of instruc- 
tion of people and animals, Between the plenary sessions, 
the work of the conference was done in sections; 

I. Automatic control, with subsections on the the- 
ory of automatic control and systems of automatic control 
II. Automatic inspection 
III, Computers 
IV. Elements and devices ofautomation and remote 
control 
V. Statistical methods in automation 
VI. Theory of relay circuits and finite automata 
VII. Automated electric drive 


I. Automatic Control Section 
A. Subsection on the of automatic control 


A report by V._N. Novosel’tsev was devoted to find- 
ing an expression for optional control in secondary relay- 





impulse systems for the case of purely relay control and 
for the case of relay control in the presence of a zone of 
insensitivity. 

I, S. Morosanov discussed in his report the effect of 
fluctuations on relay extremal systems under auto-oscil- 
lating conditions, On the basis of evaluation the princi- 
ples of harmonic and statistical linearization were as- 
sumed, The author showed that the methods of calcu- 
lating statistical transfer coefficients in the form pro- 
posed by I. E. Kazakov are inapplicable in this case, But 
a linearization equivalent in probable meaning can be 
obtained in reviewing the problems of parametric effect 
of fluctuations, The final solution for the case of linear 
members of an extreme object of any order was given in 
graphic form. As a result, on the basis of published sta- 
tistical characteristics of the fluctuations, the mean value 
of the frequency and amplitude of the auto-oscillations is 
found for a given system of extreme control. 

The joint report of V. G. Gradetskii and Yu, L 
Ostrovskii entitled "The operation of a system of ex- 
treme control (SEC) with memory of the extreme under 
noise conditions “ was devoted to the problem of raising 
the noiseproof features of systems, The value used in 
the report as an index of the noiseproof feature charac- 
terizes the number of false switchings of the executive 
organ in a unit of time. 

The task amounts to a determination of those values 
of the parameters of tuning of the regulator and noise 
filter which, at a given rate of scanning and a given loss 
in scanning, would provide the minimal amount of false 
switchings in a unit of time. This problem was solved 
approximately on account of its complexity. The for- 
mulas obtained have been verified by experiment, 

A report by N. V. Grishko presented the results of 
determination of the optimal characteristics of a single 
extreme system during occasional functioning. The sys- 
tem reviewed is a natural generalization of several very 
noiseproof and high-grade extreme systems. The opti- 
mal characteristics are determined by varied methods by 
minimization of the analytical expression obtained for 
the mean square error through the given characteristics 
of the object being controlled and the characteristics 
produced by the extreme system. 

T. G. Babunashvili devoted his report to a study of 
the triggering process of a regulating system containing 
a servomotor with a nonlinear velocity characteristic and 
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fixed feedback. The task consisted of investigating sys- 
tems with three differential first-order equations, with co- 
efficients variable in time and one nonlinearity from the 
point of view of dissipativity, largely in accordance 

with the thinking of Kh, Masser. The author of the re- 
port determined the maximal deviations and the devi- 
ations which are determined at the end of the triggering 
process, 

In the report of A. I. Moroz entitled “Qualitative 
investigation of differential equations obtained in the 
solution of the problem of synthesis," a method of de- 
signing regulating systems was explained, It was sug- 
gested by A. M. Letov and consisted of the fact that the 
equation of the regulator was defined by the task of the 
first integral of the system of differential equations of 
the object and the unknown regulator. In an occasional 
case the advantages and deficiencies of this method 
have been noted, and certain unacceptable features of 
the systems obtained have been disclosed, The report 
proposed a method of eliminating the deficiencies dis- 
covered, 

A report by M. M, Simkin was devoted to the de- 
termination of periodic conditions in nonlinear impulse 
systems, The indicated problem is linked with a review 
of nonautonomous systems, 

The author proposed an approximate frequency 
method of determining simple and complex periodic 
conditions in such systems, 

In a report of R. P. Parsheva the problem of stabil- 
ity in the Lagrangian sense of transient conditions in a 
peutametric SAR with a single nonlinearity was discussed. 

V. S. Kislyakov, in a report on “The prolonged sta- 
bility of an aircraft with an automatic pilot possessing 
retardation," investigated the stability of short- period 
motion of an airplane with an automatic pilot without 
feedback, Disturbed movement of the system is de- 
scribed by third-order differential equations with con- 
stant coefficients and retardation in the characteristic of 
the steering gear. To determine the region of asymp- 
totic stability as a function of retardation and adjust- 
ments of the automatic pilot, on the assumption that the 
characteristic of the steering gear has weak nonlinearity, 
the method of asymptotic approximation of N. M, Krylov 
and N, N. Bogolyubov was used; this was also utilized 
in studying the linear problem. The effect of retardation 
on stability of movement was explained. 

In the report of B, E, Chuprun on “The formulation 
of a law of control of linear stable objects assuring mini- 
mal displacement of the regulating organ,” a method 
was studied of controlling a constantrate servomotor, in 
which the regulating organ deviates only to the amount 
needed for compensation for the displacement, It is as- 
sumed here that the form of displacement is known.but 
its level is not. In like manner a study was made of a 
case ofa linear stable object of arbitrary order without 
and with retardation. 











B. Subsection on systems of automatic control 

B. B, Buyanov discussed in his report the application 
of the theory of systems optimal for high speed to the 
problem of controlling the flying shears of a rolling mill. 

B. G, Volik reported on an automatic optimizer 
with two channels and two limitings intended to make 
scannings of the extreme directly on the object which has 
great inertness. The problem was studied with regard to 
the automation of chemical objects. I. N. Bocharov re- 
ported on an original instrument for photographing the 
distribution curve of a random value. 

E. A. Rateeva reported on a three-channel optimizer 
for chemical production, calculated for 2-3 variables and 
two limiters. 

In the report of V. G. Sholokhov, experimental proof 
was given for convergence ‘n the process of adjusting a 
noise generator, by the filter-forming method, for the 
case of two filters,and the dependence of the time of the 
transitional period on the method of triggering was deter- 
mined, 

The report of K. B. Norkin was devoted to the prob- 
lem of automatic adjustment of a transmitter‘s output 
cascade with use of a system of automatic triggering, and 
contained an explanation of the method of optimal selec- 
tion of adjustments of the optimizator during operation 
with an inert object under noise conditions. 

In a report of V. N. Shadrin a system of program con- 
trol with frequency separation of the channels was de- 
scribed, 

Two reports made by V, V. Karibskii and A. P. 
Evseeva contained data on a universal interpolator for a 
system of digital program control and on the automatic 
selection of interpolation segments for a machine with a 
linear interpolator, 

A report by M. L Taits described an experimental- 
industrial model of an instrument for automatic program 
installation, Reports of B. N. Andreichikov dealt with 
problems of the dynamic precision of machines with pro- 
gram control and the effects of friction and slack on the 
error during reverse operation. 

The report of M, M, Khasanoy was devoted to an 
analysis of the dynamic characteristics of air-condition- 
ing equipment, 

Il. Automatic Inspection Section 

M. V. Rybashov and I. M. Ponasenko devoted their 
report to a study of several circuits for the dynamic com- 
pensation of counters, The conclusions obtained were 
confirmed by experimental investigations on models, 

V. _S, Likhoninskii, in his communication "A capac- 
itance- measuring device for digital systems of program 
control of productive processes," examined the operat- 
ing principles of capacitance-measuring devices of the 
linear and rotary types for the measurement of shifts, 

M._ A, Prusov reported on the principles of building 
and designing an instrument for measuring the tempera- 
ture of rotating parts, The instrument is based on utili- 
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zation of turbulent currents and has a number of advan- 
tages in comparison with analogous instruments based on 
other principles, 

The report of Yu, V, Gushchin was devoted to study 
of the possibility of using semiconductor crystal layers 
on a sulfur-cadmium base for the indication of radio- 
active radiations, 

A. A, Kalmakov discussed in his report the possibil- 
ity of using radioactive isotopes and the methods of nu- 
clear spectroscopy for automatic control of the content 
of metals in alloys and products of the enrichments of 
ores of nonferrous metals. 

V. A. Viktorov reviewed the effect of the high har- 
monics of a coaxial vibrator on the operation of an ex- 
tremal endovibratory level indicator. 








Ill. Computer Section 

In a report by A, G, Butkovskii, "The modeling of 
certain objects with distributed parameters," the theo- 
retical principles and a structural scheme are given for 
an electronic model of straightaway heat-treatment fur- 
naces operating in an automatic-control system. 

A. V. Shileiko, in his report on “Method of synthesis 
of an optimal structure of a digital-model,” examined 
problems linked with the selection of criteria of quality 
and method of synthesis of an optimal structure of a dig- 
ital model,for which criterions of quality he assumes ex- 
treme values, 

In a report of R. N. Chernyshev, “A solver amplifier 
with powerful output cascade,” the results were given of 
the development of a solver amplifier, completed ac- 
cording to an autocompensating scheme, The amplifier 
has a contactless modulator on silicon diodes and a pow- 
erful output cascade, 

In a report of B. A. Pereverzev entitled "Combined 
electromechanical assembly of an electronic modeling 
device" there were descriptions of the circuits, design, 
and basic technical data of an electromechanical unit 
which can operate under the conditions of a multiplier- 

















divider, a functional converter and a coefficient buncher. 


F. B, Gul'ko, in a review report entitled "High- 
speed electronic multiplying devices,” furnished data on 
new electronic multiplying devices with an elevated 
transmission band, their characteristics, and the factors 
limiting high speed, 

Zh, A, Novosel'tseva, in a report entitled "Control- 
led retardation assembly,” gave the results of the devel- 
opment of a device for the modeling of retardation, the 
time of which changes in accordance with the control- 
ling signal, The device was made for memory condens- 
ers which are switched with the help of step-by-step 
switches, 

In a report of K. B. Norkin, "On a single method of 
automatic triggering of the extreme of a function of 
some variables," a method is given for discovering the 
extreme of a function which has lumps and saddles which 
do not require calculation of the gradient of the function. 








V, A. Yakovlev, in a communication entitled "A 
discrete electrical differentiator," described the oper- 
ating principle and design of a device for the differenti- 
ation of slowly changing voltages, which permits the cal- 
culation of the difference in voltages (from 20 mv to 
2 v) at a discrete moment of time with an error not ex- 
ceeding 1% of the maximum value of the scale, 

V. A. Brik, in a report entitled “A digital computer 
for the composition of a program for the working of parts 
on a milling machine," discussed a device for preparing 
a program for the working of parts whose profiles con- 
sisted of rectilinear segments and circles. On the basis 
of the coordinates of the ends of the segments and the 
centers of the circles,the devices give the values ofthe 
linear increments of the coordinates in a binary code, 


IV. Section on Elements and Devices of 


Automation and Remote Control 











In a report of E. A. Andreeva a method was given 
for calculating the expense and power characteristics of 
a “jet-choke” element for viscous compressible and in- 
compressible liquids. 

L. A. Tenenbaum, in introducing simplified assump- 
tions and hypotheses, derived expressions for the expense 








"and power characteristics of a “jet-choke" element dur- 


ing the irregular laminar flow of a viscous incompress- 
ible liquid. 

T. K. Efremova discussed pneumatic relay assem- 
blies with temppral dependencies (guard relay, impulse 
starter,and impulse generator), In her report there was 
given a timer circuit built with these assemblies. 

In a report of V. S. Matorina entitled “Magnetic out- 
put amplifiers for alternating current electromagnets,” 
the operation of magnetic amplifiers working from the 
signals of magnetic logical elements was examined. Two 
methods of calculation for loads with high and low good- 
ness were offered. 

In M, A. Boyarchenkov's report "Standby magnetic 
amplifiers for direct current with increased efficiency,” 
he analyzed the conditions for obtaining maximal power 
under load. New circuits of three-phase standby ampli- 
fiers and full-wave transformerless standby amplifiers 
were offered, 

In a second report M, A. Boyarchenkov discussed the 
operation of a magnetic amplifier for counter emf. 

N. L. Prokhorov analyzed in his report existing 
memory circuits for magnetic logical elements from the 
point of view of continuity (the memory signal during 
exclusion and following inclusion of the power supply 
current), A new continuous memory circuit was suggest- 
ed, 

A. L. Rozovskii, in his report entitled "A contactless 
code-impulse system of telemeasurement," explained 
the principle of construction and the solving circuits of 
a system built of semiconductor and magnetic elements, 

In a report entitled “Contactless program-calculat- 
ing device for the automation of instruments for a line- 



































casting machine,” N. V. Silaev described the circuit of 
such a device, based on the use of magnetic and semi- 
conductor elements, 

V. S. Serzhers reported on the possibility of con- 
structing, from semiconductor elements, circuits of a 
proportional amplifier, a derivative unit, and an inte- 
grating unit for regulators for general use in industry, A 
typical designed assembly of units was worked out, using 
a printing arrangement, 

B. B. Gogodevskii discussed transitional processes in 
electromagnetic mechanisms and the vibration of con- 
tacts, A method was given for estimating the error in 
calculating the dynamic characteristics of electromag- 
netic mechanisms. 








V. Section on Statistical Methods in Auto- 
mation 

In a communication of L. P, Sysoev a solution was 
given for the problem of estimating the parameters and 
of detecting signals which depend nonlinearly on random 
parameters assumed for the background of gauss-corre- 
lated noise. 

M. Yu, Gadzhiev discussed the games problem 
arising in the case of detection of a useful signal mixed 
with noise during independent change of the carrier fre- 
quency. During definite assumed interference relative 
to the class, formulas have been obtained for optimal 
composite strategy of opponents and evaluation of the 
game. 

In a report of I. L. Paishev equipment with continu- 
ous-discrete functioning for the resolution of random 
functions in a canonical series was examined, Experi- 
mental results were given of the practical construction 
of canonical series for a definite class of random functions. 

E, L. Nappel*baum communicated about the form 
of the optical operator for the discovery of a suitable sig- 
nal on a background of normal noises with random distri- 
bution of dispersions. 

The report of A. I. Teiman was devoted to problems 
linked with the distribution of selections of random func- 
tions, A limiting form of distribution selections was de- 
rived,and a new method offered for the evaluation of the 
distribution of selections of random functions for a given 
level, 

In the communication of ochetkov there was 
an examination of the problem of the construction of 
theoretically and practically optimal linear integral im- 
miscible evaluations of the mathematical expectation of 
the correlational function of stationary random processes 
with continuous time. A comparison was made of the 
theoretically and practically optimal evaluations. 

T. I. Tovstukha discussed the effect of random noise 
on the operation of systems of extreme control of step- 
by-step and sloped types. Working parameters of the sys- 
tem are selected from evaluations of the values of math- 
ematical expectation and dispersions on the yield of the 
object of control. 

















In the report of V. M. Baikovskii the capacity of a 
channel with different-discrete modulation was deter- 
mined during the absence of noise. It was shown that 
potentially this method of transmission has greater ca- 
pacity than delta-modulation and impulse-code modula- 
tion. The connection of the rate of transmission of in- 
formation during a different-discrete modulation ,accord- 
ing to the channel with the fluctuating noises was deter- 
mined from the statistics of the transferable information. 

V._M, Pomazan communicated about a theoretical 
and experimental investigation of time systems of tele- 
measurement with various cycle periods and different 
methods of indication (with memory or integrating out- 
put devices). 





VI. Section on the Theory of Relay Cir- 
cuits and Finite Automata 








In the report of V, D. Kazakoy entitled "The shape 
of minimal forms of symmetrical Boolean functions of 
an arbitrary number of variables,” it was demonstrated 
that the minimal sp- and ps- forms of symmetrical func- 
tions have a completely determined form and on the ba- 
sis of this an account is given of a method of direct deter- 
mination of the minimal forms of arbitrary symmetrical 
functions of n variables. Approximated estimates were 
given of the maximal number of minimal forms there 
are in an arbitrary symmetrical function of n variables, 

In a report of Y, P, Didenko there was an explana- 
tion of a digital method of minimization of Boolean func- 
tions with calculation of unutilizable states, The method 
permits rapid obtaining of partial minimal forms and can 
be distributed on Boolean functions with a large enough 
number of independent variables, An examination was 
also made of the method of constructing bridge structures 
of relay devices on the basis of partial minimal forms, 

V. V.:Vorsheva at the beginning of her report gave 
a short review of existing work on the study of circuits 
with real contacts, By real contacts are understood relay 
contacts in which there can be arbitrary changes of se- 
quence of the interdependence of normally open and 
normally closed contacts of one relay or another. Inthe 
report there are an analysis and a synthesis of the cir- 
cuits operating on real contacts, A method was explain- 
ed for obtaining a partial minimal form giving the cir- 
cuit, when realized, contacts free of competition. 

T. M, Aleksandridj devoted his report to problems 
of the synthesis of switching circuits on the basis of two 
logical operations, the Schaeffer prime and the function 
which is double it, In the report it was proposed to real- 
ize the Schaeffer function and the function which is 
double it in elements operating on the principle of 
switching the current with transistors of type p-n-p and 
n-p-n respectively, 

Examples were given of the construction of logical 
circuits made up of the elements examined, 

In the report of L. A, Gusev, entitled "On minimi- 
zation of the structure of finite automata," a review was 
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given of the basic concepts of the theory of finite autom- 
ata and a study made of the equivalent conversions of 
finite automata and the problem linked to them of mini- 
mizing the number of states of the finite automaton. 

The report of O. P. Kuznetsov was devoted to a de- 





scription of logical networks in which the delay times 
for the various elements are unequal, The author 
showed that the relationships determining the behavior 
of such networks are more complex than the known re- 
lationships occurring in each automaton. In particular, 
there exist autonomous logical networks, whose se- 
quences of states are periodic. 

In the report of Y. D, Kazakov and Y, V, Naum- 
chenko entitled "On the realization of Boolean functions 
of a variables on contactless logical switches by the 
method of redetermination," a new approach to the min- 
imization of relay structures was explained, consisting 
of the fact that a given Boolean function leads at the 
expense of unutilized states to an ordered (in particular, 
symmetrical) form,with the help of simple matrix con- 
versions, easily mechanized, The realization of sym- 
metrical functions can be done on the contactless mod- 
els of neurons of Calberston and Makkalok-Pitts. Such 
a realization permitted a sharp reduction in the number 
of necessary elements in the circuit. 

The report of A. D. Talantsev, entitled "On the use 
of logical- algebraic operators of transition to the analy- 
sis and synthesis of finite automata of a special type," 
was devoted to an examination of the logical means of 
synthesis and analysis of potential-impulse circuits oper- 
ating under asynchronous conditions (on a nonqualitative 
time scale). It was shown that the procedure of inte- 
grating potential- impulse forms permits in some cases a 
reduction in the number of necessary logical elements 
by one-half or two-thirds, 
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Vil.Section on Automated Blectrical Drive 


O. 1, Khasaev, in his report entitled "The opera- 
tion of an induction motor from a frequency converter 
with semiconductor triodes," presented an analytical and 
experimental investigation of fixed and transitional proc- 
esses of an electrical impulse converter in the angle of 
rotation of an impulse operation with a three-stator step- 
by-step motor, Analytic expressions were found for the 
basic parameters of the converter (the angles of overlap 
and lead) and an analysis was given of some conditions 
of its operation. 
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In the report of V. D, Vershinin entitled "The use of 
the principle of invariantness in the stabilization of the 
speed of a direct current motor," there was an explana- 
tion of the conditions of complete invariantness during 
action on the voltage of a signal armature upon the cur- 
rent of the motor, and the possibility of it realization was 
shown, It was shown that the fulfillment of the condi- 
tions of complete invariantness in the actual circuit,with- 
out fail, brings upon itself instability of the system. Com- 
pensation during complete invariantness brings the sys- 
tem to astatism of an order equal to the number of com- 
pensated components (basic and productive) of the com- 
pelled error. 


In the report of O. A. Kossov, entitled "The direct 
current drive with impulse semiconductor rectifiers," 
the economy of impulse regulation of the speed of a di- 
rect-Current motor was proved in the presence of a diode 
shunting armature circuit, an analysis was given of the 
statistical characteristics of the broken system, and the 
effect of feedback was determined, The author pre- 
sented experimental data on the test of a model of the 
drive with a 0.76 kw motor. 


In the report of Chao Chou-lun entitled "Optimal ' 
control of cylindrical flying pruners with eccentrics," 
there was an examination of problems in improving the 
quality of SAR during introduction of feedback; criteria 
of the stability of a broken system were giverjand the 
conditions for formation of autooscillations were found. 


In the report of A. R. Dzhelyalov, entitled "A syn- 
chronous motor with longitudinal- transverse excitation 
as the object of automatic regulation,” the statics of 
regulation during control of excitation on both axes were 
examined in connection with the problem of slewing of 
the rotor during mutual synphasing of electrical drives 
with a reciprocating motor. The author concluded that 
the control of motors is the object of automatic control, 
described the principles of construction of systems of au- 
tomatic turning of the motor's rotor, analyzed the con- 
ditions of rotation, and also gave a qualitative character- 
ization of those conditions. 











M. Yu. Gadzhiev, F. B. Gul'ko, 

A. R. Dzhelyalov, E. E. Dudnikov, 
V. D. Kazakov, I, A. Litovchenko, 
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